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ABSTRACT 


By 

Chester  E.  Canada  and  Jerry  Ward 
DOD  Explosive  Safety  Board 


UPGRADE 

of 

"Structures  to  Resist  the  Effects  of  Accidental  Explosions" 
(Array  TM  5-1300/NAVFAC  P-397/AFM  88-22) 


For  the  first  60  years  of  the  20th  Century,  criteria  and  methods 
based  on  the  results  of  catastrophic  events  were  used  for  the 
design  of  explosive  facilities.  The  criteria  and  methods  did  not 
include  a  detailed  or  reliable  quantitative  basis  for  assessing 
the  degree  of  protection  afforded  by  the  protective  facility.  In 
the  late  1960‘s  quantitative  procedures  were  set  forth  in  the 
first  edition  of  the  present  raanual,  "Structures  to  Resist  the 
Effects  of  Acoidental  Explosions".  The  manual  was  based  on 
extensive  research  and  development  programs  which  permitted  a 
more  reliable  approach  to  current  and  future  design  requirements. 
Since  the  original  publication  of  this  manual,  more  extensive 
testing  and  development  programs  have  taken  place.  This 
additional  research  was  directed  primarily  towards  materials 
other  than  reinforced  concrete  which  was  the  principal 
construction  material  referenced  in  the  initial  version  of  the 
.An  upgrade  to  the  manual,  describing  new  design 
has  become  essential.  Design  methods  in  the  proposed 
upgrade  provide  required  structural  protection. 

This  paper  reviews  differences  and  additions  between  the  earlier 
version  and  the  proposed  upgrade.  The  planned  schedule  for 
technical  review,  Tri-Service  coordination,  and  publication  is 
presented. 
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EXPLOSIVE  FACILITY  DESIGN 
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Enhancanieiits  of  the  Prediction  of  Ground  Shock 
From  Penetrating  Weapons 


James  L.  Drake,  Elizabeth  B.  Smith  and  Scott  E.  Slouin 

Applied  Research  Associates.  Inc. 

South  Royalton,  Vermont 


ABSTRACT 

This  paper  presents  a  procedure  that  predicts 
peak  particle  accelerations,  velocities,  displace¬ 
ments  and  stresses  on  axis  beneath  bombs  exploding 
on  or  within  the  ground.  It  uses  and  reconciles 
the  high  ex.plosive  ground  shock  data  base,  results 
of  numerical  calculations,  and  theoretical  con¬ 
siderations,  including  conservation  laws,  cavity 
expansion  theory  and  similitude.  The  method 
reproduces  most  observed  contained  burst  explosion 
data  in  media  ranging  from  loose  soil  to  hard 
rock. 


BACKQROUND 

This  paper  is  an  update  to  a  paper  entitled 
“Ground  Shock  from  Penetrating  Conventional 
Weapons"  authored  by  Jones  1.  Drake  and  Charles  0. 
tittle,  Jr,  The  original  paper  presented  an 
empirical  method  to  predict  the  ground  shock 
environment  in  soil  from  conventional  weapons  as  a 
function  of  burst  position,  soil  indices  and 
burster  layer  thickness.  The  analysis  presented 
here  describes  certain  enhancements  to  that  pre¬ 
diction  method  especially  the  prediction  of  the 
near  source  ground  shock  environment. 

SOIL  PROPERTY  EFFECTS 

While  no  single  material  property  index  or 
combination  of  indices  can  be  used  to  fully 
prescribe  ground  shock  propagation,  the  seismic 
velocity  cj,  the  primary  loading  wavespeed  Cq  and 
th-j  comractible  air  void  volume  to  are  valuable 
ir.iices  for  assessing  ground  shock  magnitudes. 

The  seismic  velocity  is  associated  with  propags- 
tion  of  low  amplitude  waves  through  the  in  sttu 
eate'-ial  and  is  not  always  indicative  of  the  beha¬ 
vior  of  in  situ  materials  at  higher  stress  levels. 
The  primary  loading  wave  velocity  >s  controTitd  by 
the  in  situ  material  response  at  stress  levels  of 
engineering  interest  and,  as  explained  below,  is 
generally  slower  than  the  seismic  velocity  and  a 
oettcr  indicator  of  the  dynamic  in  situ  response. 
Unfortunately,  the  primary  loading  wave  velocity 
is  considerably  more  difficult  to  determine  than 
is  the  seismic  velocity.  Qrcurtd  shock  attenustion 
rates  with  depch  are  controlled  by  the  compaction 
of  wterial  during  the  passage  of  the  stress  wave. 
Dry  week  rocks  such  as  sandstones  attd  tuff  have 


crushable  air  voids  and  attenuate  the  shock  more 
rapidly  than  hard  rocks  that  contain  few  or  no 
voids.  Alluvial  materials  and  soil  have  high 
crushable  air  voids  and  attenuate  the  shock  more 
rapidly  than  all  but  the  weakest,  most  dry  porous 
rocks . 

The  seismic  velocity  and  primary  loading  wave 
velocity  provide  a  fundamental  relation  between 
space  and  time  lar  scaling  ground  motions.  Since 
the  characteristic  time  for  a  given  event  is 
inversely  proportional  to  the  propagation  velo¬ 
city,  explosions  in  strong  competent  materials 
(high  seismic  and  loading  wave  velocities)  will 
produce  much  shorter  duration  ground  motion  pulses 
than  like  bursts  in  soft  rock  and  soils.  The  pri¬ 
mary  loading  wavt  velocity  provides  a  measure  of 
the  modulus  of  the  medium.  Thus,  explosions  in 
strong  competent  materials  with  high  propagation 
velocities  will  produce  higher  accelerations  and 
lower  displacements  than  corresponding  explosions 
in  soft  madia  having  lower  propagation  velocities. 

OROUkO  SHOCK  ENVIRONMENT 

Approximate  analytical  solutions  for  an  ex¬ 
panding  sphtrical  shock  wave  in  a  nearly  incom¬ 
pressible  media,  were  used  as  the  basis  of  the 
near  source  region  analysis.  Considering  conser¬ 
vation  of  mass  between  the  expanding  cavity  and 
the  Observation  point,  it  was  found  that  the  flow 
field  can  be  estimated  by 


where  v  is  the  particle  velocity,  Pglt)  is  the 
expanding  cavity  radius  and  n  is  a  constant  of 
approximately  2.  The  magnitude  of  the  particle 
velocity  and  the  attenuation  tn  the  high  pressure 
large  flow  region  near  the  source  are  controlled 
by  the  cavity  expansion  mie  and  tha  geometric 
spreading  of  the  shock  front. 

The  kinetic  energy  in  the  source  region  can 
be  estimated  by 

J  ^  <ir  •  (l  (2) 
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assuKiog  n  =  2  and  Pq  is  the  initial  density. 
Assuming  that  within  the  source  region  the  kinetic 
energy  in  the  flow  field  is  approximately  equal  to 
1/2  of  the  total  weapon  yield  and  also  noting  that 
for  TNT  1kg  =  4.6186E6  joules,  then 

u 

Vjj  »  606.2  (-^)  (3) 


where  Vp  is  the  peak  particle  velocity  and  w  is 
the  contained  yield  of  the  weapon  in  kg. 

For  incompressible  flow  n  =  2.  However,  it 
has  been  shown  that  n  can  oe  generalized  to 
account  for  volume  changes  in  compactible  mat¬ 
erials  by  the  expression 


where  e  is  a  constant  <<1,  that  relates  the  volume 
change  to  the  shear  strain  in  the  flow  field. 
Attenuation  rates  observed  in  the  data  base  show 
that  €  a  Co  where  Cp  is  the  compactible  air  voids. 
For  saturated  material •  n  is  less  than  2  and 
depends  on  the  strength  of  the  material. 

Empirical  fits  were  made  to  peak  ground 
Mtion  and  stress  data  obtained  from  contained  he 
detonations  in  soil  and  to  results  of.  finite  dif¬ 
ference  calculations  in  various  soil  and  rock 
geologies.  The  resulting  expressions  represent 
the  best  estieate  of  the  contained  explosion  data 


base  and  are  consistent  with 
and  near  field  analysis. 
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Terms  in  the  above  equations  are  defined  as 

ap  =  peak  radial  acceleration  (g) 
dp  =  peak  displacement  (m) 

<Tp  =  peak  stress  (Pg)  (1  Kbar  =  10^  pg) 
r  =  radial  distance  from  the  explosion  (m) 
g  =  gravitational  acceleration  (m/s^) 
n  a  peak  velocity  attenuation  exponent 
S  a  equation  of  state  factor  (si. 5  for 
geologic  media) 
t^  =  rise  time  (s) 

The  peak  velocity  and  stress  attenuations  from  the 
above  equations  and  the  soil  parameters  from  Table 
1  are  compared  to  the  data  in  Figure  1  for  both 
dry  and  saturated  media.  The  previous  fits  to  the 
data  are  also  included  in  the  plots.  The  new 
methodology  extends  the  prediction  into  much 
higher  stress  regimes. 

The  peak  particle  velocity  in  the  source 
region  are  included  in  the  above  equations.  The 
magnitude  and  attenuation  rate  are  controlled  by 
the  cavity  expansion  rate  and  the  geometric 
spreading  of  the  shock  front.  This  region  extends 
to  a  range  of  .155  m/kgl/3  which  corresponds  to 
the  approximate  size  of  the  expanded  cavity  in 
most  rocks.  Within  this  region,  the  material 
displacements  are  largu  with  ^-espect  to  the  ini¬ 
tial  position,  so  that  material  is  pushed  into  a 
relatively  thin  shell  with  an  approximately 
constant  kinetic  energy  density.  The  resulting 
peak  particle  velocity  is  only  weakly  dependent  on 
the  material  properties. 

The  loading  wave  velocity  is  a  function  of 
the  loading  intensity.  It  generally  starts  from 
the  primary  loading  wavesoaed  and  monotonically 
inereasaa  with  increasing  stress.  For  most  geo¬ 
logic  materials,  the  propagation  speed  can  be 
cpproximated  by  relating  it  to  the  initial  loading 
wavespead  and  the  peak  particle  velocity  through 
the  use  of  the  factor  S  which  is  a  function  of  the 
overall  compressibility  of  the  material. 

The  material  model  for  all  rocks  and  soils 
has  been  simplified  to  relate  the  loading 
wavespead  to  the  peak  particle  velocity  with  the 
parameter  S  a  t.S.  This  has  been  done  to  allow 
the  user  to  employ  this  methodology  with  relati¬ 
vely  little  knowledge  of  the  actual  material  model 
of  the  geology  of  interest,  if  the  actual 
uniaxial  stress  strain  response  is  known,  it  can 
easily  be  substituted  to  relate  the  loading 
wavespaed  to  the  peak  particle  velocity. 

The  magnitude  of  the  stress  and  ground  motion 
will  be  greatly  enhanced  as  the  weapon  penetrates 
sore  deeply  into  the  soil.  The  concept  of  a 
coupling  factor  was  introduced  in  the  previous 
paper  and  is  summarized  here.  This  factor 
accounts  for  the  effect  of  weapon  penetration  on 
tfie  ground  shock  parameters  and  is  defined  as  the 
ratio  of  the  ground  shock  magnitude  from  a  par¬ 
tially  burieii  weapon  to  the  ground  shock  from  a 
fully  contained  burst  in  the  same  media.  The 
coupling  factor,  f.  can  be  determined  from  Figure 
i  and  the  acceleration,  velocity,  stress  and 
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Flgup#  lb.  Comparison  of  paak  vmloclty  and  stress  attenuation  prediction  to  explosive  data. 
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dlsplacemont  obtained  from  the  above  equations 
should  be  multiplied  by  this  factor  for  shallow 
buried  bursts. 


Scaled  Depth  of  Bunt  (m/kg‘^’ ) 
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Figure  3.  Particle  velocity  time  history. 


CONCLUSIONS 


Figure  2.  Ground  shock  coupling  factor  as 
a  function  of  scaled  depth  of 
burst. 

Stress  and  notion  time  histories  can  be 
characterized  by  exponential  tine  histories  that 
decay  rapidly  in  aaplitude  and  broaden  as  they 
propagate  outward  fron  the  explosion  (Figure  3). 
The  characteristic  time  for  these  time  histories 
can  be  measured  in  arrival  tine  fron  the  source, 
t*  ■  rh\,  the  tine  of  peak,  tp  •  t/cl,  and  the 
tine  of  positive  phase  duration  t«  where 


U  ■  2*T2  &  -  0.36  t>  (13) 

^P 

Assuming  a  linear  rise  to  peak  velocity,  the  velo¬ 
city  waveform  can  be  approximated  by 


Empirical  expressions  were  derived  from  a  fit 
to  a  large  body  of  ground  shock  data  from  buried 
bursts  in  soil.  These  expressions  were  derived 
using  a  more  sophisticated  material  modal  than 
originally  used  by  Drake  and  Little.  The 
resulting  equations  more  accurately  predict  the 
near  source  environments.  These  n^ifications  of 
Drake  and  Little's  work  do  not  negate  the  previous 
study.  The  extend  it  over  a  broader  range  of  peak 
stresses  and  over  a  broader  range  of  geologies. 


The  method  presented  here  is  not  necessari ly 
rocoaaended  over  the  previous  method  but  is  more 
valid  at  regions  closer  to  the  source  and  is  based 
on  a  aort  rigorous  theoretical  approach  which  clo¬ 
sely  approximates  actual  test  data  in  virtually 
all  types  of  geologic  materials.  The  two  methods 
are  in  close  agreeaent  at  farther  ranges. 


v(t)  »  Vp(l.^)  t,  <  t  «  tp  (14) 

-(*  -  tp  j 

v(t)  •  Vp(*l-^J-)e  'P  t  <  tp  (15) 


where  t(j  is  the  absolute  time  at  the  end  of  the 
positive  phase  duration.  Since  the  time  domain 
waveform  features  are  inversely  proportional  to 
the  propagation  velocity,  explosions  in  competent 
aaterial  (high  o^)  will  produce  shorter  pulses 
with  higher  ncceltrstions  and  lower  displeccments 
than  corresponding  explosions  in  less  competent 
material  (low  c^).  Consistent  scaling  rela¬ 
tionships  also  rtquire  higher  aocalerations  to  be 
easooiated  with  hightr  atressas  and  lower  disple- 
cemants.  Conversely,  low  accelerations  will  be 
aasociated  with  lower  stresses  end  higher  disple- 
ceoMnts. 
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SOIL  CATEGORY  TYPES  AMO  EXAMPLES 


0«ns*  Dry  Soil  alluvlua 

dans*  sand 


Htdlu*  Oanslty  sand,  loaa,  alluvlu* 
Dry  Soil 


TYPICAL 
IN  SITU 
DENSITY 
Po 

(kg/»3) 


TYPICAL  TYPICAL 

SEISMIC  INITIAL  VELOCITY 

WAVESPEED  LOADING  ATTENUATION 
Co  EXPONENT 

(«/s)  (m/s)  n 


2.5  to  3.0 


LoM  0«ns1ty 
Dry  Soil 


looso  sand,  loam 


Saturatad  Soil  all  typas 


.  . .  — .1 

High  Strength 

Rock 

quartxlta,  diabase, 
baaalt,  granit* 

2650 

— 

5500 

4500 

2 

Nadlua  Strength 
Rook 

shales,  parous  send* 
stones  and  llasstonaa, 
ohlorlta 

2500 

3400 

2750 

2 

LoM  Strength 

Rock 

porous  tufP.  clay, 
shales,  schist 

2000 

2050 

1550 

2 

vary  Lom 

Strength  Rock 

vary  porous  and  tri¬ 
able  Msathered  rock 

1400 

1400 

1050 

2 

I 


Taol*  i,  Ganaric  soil  and  rock  proparties. 
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ABSTRACT; 

This  paper  presents  the  results  of  current  research 
aimed  at  developing  an  unproved  single*degree-of-freedom 
(SDOF)  analysis  model  for  buried  arches.  Several 
refinements  in  the  SDOF  analysis  procedures  for  buried 
arches  are  derived,  including  (1)  development  of  a  general 
soil-arch  interaction  model  for  defining  the  structure-media- 
interaction  (SMI)  loads  on  the  arch,  (2)  development  of 
equivalent  SDOF  response  parameters  for  arches  of  arbitrary 
half-central  angle,  and  (3)  development  of  a  resistance 
function  that  accounts  for  the  moment-thrust  load  path  the 
arch  follows  during  deformation.  The  work  assumes  a 
cylindrical  arch  cross-section  of  arbitrary  arch  angle  subjected 
to  a  vertically  propagating  planar  ground  shock.  The  elastic 
and  plastic  deformed  shapes  are  derived  for  a  characteristic 
load  distribution  and  used  to  define  the  equivalent  parameters 
for  the  SDOF  model.  This  characteristic  load  is  replaced  in 
the  response  calculation  by  a  simplified  soil-arch  interaction 
model  which  more  accurately  reflects  the  variation  of  the 
interface  loads  on  the  arch  as  it  is  engulfed.  While  the 
resent  model  assumes  a  planar  ground  shock  from  a  nuclear 
urst,  the  methodology  is  ^ptablo  to  conventional 
munitions. 

INTRODUCTION; 

The  buried  arch  is  a  popular  geometry  for  protective 
structure  design:  its  curved  geometry  is  capable  of  resisting 
high  pressures  through  a  combination  of  bending  and  iltrust 
while  providing  a  large  interior  span  that  can  be  used  to 
store  and  protect  military  equipment  such  as  aircrafL 
However,  because  of  its  geometry,  analysis  of  the  ground 
shock  loads  acting  on  the  arch  and  its  dynamic  response  is 
extremely  complex  and  not  well  understood.  As  a  result, 
current  methods  for  analysis  and  design  of  buried  arches  are 
generally  crude,  over  simplistic,  and  inadequate  for  predicting 
Ste  dynamic  toads  tuid  l^ge  deformation  response  of  these 
structures. 

One  approach  to  modeling  the  complex  response  of 
buried  arches  is  the  finite  element  method  wherein  the 
complete  arch  geometry  and  structure-medla-interaciion  can 
be  modeled  explicitly.  The  fast  computational  speeds  and 
large  memory  c^utciiy  of  modem  computers  allow  one  to 


model  the  arch  in  great  detail  without  too  great  a  penalty  in 
tum-around  time.  However,  finite  element  models  are  still 
time  consuming  to  prepare  and  the  results  ate  often  difficult 
ia  interpret.  In  addition,  finite  element  models  are 
cumbersome  to  use  in  preliminary  design  calculations, 
survivability/vulnerability  assessments,  and  cost-trade 
studies.  Hence,  there  is  sdll  a  need  for  more  simplified 
analysis  procedures  that  can  accurately  predict  the  loads  and 
response  of  buried  structures.  Furthermore,  the  high  cpu 
speeds  of  modem  computers  alleviate  the  need  to  make  many 
of  the  simplifying  assumptions  in  SDOF  models  so  that  tlm 
accuracy  of  these  models  can  be  improved  and  potentially 
^proach  that  of  more  sc^ihisticated  finite  element  models. 

The  research  reported  herein  is  addressing  this  need 
through  the  development  of  refined  SDOF  analysis 
procedures  for  buried  arches.  This  work  attempts  to  include 
as  much  of  the  appropriate  physics  as  is  possible  while 
retaining  the  simplicity  of  the  SDOF  modeling  approach. 
While  the  work  is  currently  ongoing,  the  model  derivation  is 
complete  and  is  presented  in  this  paper.  Comparisons  with 
test  diua  will  be  presented  in  a  ftiture  paper. 

SDOF  MODEL  FORMULATION: 

A  schematic  of  the  loading  and  response  of  a  buried 
arch  os  it  is  engulfed  by  a  ground  shock  wave  is  shown  in 
Figure  1.  As  the  shock  wave  impinges  on  the  arch,  the 
ground  shock  stresses  ore  reflected  due  to  the  impedance 
mismatch  between  the  soil  and  die  structure.  Inidally,  only 
asmall  portion  of  the  arch  perimeter  is  loaded  by  the  ground 
shock.  As  the  ground  shock  propagates  deeper,  the 
perimeter  of  the  arch  that  is  loaded  increases,  until  at 
enguliment  the  entire  arch  perimeter  is  loaded.  Response  of 
the  arch  is  primarily  in  the  first  symmetric  bending  mode 
and  analysis  of  recent  test  data  and  results  of  finite  element 
calculations  tl]  show  that  the  peak  flexural  response  occurs 
during  this  engulfment  ptiase.  Those  studies  have  also 
shown  that  the  response  of  the  arch  is  not  sensitive  to  the 
details  of  the  loading  and  that  the  primary  loading  on  the 
arch  is  due  to  the  normal  pressures  on  the  arch.  Inierfhce 
shear  stresses  developed  on  the  arch  interface  ore  generally 
small  and  can  be  neglected  [1]. 
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Figure  1;  Planar  Ground  Shock  on  a  Buried  Arch 

Crawford,  et  al.  [2],  have  observed  that  the  ground 
shock  loads  on  buried  cylinders  and  arches  can  be  described 
using  a  characteristic  q + pcos20  stress  distribution.  This 
load  distribution,  shown  in  Figure  1,  consists  of  a  uniform 
component  given  by  ^  and  a  nonuniform  component  given 
given  by  pcos2^.  However,  the  magnitude  and  distribution 
of  the  interface  stresses  on  the  arch  vary  as  the  arch  is 
engulfed  so  that  the  characteristic  load  distribution  is  not 
constant  in  time,  rather,  the  relative  magnitude  of  the  two 
components  varies.  However,  as  will  be  demonstrated  in  the 
derivation  of  the  deformed  shapes,  the  q  term  does  not 
influence  the  deformed  shtqre  or  capacity  (except  through  the 
moment'thnist  interaction  diagram)  so  that  the  temporal 
variation  of  the  relative  mapiui^  of  the  two  components  is 
not  important  to  the  derivation  of  the  SDOP  model 
parameters.  It  is  Important  in  determining  the  moment 
capacity  and,  consequently,  the  arch  dynamic  re^Kinse. 

Based  on  these  observations,  on  equivalent  SDOF 
model  of  the  arch  con  be  developed  os  illustrated  in  Figure  2. 
The  interface  loading  on  the  arch  is  represented  by  a 
characteristic  q  +pcQs2^  load  distribution.  The  response  is 
assumed  to  be  in  the  first  symmetric  bending  mode  and  the 
assumed  mode  shape  is  utk^  as  the  deformed  sltape  of  the 
arch  under  this  load  distribution  applied  statically.  The 
equation  of  nuMkm  for  the  equivalent  SDOF  system  is  given 
by 


where  Mt,  Rg,  and  Fg  are  the  equivalent  moss,  resistance, 
and  oppUed  force  for  the  SDOF  system  and  is  the 
vertical  acceleration  of  the  crown.  The  derivation  of  these 
parameters  and  the  ioadiag  function  are  discussed  in  the 
following  paragraphs. 

Loading  Function. 

White  the  characteristic  q-(-pcos2p  load  distribution 
is  used  to  derive  the  deformed  sltape  and  equivalent  SDOF 
parameters,  a  more  accurate  description  of  die  loading  is 
desired  in  solving  for  the  te^ixmse  of  the  arch.  Theappro^ 


followed  uses  a  simplified  SMI  model  based  on  linear  wave 
theory  to  defme  the  interface  stresses  on  the  arch.  This 
^proach  has  been  used  with  good  success  in  SDOF  models 
for  buried  slabs  subjected  to  ground  shock  from  conventional 
and  nuclear  munitions  P].  The  interface  stresses  on  the  arch 
are  determined  using  a  combination  of  simple  wave 
propagation  and  rigid  body  mechanisms  to  pose  the 
boundary  conditions  between  die  soil  and  the  structure.  The 
rate  of  external  work  done  by  the  interface  stresses  is  then 
calculated  at  each  time  increment  and  used  to  determine  the 
response  of  the  SDOF  system  and  the  equivalent  q  and 
pcos2p  interface  stress  components. 


Figure  2,  Arch  SDOF  Model 

Assuming  a  fulUslip  interface  condition  (zero 
interface  strength),  the  interface  load  consists  of  normal 
stresses  only  which  are  given  by  (see  Figure  3) 

cfii«crjf,+pCt(v>,-ii#icosp-w) 

where  *  fnefleld  ground  shock  stress  normal  to  the  arch 
interface.  pCi  a  soil  impedance,  Vff^  >  freofleld  velocity 
normal  to  the  arch  interface,  iingm  vertical  rigid  body 
velocity  of  the  arch,  p  ■  angle  to  point  of  Interest  on  the 
arch,  and  w  R  radial  velocity  of  the  arch  at  the  point  of 
Interest  A  planar,  vertically  propagating  ground  shock  is 
assumed  and  the  Creefield  ground  shock  velocity  and  stress 
are  assumed  constant  behind  the  shock  front  The  component 
of  the  freefleld  velocity  normal  to  the  interioco  is  given  by 

Vff^  m  VffCOS  p  (3) 

The  freefleld  stresses  normal  to  the  arch  interface 
ore  determined  by  transforming  the  freefleld  vertical  and 
horizontal  stresses  to  the  plane  nonnol  to  the  arch  Interface. 
Representing  the  horizontal  stress  by  of/i^  >  Kq  off^  and 

assuming  the  vertical  and  horizontal  freefleld  stresses 
represent  principal  stresses,  the  freefleld  stress  normal  to  the 
or^  interface  IS  given  by 

(4) 
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(8) 


PfS, 


W^, 

wco^^rd(t) 


Vff-liu 


Figures.  Planar  Wave  Loading  Model 


and  the  equivalent  uniform  component  (evaluated  at  the 
crown)  is  given  by 


q«=o),(0)-Pe 


(9) 


Deformed  Shape 
Elastic  Response 

The  mode  shape  for  Lhe  arch  in  the  elastic  response 
range  is  taken  as  the  deformed  sliape  for  the  arch  undtf  the 
chmacterisdc  q + pcos2^  load  distribution.  The  solution  for 
the  deformed  shape  follows  the  general  approach  outlined  by 
Timoshenko  and  Gere  in  Reference  [4]  and  has  been  derived 
for  Qxed  and  pinned  (2-lunged)  support  condititms.  Referring 
to  Figure  4.  the  elastic  force  r^tants  at  any  arbitrary  angle 
from  the  crown  are  given  by 


where 

qfm  and 


V a  ^ {{a vjsin  0  •  6  sin  20) 

(10) 

Af  a  ^  ({ A  V  )cos0  +  3  cos  20)  - 

(11) 

The  rate  of  external  work  done  by  the  interface  stresses  in 
deforming  the  arch  is  given  by 


A/ •«>  3  cos2^  •  {Av}  cos  0) 


(12) 


-♦* 


uicos20rd0 


where  the  bracketed  terms  (Am)  and  (Ay)  depend  on  the 
(Q  support  conditions  and  ate  of  the  general  form 


The  constants  for  Equations  9  •  12  axe; 


f 

’*1 

*pCl  ' 

[Vff^UKt) 

wcos0rd0- 

*» 

tt»*rd0 

H 

where  the  limits  of  integration  are  0r  ■>  arch  angle  to 
complete  unloading  (tTn  si  0)  and  05  »  arch  angle  to  the 
cunem  position  of  the  shock  fraiu. 

This  external  work  rate  is  then  used  to  deTmo  the 
equivalent  load.  where 

m  tm.  (7) 

Wo 

Likewise,  the  equivalent  nonunifonn  loading  componont  is 
given  by 


Av  »  (cos  309  •  cos  0)  ) 

^  Ipinned:  •?  sin  30a  9  sin0o 

^  Irwed:  -2  sin  300  •  6  sin  09 
‘  Ipinned:  6  (cos  30)  -f  cos0} ) 

Av}**  /fixed;  0.0 
Ipinned'  0.0 

^  ,  /fixed:  025  (cos  40)  +  8  cos20)  -  9) 
Ipinned:  sin  4^  sin  20) 

/fixed'  3  sin  209 
Ipinned:  -6005  20) 

Aiti"  /fixed:  0.0 
Ipinned*  0.0 
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nl.Bxed:  2-2cos2^ 
^  Ipinned:  -3  sin  2^ 

„  (fixed:  -  sin  2^ 

^  Ipinned:  2  (2  +  cos  2^) 


aid 


where  =  deflecdon  at  the  crown  (^=  0). 


(16) 


»  (fixed:  -2.0 
^  Ipinned:  0.0 


a.  eUSTtC  FORCE  RESULTAWTS 


t>.  STATIC  C0I1APS8  MODE 


Figure  4.  Arch  Parameter 
Definitions  and  Static  Collapse  Mode 


Finally,  I 


w  ■  ^  {{A«)  cos  4  +  cos  2#  •  0.5  (Av)  4  Sin^  +  {An]] 


Note  that  for  both  the  pinned  and  fixed  arch  the 
the  work  done  by  the  uniform  load  component  is  zero.  This 
results  in  an  inconsistency  in  the  case  of  a  semi-circular 
arch,  where  ^:f9iy*aad  sin  2^  =  0.  In  this  case,  the  load 
component  performing  work  is  normalized  to  the  load 
component  not  doing  work  and  the  load  transformation 
factor  does  not  have  a  physical  meaning.  Due  to  this 
inconsistency  the  more  general  statement  of  the  equation  of 
motion  has  been  retained. 

Plastic  Response 

The  limit  capacity,  deformed  shape,  and  plastic 
response  parameters  for  the  arch  are  derived  using  the  limit 
plasticity  assuming  the  five-hinge  static  collapse  mode  as 
shown  Figure  4.  From  kinematic  considerations,  the 
relationships  between  the  vertical  velocity  at  the  crown  and 
the  rotation  of  the  arch  segments  1  and  2  are  given  by 


»  • 

where  A  and  ft  are  the  rotation  rates  of  segments  1  and  2, 
reflectively;  and  ^0  is  the  angle  to  the  haunch  hinge. 


what 

Awo*  5(0i  sin  4^0  •  sin 

"  ipitmed:-Scos4^-»-8cos2^-3 

fixed:  0.3  (-5  cos  4^  +  4  cos  2^ -f  l) 
pinned:  0JS(>11  sin  4^  -  2  sin  2^ ) 

.  (fixed:  -  sin  4^0  +  4  sin  2^ 

‘  Ipinned:  3  (cos  4^  -  2  cos  2do  1 ) 

This  deformed  shape  is  used  as  the  mode  shape  for 
the  equivalent  SDOF  and  can  be  used  to  derive  the 
transformation  factors.  K14  and  Ki,  resulting  in 

(*)  ftO 

wd^  ■¥  p  w  cos 

•  - Jo - ^  (15) 

Hio(q^-fO.Spsin2#) 


The  capacity,  p*,  of  the  arch  is  derived  by  equating 
the  internal  and  extomal  work  rates  (Note:  the  wo^  done  by 
the  q  componeuc  is  zero  os  in  the  ela^  case)  giving 

p-« _ - \ 

2r*(l-co8^)|(^o  +  CicosqiD  +  Cicos*t(to/  '  * 
where 

(fixed:  2(l-cos^) 

Ipinned:  1  -2003^  Coa-costto 

A, -(fixed:  0.0  Ct-Ucos<o 

(pinned:  1.0  Cza-l 

and  fin  ■  plastic  moment  capacity  of  arch  cross-section. 
Minimualioa  with  respea  to  ^  gives 
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pinned;  cos^  =(V2-l}  +  (2-V2lcos^ 

fixed:  cos^  =-i-{l  +  cos^) 

2 

This  solution  is  a  upper  bound  in  that  the  internal 
work  done  by  the  membrane  stresses  is  neglected.  The  effect 
of  the  existence  of  membrane  forces  is,  however,  accounted 
for  in  determining  the  plastic  moment  capacity,  Mp.ns 
discussed  in  tt»  foUowing  section. 


Resistance  Function 


ADAPTATION  TO  CONVENTIONAL 
MUNITIONS 

The  ixoposed  SDOF  methodology  can  be  adapted  to 
conventional  munitions  for  the  case  of  an  overhead  burst  and 
two-dimensional  response  (for  example,  a  longitudinally 
segmented  arch).  This  is  accomplished  by  replacing  the 
planar  ground  shock  with  a  spherically  propagating  shock 
wave  as  illustrated  in  Figure  6.  Assuming  that  the  radial 
stress  and  velocity  in  the  freefield  are  constant  behind  the 
shock  front,  the  normal  freefield  ground  shock  stress  and 
velocity  are  given  by 


The  resistance  of  the  arch  is  modeled  using  an 
elastic-perfectly  plastic  resistance  function  as  shown  in 
Figure  S,  where  the  stiffness  and  maximum  resistance  are 
determined  using  the  solutions  presented  earlier.  The  major 
refinement  here  is  that  the  maximum  resistance  is  allowed  to 
vary  depending  on  the  moment-thrust  combination  existing 
in  the  arch  at  any  point  in  time.  Observations  from  test  data 
[S,6]  hawe  shown  that  the  response  of  the  arch  is  generally 
in  the  compression  region  of  the  moment-thrust  diagram  and 
that  the  eccentricity  (ratio  of  moment  to  thrust)  decreases  as 
the  arch  is  engulf^  by  the  ground  shock.  This  response  is 
illustrated  in  Figure  S  and  is  modeled  by  allowing  the 
moment  capacity  of  the  arch  to  vary  based  on  the  weighted 
average  of  the  elastic  moment  and  thmst  at  the  arch  critical 
cross-sections  (crown,  0x)  and  ^},  as  calculated  using  the 
equivalent  and  pe  components.  Since  as  the  arch  is 
engulfed  the  magnitude  of  the  uniform  component,  qg, 
increases  relative  to  the  nonunifbrro  component,  pe,  this 
will  naturally  increase  the  magnitude  of  the  thrust  relative  to 
(he  moment,  replicating  test  data  observations.  This  can 
potentially  result  in  the  arch  yielding  at  a  moment 
magnitude  less  than  the  peak  moment  due  to  the  decreasing 
moment  capacity  above  the  balance  point 


RESISTANCE  FUNCTION 


MOMENT-THRUST  DIAGRAM 


q)5^rt  =  '7^+Pjycos2(0  +  j3) 

=  cos(^  +  ^) 

so  that  the  external  work  rate  is  given  by 


WtM^qjj 


*s 

wcos2(p4'^)rdp 


J*M 


(21) 


(22) 


•itu 


wcos^rdf^  •  w/dp 
<  h 


Since  ihe  angles  ^  and  are  related  (as  determined  by  the 
standoff  distance),  ^  can  be  solved  for  in  terms  of  0  giving 


(23) 


tnl 


(24) 


Simitariy, 


cos 


'  2  cos  2#  ' 

\  (  2r  sin  #  sin  2e 

1 

(25) 


Substitution  of  Eqns  24  and  £q.2S  into  Eqn  22  for  the 
external  work  rale  results  in  a  radical  integral  considerably 
more  complex  than  those  for  the  nuclear  case.  Nonetheless, 
closed  form  solutions  for  these  integrals  exist  so  that  a 
computer  efficient  model  can  be  deveUgred. 


Figures.  Resistance  Function. 


The  loading  for  a  conventional  weapon  is  more 
localized  than  that  for  a  nuclear  weapon  so  that  the  deformed 
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shape  and  collapse  mechanism  may  also  require 
modiHcation.  This  can  be  modeled  by  assuming  a 
concentrated  load  at  the  crown  for  the  characteristic  load 
distribution.  The  deformed  shape  and  plastic  limit  capacity 
can  then  be  derived  for  this  chaiacterisdc  loading. 


STATUS  AND  FUTURE  EFFORT 

The  research  is  cuirendy  in  the  final  stages  of  the 
model  development  and  programming.  In  addition  to  the 
refinements  reported  herein,  models  have  been  derived  for 
including  the  effects  of  rise  tutie  in  the  loading  funcdon  and 
for  calculating  the  rigid  body  response  of  the  arch.  The 
elastic  response  is  currently  being  expanded  to  include  a 
three  hinge  turch  (crown,  supports).  After  programming  and 
clieckout  of  the  model,  comparisons  of  the  model  against 
mst  data  from  the  AFWL  idichina  test  series  [7]  and  the 
DNA/VVES  Dynamic  Arch  Test  [8].  Comparisons  against 
finite  element  calculations  [1]  and  other  arch  SDOF  models 
[5.9]  will  also  be  conducted. 


Figured.  Adaptation  to Convcutionol  Weapon 
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ABSTRACT 

A  study  has  been  inducted  to  investigue  expedient  methods 
for  hardening  existing  smictures  and  expos^  assets  against 
the  effects  of  conventional  weapons.  The  expedient 
hardening  methods  covered  in  this  paper  are  soil  benos, 
sandbag^g,  sand  grids,  concrete  modular  revetments,  bin 
revetments,  and  sadrincial  panels.  For  each  of  the  methods 
the  paper  gives  a  short  description  of  the  method  including 
thrtiat  protection  afforded,  a  summary  of  available  test 
resutis,  limitations  of  the  method  and  spwial  consideratioas 
requited  for  deployment. 

1.  INTRODUCTION 

In  peacetime  or  wartime  situations,  needs  can  arise  to 
upgrade  the  hardness  of  an  existing  facility  in  an  expedient 
Cii^ioit.  'this  may  be  due  to  a  change  in  funcdon  of  the 
facing  (Ovg.,  retooiion  of  imptatant  assets  to  the  faciUty), 
an  increase  in  die  perceived  Impottaoco  of  the  facility,  or  an 
tnereafe  in  (he  mtpected  threat.  A  study  has  been  conducted 
to  invesdiate  bom  traditional  and  newly  develr^  methods 
of  expedietu  hardening  and  to  develop  guidance  for  designers, 
the  expedicm  harde^  methods  coveted  in  this  paper  are: 
soil  berms,  srmdbagging.  sand  grids,  concrete  modular 
revetments  (aircrafl,  Biibwg,  StFd)N),  bin  rtvetmen^,  and 
sacrindal  panels.  For  ^uh  of  the  rnethods  investigated,  the 
paper  presents  a  short  description  of  the  method  Inclo^ 
threat  protection  afforded,  a  Siimmary  of  available  test 
results.  Umitatrorae  of  the  method  and  special  oonsideiaiions 
resulted  for  deploymeni.  the  comp^  results  of  (his  study 
can  be  found  in  the  expediem  hardening  adtienduffl  to 
Force  Manual  for  the  Design  of  Fit^Uve  Structures  for 
Co^'*^..lionai  Weapons  Effects  (ARA,  19SS).  The 
addendum  includes  methods  not  covered  herein,  design 
guidioce.  and  a  sekctkM  j(oadm.«;i  to  aid  in  making  optimal 
selections. 

2.  SOIL  BERMS 

Cencral.  Bettns  are  employed  as  free  standing  structures  or 
constructed  against  walls.  A  benned  wall  catftguiauoo  is 
shown  in  Figure  1. 


Threat  Protection.  Berms  provide  protection  against  near 
miss  general  purpose  bombs,  high  explosive  rounds,  and 
ballistic  penetration  (USOCD,  1941),  As  free  staiiding 
walls,  berms  can  be  to  deny  a  dinset  line  of  sight  to  a 
protected  asset  or  a  vulnenble  area  such  as  a  door  opening. 

Recou  tests  on  berms  include  the  CHEBS  series  conducted 
at  Kinland  AFB  (Hyndman,  1987)  and  the  NATO  tesu 
conducted  at  TyodaU  AFB  (Hyde,  1989),  The  CHEBS  tests 
included  free  standing,  full  scale  reinforced  concrete 
revetments  that  were  both  bermed  and  unbermed.  The 
environment  was  provided  by  a  nose  tangent  MARK  83 
general  purpose  bomb.  suUcally  detonated  at  a  standoff  of 
SO  feet  The  NATO  test  examined  the  effectiveness  of  a 
bermed  semi-hardened  reinforced  concrete  wall  subjected  to 
the  NATO  OP  bomb  threat  In  both  tests,  no  dan^  was 
reported  to  the  bermed  structures,  while  the  uubeitsed 
structures  suHenM  sigttincant  front  face  crateHng  and 
backface  spall.  Significant  aiibbut  protection  was  also 
demonstrated  in  these  tests.  The  CHEBS  tests  ^wed  i 
(cducdon  of  idwut  93%  fo  pok  waU  prtsssiues  (at  the  beim- 
wall  intetfsce);  the  NATO  test  showed  similar  reductions  in 
peak  wall  pevtsmi.  Because  of  the  reduced  wall  pressures, 
dte  risk  of  ^  and  breach  are  ipeaUy  reduced  (even  ihou|^ 
total  impulse  delivered  remains  relatively  unchanged).  It 
must  be  cautioned,  however,  that  if  the  weapon  penetrates 
the  berm  prior  to  detonation,  the  coupling  between  the 
explosion  and  the  surrounding  soil  can  result  in  wall 
pressures  gieater  than  for  an  unbermed  walL  Also,  in  a  Bee 
standing  conriguration,  »  berm  provides  little  airblast 
protection.  This  U  because  me  blast  wave  reforms  behind  it 
alter  passfog  over  tlte  top.  this  was  dernonsttated  by  studies 
made  during  the  CHEBS  tests.  The  tests  demonstrate  that 
berms  provide  exceltem  second  or  taullipfo  attack  protection. 

Limitationa/Spfcial  Considerations.  The  main 
disadvantage  of  free  standing  berms  ud  to  a  lesser  extent 
termed  walte  is  theu  targe  space  requirements.  Bernung  may 
not  be  a  practical  hardening  option  for  stiuctures  sited  in 
very  rocky  terrains  or  where  grading  equipment  is  not 
available.  At  air  base  facilities,  berois  sited  nesr  taxiways 
and  runways  may  exacerbate  problems  related  to  blowing 
dust  and  foreign  debris.  Erosion  control  measures  are 
particularly  tmpottani  under  these  circumstancei  Typical 
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facings  used  to  control  erosion  include  sod,  sandbags,  and 
asphalt  cutbacL  Generally,  structures  that  were  originally 
designed  to  withstand  loadings  associated  with  weapon 
effects  have  ample  capacity  to  support  additional  dead  loads 
associated  with  berming.  This  is  not  the  case  for 
ccnverdonal  structures  and  analysis  is  required  to  determine 
the  need  for  additional  support 

3.  SANDBAGGING 

General  Sandbagging  is  a  traditional  method  of  providing 
effective  protection  to  walls,  overhead  structures,  and 
revetments  (DOA,  1985;  Hoot,  et  al.,  1974).  Sandbags  can 
also  be  used  to  construct  free  standing  wails  and  wall 
structures  for  protection  of  otherwise  exposed  assets.  Figure 
2  shows  a  sand  bag  upgrade  of  an  existing  structure. 

Threat  Protection.  Several  test  series  conducted  by  the 
U.S.  Army  Waterways  Experiment  Station  (WF.S)  have 
demonstrated  that  sandbags  placed  against  walls  and  over 
roofs  ate  successful  in  protecting  against  near  miss  and  direct 
hit  high  explosive  rounds,  and  direct  ballistic  impacts  (HcK't, 
et  al.,  1974;  Bucci  and  Mlakar,  1976;  Hamlin,  1986).  Placed 
against  walls,  sandbags  provide  protection  similar  to  that 
obtained  from  berming  with  soil. 

In  one  scries  of  tests  (Hoot,  et  al.,  1974),  a  timber  framed 
roof  structure  protected  by  4  layers  of  sandbags  (16  inches  of 
cover)  tvas  not  breached  by  contact  detonations  of  an  82  mm 
mortar,  107  mm  rxKket,  or  a  122  mm  r(x;ket.  Similar  tests 
of  6  inch  thick  precast  reinforced  concrete  roof  panels  also 
demonstrated  the  effectiveness  of  sandbagging.  An  82  mm 
mortar  round  detonated  directly  against  iho  slab  caused  major 
spalling  on  the  interior  surface.  A  107  mm  rocket  round 
caused  massive  spalling  and  breached  the  slab.  With  two 
layers  of  sandbags  covering  the  roof,  the  test  using  the  107 
mm  rocket  was  repeated  with  only  minor  cracking  occurring 
on  the  interior  face.  A  15  inch  layer  of  sandbags  stacked 
against  a  6  inch  precast  wall  provided  good  protection 
against  a  155  mm  HE  artillery  round  detonated  at  a  standoff 
of  5.0  ft.  A  similar  unprotected  wall  was  decimated  by  the 
round.  Sandbag  berms,  wail  and  roof  coverings  provide 
multiple  attack  protection  to  near  miss  high  explosive 
tlireats  similar  to  soil  berms  or  soil  covers. 

Limitations/Special  Considerations.  The  main 
limitation  of  sand  bags  is  aesthetics  and  for  this  reason, 
sandbagging  is  generally  not  considered  an  acceptable 
approach  for  more  permanent  upgrades.  Sandbags  have  had  a 
history  of  susceptibility  to  rot,  however,  newer  bag 
materials  made  of  an  acrylic  fabric  are  available  and  are 
reported  to  remain  serviceable  for  over  2  years  with  no  signs 
of  deterioration  under  all  climatic  conditions  (DOA,  1985). 
Improved  performance  against  ballistic  penetration  is 
obtained  by  mixing  the  fill  material  with  dry  portland 
cement  (1  pan  cement  to  10  parts  soil)  or  dipping  the  filled 
bags  in  a  cement-water  sloiiy  (DOA,  1985). 


4.  SAND  GRIDS 

General.  Sandgrids  were  originally  developed  as  a  soil 
confining  system  for  use  in  roadway  constraction  over  loose 
soils.  Sandgrids  are  constructed  by  filling  a  prefabricated 
plastic  form  shaped  like  cells  of  a  honeycomb  with  a 
granular  material  such  as  sand  or  gravel  (Figures  3, 4).  The 
sand  grids  are  available  in  a  standard  configuration  or  a 
newer,  notched  configuration.  The  notched  configuration 
allows  for  development  of  a  lapped  joint  between  layers  that 
prevents  leakage  of  the  fill  material.  The  principal 
advantages  over  the  more  traditional  methods  of  hardening 
with  soil  are  ease  of  construction  and  reduced  space 
requirements.  Currently  available  sandgrids  are  38  inches 
wide  in  place. 

Threat  Protection.  Testing  ha.s  demonstrated  that 
sandgrids  can  provide  efficient  and  effective  protection 
against  near  miss  general  purpose  bombs,  high  explosive 
artillery  rounds,  shoulder  launched  rockets,  and  machine  gun 
fire  (White,  1983;  Wood,  1985;  Hamlin,  1986;  Hayes, 
1987, 1988;  Hyde,  1989). 

Recent  tests  conducted  on  the  NATO  test  facility 
demonstrated  that  sandgrids  are  very  effective  in  protecting 
semi-hardened  walls  ffom  near  miss  general  purpose  bomb 
fragments  and  airblast  (Hayes,  1988;  Hyde,  1989).  In  this 
test  a  65  cm  wall  protected  by  a  sandgrid  shield  (Figure  4) 
suffered  no  damage  from  a  GP  bomb  detonation.  The  same 
waU,  unprotected,  suffered  severe  spall  when  subjected  to  the 
same  weapon  at  the  same  standoff. 

Wood  (1985)  reports  the  results  of  tests  on  free  standing 
sandgrid  revetments  against  105  mm  flechette  artillery 
roun^,  155  mm  HE  airbursts,  and  155  mm  surface  bursts. 
Six  105  mm  flechette  rounds  were  detonated  120  ft  from  the 
target.  None  of  the  fl;\hette$  completely  penetrated  the 
revetment.  Six  sta.idard  D544  fragmentary  EC  rounds  were 
detonated  at  a  stand  off  of  40  ft  and  an  elevation  of  20  ft 
The  revetment  remained  stable  ind  erect  and  showed  slight 
deterioration.  There  was  no  evidence  of  fragments  passing 
through  the  sandgrid.  Six  155  mm  rounds  each  were 
statifMy  detonated,  nose  tangent,  at  ranges  of  5, 10,  and  15 
ft  from  the  test  revetments.  At  the  5  ft  range,  the  six  rounds 
totally  collapsed  the  revetment.  Good  fragmentation 
protection,  however,  was  provided  for  the  first  two  rounds. 
At  the  10  ft  range,  after  six  rounds,  the  sandgrid  revetment 
remained  inv4Ct  and  erect  but  the  top  three  layers  were 
significantly  damaged.  No  fragments  penetrated  through  the 
revetment.  The  six  rounds  fired  at  the  IS  ft  range  caused 
only  superficial  damage.  Sandgrid  revetments  have  also 
demonstrated  effectivene.ss  in  defeating  penetration  by  small 
arms  fire  up  to  50  caliber  rounds  (White,  1983). 

Gravel  filled  sandgrid  revetments  have  been  tested  against 
U.S.  LA^V  and  Soviet  RPG-7  antitank  weapons  (Hayes, 
1988).  Crt  hed  limestone  (3-4  in,  diameter)  was  used  for  the 
revetments  tested  against  the  LAW,  and  rounded  river  gravel 
(1.5-2  in.  diameter)  was  used  for  the  revetment  tested  against 
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the  RPG-7.  Two  consecutive  hits  by  LAW  rounds,  in 
approximately  the  same  Ircadon,  cratered  the  front  face  of 
the  revetment  but  the  revetment  was  not  perforated.  The 
revetment  hit  by  the  RPG-7  w?s  perforated  by  the  rocket 
motor,  but  the  residual  velocity  of  the  motor  was  not 
sufficient  to  cause  damage  to  a  projected  wall  or  lightly 
armored  assets  placed  behind  the  revetment. 

Sandgrids  filled  with  sand  and  gravel  have  also  been  tested 
against  a  wire  guided  warhead  (TOW)  placed  against  the  side 
of  the  revetment  and  statically  detonat^  (Hayes,  1988).  The 
sand  filled  revetment  was  breaclied  by  die  weapon  and  the 
steel  witness  plate  behind  the  revetment  was  perforated.  The 
gravel  filled  revetment  was  breached  and  fragments  from  the 
weapon  perforated  the  revetment;  however,  the  witness  plate 
was  not  damaged. 

Sangrid  revetments  provide  limited  second  attack  protection 
from  near  miss  GP  bombs  based  upon  evidence  from  the 
NATO  test  facility.  The  sandgrid  form  was  severely 
dama;»-:d,  however,  much  of  the  soil  remained  in  the  shape 
of  a  small  berm.  In  other  applications  against  less  severe 
threats  such  as  near  miss  HE  artillery  rounds,  the  .sandgrid 
revetment  provides  excellent  .second  attack  protection.  As 
overhead  protection,  sandgrids  provide  limited  second  attack 
protection  against  direct  hit  mortar  and  HE  artillery  attack 
but  can  be  repaired  with  loose  soil  oi  sand  bags. 

Limitations/Special  Considerations.  The  maximum 
fiee  standing  height  for  use  as  a  protective  structure  against 
conventional  weapons  is  about  8  ft. 

5.  CONCRETE  MODULAR  REVETMENTS 

General.  The  discussion  in  this  section  focuses  on  several 
modular  designs  constructed  of  conventionally  reinforced 
concrete  (R/C)  and  SIFCON  (Slurry  Infiltrated  Fiber 
reinforced  Concrete)  that  have  been  tested.  Figure  S  shows 
one  common  type  of  modular  unit  The  revetments  ore 
often  bermed  by  soil  to  improve  overall  performon  .e  (see 
Section  2  and  Figure  1). 

Threat  Protection.  These  revetments  provide  protection 
against  fragments  and  airolast  from  near  miss  general 
purpose  bombs  arid  other  lesser  threats  such  as  HE  artillery 
rounds,  rockets,  and  mortars.  Also,  they  ore  often  used  to 
deny  iine-of-sight  to  doors  and  other  vulnerable  openings. 

Recent  tests  on  R/C  revetntents  include  (he  CHEBS  series 
conducted  at  Kiiiland  AFB  (Hyndman  and  Bulunan,  1987; 
Carson  and  Morrison,  1986)  and  the  NATO  tests  conducted 
at  Tyndall  AFB  (Hyde,  1989).  The  CHEBS  tests  included 
revetments  that  were  both  bemted  and  unbermed  (see  Section 
2).  In  the  CHEBS  9  and  10  tests  four  standard  revetment 
designs  were  tested  in  a  variety  of  configurations:  (1)  the 
Bitburg  Revetment  with  a  wall  thickness  of  0.30  m,  (2)  the 
3-me(er  Aircraft  Roveur.ent  with  a  wall  thickness  of  0.24S 
m,  (3)  the  a-meter  .\ircraft  Revetment  with  a  wall  thickness 
of  0.24S  m,  and  (4)  the  4-roeter  Au'crafl  Revetment  with  a 


wall  tapering  fre.m  0.245  m  at  the  base  to  0.085  m  at  the 
top.  A  second  scries  of  tests  (CHEBS  16)  was  conducted 
using  similar  designs  constructed  using  SIFCON  (Carson 
and  Morrison,  1986).  The  NATO  test  examined  the 
effectiveness  of  unbermed  Bitburg  revetments  subjected  to 
tl<e  NATO  GP  bomb  threat.  R/C  revetments  have  dso  been 
tested  against  a  variety  of  other  conventional  weapons 
mciuding  HE  rockets,  mortars,  and  machine  gun  fire  (Hoot, 
etal.,  1974). 

The  results  of  the  CHEBS  and  NATO  tests  demonstrated 
that  unbermed  R/C  revetments  can  provide  good  first  strike 
protection  but  limited  second  attack  protection  due  to  tlieir 
susceptibility  to  fragment  damage  and  movement  (i.e.,  rigid 
body  shifting).  As  bermed  structures,  R/C  revetments 
provide  excellent  second  attack  protection. 

For  the  revetments  constructed  of  SIFCON,  the  general 
damage  characteristics  resulting  from  fragment  impact 
differ^  significantly  from  the  damage  characteristics  of  the 
P/C  revetments,  llte  damage  on  the  front  face  or  the 
SIFCON  revetments  tended  to  be  localized  immediately 
around  the  impact  area  and  not  cratered  as  was  the  case  with 
the  R/C  revetments.  Little  or  no  rear  spall  occurred  with  the 
SIFCON  revetments.  When  the  revetments  were  perforated 
by  fragments,  the  penetrations  were  clean,  nearly  cylinuiical 
and  easily  repaired. 

R/C  revetment  panels  6  inches  thick  and  supported  by 
precast  blocks  (Hoot,  et  al.,  1974)  were  demonstrated  to  be 
effective  against  81  mm  mortars  and  120  mm  rockets 
detonated  at  a  standoff  of  5.0  ft  Against  155  mm  HE 
artillery  rounds  the  6  inch  panel  was  considered  to  be 
effective  only  for  standoffs  greater  than  30  fL  Against  122 
mm  rockets,  the  revetment  was  considered  .satisfactory  for 
standoffs  of  iO.O  ft;  and  with  the  addition  of  a  IS  inch 
sandbag  lacing,  the  revetment  defeated  the  frapient  effect  at 
a  5  ft  standoff. 

Limitattorts/Speciai  Considerations.  R/C 
revetments  ore  limited  by  the  resources  necessaiy  to  fabricate 
and  deploy  them.  Their  use  as  on  expedient  measure  requires 
that  they  be  prefabricated  and  available  for  deploymenu 
Modular,  prefabricated  revetments  reqtiire  a  reiotivdy  smooth 
base  surface  for  deployment  Bermed  R/C  revetments  have 
faLiy  large  space  requirements  and  other  limitations 
associated  with  soil  berms  or  bermed  walls  (see  Section  2). 
The  geometry  of  the  Aircraft  Revetment  dou  not  allow  for 
the  formation  of  ninety  degree  turns  In  the  layout  of  an  array 
of  revetments.  The  base  of  the  Bithuif  revetment  has  comers 
mitered  k  45  degrees  which  allow  ioi  ninety  degree  turns, 
but  these  comers  ore  vulnerable  to  fragment  penetration 
because  tlte  wall  thicknesses  of  the  two  revetments  just  meet 
and  do  not  lap.  These  comers  should  be  protected  with 
sandbags. 
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6.  BIN  REVETMENTS 

General.  Bin  revetments  refer  to  any  of  a  variety  of 
methods  used  to  create  vertical  walls  of  sand,  soil,  gravel,  or 
rock  rubble.  These  systems  combine  the  protective  qualities 
of  soil  structures  with  an  efficient  use  of  space.  The 
thickness  of  the  soil  wall  is  the  primary  means  of  providing 
protectfon.  The  structural  system  is  designed  to  confine  the 
soil  and  can  be  constructed  of  reinforced  concrete,  steel, 
wood  or  wire  caging.  Typical  configurations  are  two  parallel 
walls  with  fill  between  them  or  prefabricated  containers 
filled  with  soil  and  arranged  into  a  revetment  (Figure  6).  In 
addition,  bin  revetments  can  take  the  form  of  planters  for 
aesthetic  permanent  upgrades. 

Threat  Protection.  These  revetments  provide  essentially 
the  same  {sotection  as  bermed  walls  or  revetments.  They  can 
be  employed  against  fragment  and  airblast  threats  fhim  near 
miss  general  purpose  bombs,  other  lesser  conventional 
munitions,  and  ballistic  threats.  They  can  also  be  used  to 
deny  line-ofsight  to  doors  and  other  vulnerable  openings. 

Soil  bin  reveunents  with  a  soil  thickness  of  12  inches  were 
tested  by  WES  for  protection  against  near  miss  mortar  and 
rocket  detonation  in  support  of  the  U.S.  Army  during  the 
Vietnam  conflict  (Carre,  1969,  1972).  Soil  bin  structures 
constructed  of  plywood,  18  gage  corrugated  metal,  and 
M8A1  landing  mats  were  tested.  The  M8A1  soil  bin  defeated 
aU  fragmenu  from  81  to  120  mm  mortars  and  107  mm 
rocket  detonations  at  a  stand  off  of  S  ft.  The  plywood  and 
corrugated  metal  bin  revetments  were  also  dfective  at  a 
range  of  S  ft  but  suffered  greater  damage. 

As  a  part  of  the  CHEBS  9  and  10  tests  (Hyndman  and 
Bultman,  1987),  a  soil  bln  revetment  was  constructed  using 
two  parallel  rows  of  3  m  Aircraft  revetments  (see  Section  S). 
The  soil  thidmess  was  approximately  S  ft.  The  soil  bin  was 
positioned  SO  ft  from  a  1000  lb  CP  bomb.  The  face  of  the 
front  Aircraft  revetment  was  severely  damaged  by  the 
fragment  impacts  but  no  damage  to  the  rear  revetments 
forming  the  soil  bin  was  reported. 

C  AC  Chamber  revetments  were  tested  as  a  port  of  the 
NATO  facility  wsts  (Hyde,  1989).  This  soil  bin  structure 
makes  use  of  off>ihe>shelf  precast  concrete  manhole  liners 
(rectangular  tubes  3.0  x  4.5  ft  in  section,  2.46  ft  long,  with 
a  wall  thickness  of  3  inches).  The  tubes  allow  for  a  soil 
thickness  of  2.S  ft.  Several  reinforcing  methods  were  tested 
including  wire  mesh,  synthetic  and  steel  fibers  and  various 
combinations.  Revetment  walls  were  constructed  by  bolting 
the  tubes  together  with  steel  straps  and  filling  the  units  with 
sand.  Walls  one,  two,  and  duee  layers  high  were  tested.  The 
tests  demonstrated  that  this  configuration  can  provide 
fragment  protection  against  general  purpose  bombs.  If 
reinforced  with  steel  mesh,  multlplo  nttiick  protection  con  be 
expected. 

Liotitatioa/Special  Consideration.  Soil  bin 
revetments  generally  require  significant  construction 


resources.  If  built  of  expedient  materials,  they  tend  to  be 
temporary  measures.  Constructed  of  reinforced  concrete  or 
masonry,  they  are  usually  part  of  permanent  upgrades  or  new 
facility  construction. 


7.  SACRIFICIAL  PANELS 

General.  Sacrificial  panels  are  panels  attached  to  the 
exterior  of  a  main  structural  wall  so  that  an  air  space  is  left 
between  the  main  wall  and  the  panel.  The  panels  are  not 
expected  to  survive  a  weapons  blast  but  still  provide 
si^ificant  protection  to  the  main  wall  (Figure  7),  Sacrificial 
panels  constmeted  of  a  variety  of  materials  such  as  reinforced 
concrete,  steel  plate,  plywood,  or  a  layered  combination  of 
materials  have  been  tested.  Sacrificial  panels  are  used  for 
both  expedient  upgrading  of  existing  facilities  and  as 
inexpensive  hardening  methods  for  new  construction. 

Threat  Protection.  Sacrificial  panels  have  demonstrated 
effectiveness  in  proteedng  against  fragments  and  airblast 
from  near  miss  general  purpose  bombs  and  HE  rounds. 
Recent  tests  include  a  series  of  19  scaled  tests  on  a  variety  of 
panel  types  (reinforced  concrete,  steel,  plywood,  and 
composite)  conducted  by  the  AFESC  and  the  WES  (MeVay, 
1988),  the  full  scale  NATO  tests  conducted  at  Tyndall  AFB 
(MeVay,  1988;  Hyde,  1989)  and  a  series  of  tests  on  concrete 
panels  conducted  by  WES  (Colthoq),  1937). 

In  all  of  these  cases  the  damage  to  the  walls  protected  by  the 
sacrificial  panels  was  greatly  reduced  from  the  damage 
observed  for  the  unprotected  walls.  For  example,  for  a  w^ 
protected  by  a  6  inch  precast  panel  with  a  nominal  1  inchair 
space  at  the  full  scale  NATO  test,  only  minimal  scabbing 
and  cratering  occurred  on  the  front  face  of  the  main  wall  and 
no  spall  occuned  on  the  backfoce.  The  same  wall  tested  bare 
suffered  severe  spall  damage  that  would  have  been  lethal  to 
equipment  and  personnel  inside.  Ihe  mechanism  behind  the 
protection  afforded  is  due  to  the  greatly  reduced  number  and 
momentum  of  the  fragments  that  impact  the  main  wall  and 
greatly  reduced  peak  pressures  on  the  main  wall.  This 
results  in  a  much  lower  potential  for  spall  and  breach  even 
though  load  duration  from  airblast  may  be  greatly  increased. 

Limitatloos/Speciai  Considerations.  The  main 
limitation  of  sacrificial  panels  is  their  single-bit*oniy 
capability.  If  multiple  hit  capability  is  required  a  system 
must  be  in  place  to  provide  for  expedient  teplocoment  of  the 
panels. 

8.  CONCLUSIONS 

A  wide  range  of  methods  may  be  used  to  provide  expedient 
structural  hardening  against  the  effects  of  conventional 
weapons.  These  method  vary  in  the  degree  of  protection 
afforded,  multiple  strike  capability,  cose  of  construction,  and 
ease  of  deployment.  Optimal  selection  of  a  hardening 
method  reqtiiies  careful  consideration  of  each  of  these  factors 
os  they  relate  to  the  particular  needs  of  a  given  situation  and 
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the  resources  available.  From  the  test  data  reviewed  the 
most  effective  protection  is  provided  by  earth  structures. 
Earth  structures  can  significantly  reduce  the  effects  of 
fragments  and  airblast  and  can  provide  excellent  multiple 
strike  protection.  Confining  the  earth  through  the  use  of 
bags,  sand  grids,  or  bins  eliminates  many  of  the  limitations 
associated  with  traditional  berms.  Sacrificial  panels  and 
unbermed  concrete  revetments  also  provide  effective 
protection  but  do  not  provide  the  multiple  hit  capability  of 
earth  structures. 
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Figure  1.  Benued  Wall 
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Figure  2.  Sandbag  Upgrade 


Figure?.  Sacrificial  Panel 
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ABSTRACT 

High  explosive  threats  directed 
against  U.S.  facilities  throughout  the 
world  have  Identified  a  need  to  blast 
harden  these  facilities.  Lose  of 
windows  during  blast  loading  poses  a 
serious  threat  to  personnel  and  facility 
security.  Research  has  shown  that  a 
blast  resistant  window  Is  feasible  and 
that  a  slngle>degree>of*freedon  (SOOF) 
analysis  can  predict  the  window 
response . 

The  WES  generic  blast  resistant 
window  was  designed  to  withstand  a  peak 
reflective  pressure  of  80  pel.  This 
window  survived  the  test  with  no 
damage .  The  SDOF  predictions  for  the 
26>lnoh  by  26 'Inch  windows  were  In  good 
agreement  with  test  data,  with  the 
analysis  being  oonaervatlve.  This  good 
agreement  resulted  from  the  use  of  a 
resistance  function  developed  from  a 
static  test. 


INTRODUCTION 

High  explosive  threats  directed 
against  U.S.  facilities  throughout  the 
world  have  identified  a  need  to  blast 
harden  these  facilities.  One  of  the 
most  vulnerable  elements  in  these 
structures  is  the  window  systems.  Lose 
of  windows  during  blast  loading  poses  a 
serious  threat  to  personnel  and  facility 
security.  Research  has  been  funded  for 
the  purpose  of  designing  a  generic 
window  system  to  withstand  high  blast 
pressures.  From  this  research,  it  was 
shown  that  a  window  systsm  can  be 
dosigned  to  resist  high  blase  pressures 
and  that  the  dynamic  response  can  be 
predicted  with  a  Bingle>dagraa-of< 
freedom  (SDOF)  analysis. 

UBS  BLAST  UINDOU 

The  UBS  generic  blast  resistant 
window  was  designed  to  withstand  a  peak 
reflaotive  pressure  of  80  pel.  The 
window  was  designed  in  three  parts;  the 
anchoring  frame,  which  is  oast  into  the 


reinforced  concrete  test  wall,  a  rigid 
window  frame,  and  1-1 /4-inch-thick 
polycarbonate  glazing.  Figure  1  shows  a 
schematic  view  of  the  design  concept. 


The  critical  element  in  designing  the 
anchoring  frame  was  to  develop  a  large 
shear  area  across  the  thickness  of  the 
concrete  test  wall.  The  frame  was 
designed  to  be  positioned  in  the  wall 
and  cast  Ln  place,  thereby,  developing  a 
shear  area  l.S  times  larger  chan  the 
thickness  of  the  wall  (12-inohea  thick). 
The  anchoring  frame  was  fabricated  using 
A •inch  by  8 'inch  by  1/2 'Inch  thick  steel 
angle  with  a  minimum  yield  stress  of 
36,000  psi.  To  anchor  the  frame  into 
the  concrete,  5/8 -inch  diameter  by  8- 
Inoh'long  concrete  anchors  were  welded 
to  the  frame.  These  anchors  were 
attached  to  the  frame  at  a  cantarllna 
spacing  of  4 •inches.  The  frame  was 
positioned  and  #3  steel  reinforcing  bar 
stirrups  were  added  in  the  high  shear 
region  of  the  slab  adjacent  to  the 
frame. 

The  test  wall  consisted  of  two 
layers  of  S\8-inoh'-diaaeear  reinforcing 
bars  placed  5- inches  on  center  in  both 
directions  and  in  both  faces  of  the 
slab.  The  shear  stirrups  were  added  for 
a  distance  of. 1  foot  from  the  anchor 
frame.  The  concrete  design  used, in  the 
test  wall  had  a  minimum  yield  (f  q)  of 
SOOO  psi  at  28  days. 
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Tho  window  frame  waa  designed  to 
be  structurally  rigid;  therefore,  the 
remaining  critical  elements  In  the 
design  of  the  window  were  the  bite 
placed  on  the  glazing  and  the  type 
glazing  used.  The  frame  was  fabricated 
using  4-  by  3-  by  1/2 -Inch  steel  angle, 
with  the  4- Inch  edge  providing  the  bite 
of  the  glazing.  Figure  2  shows  the 
design  of  the  window  frame.  Vlth  this 
design,  a  bite  of  1-1/4 -Inches  was 
placed  on  the  glazing. 


Fig.  2 

The  glazing  selected  for  the 
window  was  1-1/4- Inch  lamlnaead 
polycarbonate  (SP-12S0).  This  material 
was  selected  for  its  ability  to  undergo 
large  deflactlona  and  remain  undamaged. 
The  SF'1230  contains  four  layers  of 
polycarbonate  consisting  of  two  -  1/2- 
inch  layers  placed  in  the  center  and  two 
•  1/8 -Inch  layers  placed  on  the  outside. 


Fig.  3a 


ky 


ret) 


Fig,  3b 


AtlALYSlS 

The  determination  of  the  dynamic 
response  of  a  sinpla  structural  system 
using  numerical  procedures  is  presented 
in  detail  In  References  1,  2,  and  3. 

More  complex  systems,  such  as  the  window 
system  reported  herein,  can  also  be 
analyzed  with  an  SDOF,  provided  an 
accurate  represontatlon  of  the  load 
function  P(t),  reslstanca  function 
(load-daf lection  curve),  and  mass  <M) 
can  bo  obtained, 

SDOF  Analysis 

The  selection  of  the  idealized 
spring-mass  system  in  Figure  3a  is  such 
that  the  deflection  of  the  mast,  y,  is 
the  aasM  as  the  centerline  deflection  of 
the  window  glazing.  From  the  freebody 
diagram  shown  in  Figure  lb,  the  equation 
of  motion  is  derived  and  then  solved 
numerically. 


To  solve  for  the  dynamic  response 
of  the  window,  a  static  test  was 
conducted  to  determine  the  resistance 
function,  the  P(c)  was  calculated  from 
formulas  developed  in  ftaferenee  2,  and 
the  mass  of  the  polycarbonate  glazing 
was  determined.  Once  these  three  values 
were  obtained,  the  use  of  a  computer 
code  <SDOF)  developed  at  UBS  waa 
utilized  in  preforming  the  analysis. 

Static  Test  Device 

To  conduct  the  static  testing,  UBS 
designed  and  constructed  a  static  test 
device  as  shown  in  Figure  4,  The  test 
device  was  designed  to  subject  a  test 
article  to  a  maximum  hydrostatic 
pressure  of  200  psl.  The  parts  of  the 
device  consist  of  a  U-framad  base,  a 
test  wall  slab,  and  the  hold-down  slab. 
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This  Cttt  d«vlc«  «llows  for  the 
testing  of  a  window  as  If  the  test 
article  were  Installed  In  a  building 
wall.  Once  the  window  Is  Installed,  the 
test  wall  Is  positioned  on  the  U-frane 
base  and  the  hold-down  slab  positioned 
on  top.  The  two  slabs  are  then  bolted 
together.  In  this  configuration,  a 
chaaber  Is  fomed  becwaan  the  window  and 
the  bottom  of  the  hold-down  slab.  It  Is 
Into  this  chamber  chat  water  Is  pumped 
to  produce  the  uniform  hydcostatlc 
pcataure.  A  maximum  of  300  psl  has  been 
achieved  with  this  setup. 

By  utilising  the  U-shapad  base, 
access  to  the  underside  of  the  window  Is 
possible  and  placement  of 
Instrumentation  can  be  accomplished. 

The  design  also  allows  clearance  for 
vldaotaps  equipment  to  be  used  In 
recording  the  response  of  the  test 
article . 

Static  Testing 

The  window  system  was  subjected  to 
a  hydrostatic  pressure  chat  produced 
failure  of  the  window.  Failure  Is  taken 
to  be  Che  point  where  the  glaslng  or 
window  system  cannot  sustain  additional 
load.  In  the  test  oonduoted,  the 
failure  mechanism  for  the  window  was  the 
glaslng  slipping  from  the  frame  bits. 
Ulth  the  data  collected  on  the  pressure 
and  oanterpolnt  dafleotlon,  a  load- 
deflection  curve  was  derived  for  the 
analysis,  figure  3  shows  the  resistance 
function  developed  from  Che  test  data. 


i«mHi  tut 

Fig.  5 


SDOF  Code 

The  procedure  used  In  the  SDOF 
coda  Is  described  In  Reference  1.  The 
procedure  Is  referred  to  as  the 
constant-velocity  or  lumped- Impulse 
procedure.  The  code  requires  the 
following  Information  to  perform  the 
calculation;  mass  (H) ,  area  <A) ,  load- 
mass  factor  (K^t)  I  percent  of  doping 
(c),  elms  step  (d^) ,  resistance  function 
f(R),  and  forcing  function  F(t>. 

In  ths  analysis  of  the  window 
system,  the  mass  was  determined  by 
weighting  Che  glazing.  The  area  was 
selactsd  as  Che  clear  span  area  of  the 
window  <26  Inches  by  26  Inches) .  The 
values  for  the  load-mass  factor  and 
percent  of  damping  were  chosen  to  be 
0.67  and  0.03,  raspeotlvaly.  The 
numerical  Iterative  time  step  of  0.0001 
seconds  was  used  for  the  SDOF 
calculation. 

Ths  analysis  was  performed  by 
using  a  reflsccsd  pressure-time  history, 
with  clearing  time,  for  a  hemispherical 
high-explosive  <HE)  charge  producing  a 
peak  pressure  of  52. 6  psl.  The  negative 
phase  peak  pressure  was  -4.0  psl. 

Figure  6  shows  Che  calculated  pressure - 
time  history. 


Utilizing  the  Input  data,  the  SDOF 
analysis  was  oonduoted. 
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DYHAMIC  TEST 

To  v«rl£y  the  analysis,  a  dynaaic 
taaC  was  conducted  on  five  WES  blast 
windows .  The  window  clear  span  openings 
were  two  -  26-lnch  by  26-lnch,  two  -  36* 

Inch  by  36-lnch,  and  a  40*inch  by  40- 
Inch. 

The  test  was  conducted  and 
subjected  the  two  36-lnch  by  36-lnch  to 
an  average  peak  reflective  pressure  of 
49.2  psl  and  the  two  26-inch  by  26-lnch 
and  the  40 -Inch  by  40- inch  windows  to  an 
average  peak  reflective  pressure  of  S7.4 
psl. 

COHFARISON  OF  CALCUL4TED  RESPONSE  AND 
TEST  DATA 

In  deceralnlng  the  dynanlc 
response  of  the  three  sizes  of  blast 
windows,  an  analysis  was  perfomed  using 
Che  static  load-deflection  curve 
generated  for  the  26-lnch  by  26-lnch 
window  and  a  theoretical  resistance 
function  calculated  for  the  36 -Inch  by 
36-lnch  and  40-lnch  by  40-lnch  windows. 

The  theoretical  resistance  functions 
were  generated  froa  equations  developed 
In  Reference  1. 

The  analysis  pcedleted  a 
cancerpolnc  deflection  of  1.81  Inches 
for  the  two  26<lnch  by  26-lnch  windows, 

6.1  inches  for  the  two  36-lnch  by  36- 
Inch  windows,  and  graatsc  than  10  Inchaa 
fee  the  40-lnch  by  40-lnch  window,  A 
ooapacliion  of  the  pcecsst  predlotlon  and 
aotual  test  data  la  shown  In  Table  1. 

Table  1  Suwscy  of  Oofleotlon  Data  for  USS 
Blast  Windows 


Window 

Sizes  and 

Type  PC 

Predicted 
Deflection 
with  SDOF 
(In) 

Test  Data 

Summary 

(In) 

26x26  New  PC 

1,84 

1.30 

1.35 

26x26  Old  PC 

1.84 

1.32 

1.34 

36x36  New  PC 

6.1 

3.25 

3.10 

36x36  Old  PC 

6.2 

2.92 

3.15 

40x40  Hew  PC 

10 

4.00 

3.98 

From  comparison  of  the  results  from  the 
test.  It  appears  that  revision  of  the 
calculated  resistance  functions  will  be 
required.  The  test  data  shows  that  the 
polycarbonate  glazing  has  a  higher 
stiffness  (resistance)  than  predicted  by 
the  theoretical  resistance  function. 

A  posttest  analysis  was  performed 
using  the  data  from  the  26-lnch  by  26- 
lnch  window.  The  pressure- time  history 
used  In  the  analysis  was  the  average  of 
the  records  from  the  six  pressure  gages 
used  In  the  test.  From  this  analysis 
the  peak  centerpolnt  deflection  was 
calculated  to  be  1.58  Inches.  Although 
this  comparison  Is  based  on  one  test.  It 
Is  concluded  that  the  proposed  procedure 
can  be  used  to  verify  or  determine  the 
blast  resistance  of  a  window  system. 

The  good  agreement  between  the  predicted 
value  and  actual  test  data  suggests  that 
the  resistance  function  developed  from 
static  testing  correctly  described  the 
spring  (k)  of  the  SDOF  system. 

The  poor  agreement  of  the 
predicted  spring  (k)  and  test  response 
data  Indicates  that  the  boundary 
conditions  are  leore  complicated  than 
assumed  for  the  theoretical  resistance 
function.  This  Is  one  area  where 
further  studies  are  nsedad. 
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ABSTRACT 

Finite  element  modeling  Is  an  analysis 
technique  which  is  commonly  used  In  the  design  and 
analysis  of  structural  systems.  The  approach  has 
been  applied  In  conventional  weapon  effects 
research  particularly  In  the  detailed  analysis  of 
structural  subsystems  and  components.  Limited 
studies  have  been  conducted  using  finite  element 
models  to  simulate  the  complete  structural  system 
subjected  to  conventional  weapon  effects.  This 
paper  addresses  the  analysis  of  three-story  pre¬ 
cast  and  reinforced  concrete  structures  subjected 
to  Interior  detonations  of  general  purpose  bombs. 


TEST  SERIES  DESCRIPTION 

In  support  of  OSD  tri -service  requirements, 
AFATL  and  3246  Test  Wing  conducted  a  test  series 
which  addressed  the  effects  of  general  purpose  air 
delivered  munitions  against  standard  office 
building  construction.  Specifically,  the  tests 
Involved  three-story  structures!  two  of  which  were 
constructed  using  pre-cast  concrete  panels  and  the 
third  representing  cast-ln-place  reinforced 
concrete  construction.  Figure  1  shows  an  Isometric 
view  of  the  test  structures  while  Figure  2 
Illustrates  the  floor  layout  and  building 
dimensions.  All  wall  elements  consisted  of  6-inch 
thick  reinforced  concrete.  The  pre-cast  panels 
were  assembled  and  fastened  using  standard  pre-cast 
construction  techniques  including  the  use  of  weld 
plates  as  the  connection  type  between  panels. 
Figure  2  also  shows  that  each  building  was 
constructed  In  two  segments  with  a  dividing  wall  in 
the  center.  The  only  structural  connection  between 
the  two  building  halves  was  the  foundation  system. 
Separating  the  buildings  Into  two  separate  sections 
allowed  testing  to  be  conducted  on  each  building 
section. 

The  tests  were  conducted  using  dynamic  and 
static  placement  of  the  test  weapon.  For  dynamic 
placement  tests,  a  rocket  sled  track  propelled  the 
weapon  Into  the  first  floor  of  the  structure.  The 
weapon  fuzing  was  designed  such  that  the  weapon 
penetrated  the  building  and  detonated  in  the  first 
room  as  shown  by  the  weapon  trajectory  line  In 
Figure  2.  In  the  static  tests,  the  event  consisted 
of  placing  the  munition  in  the  desired  room  and 
then  remotely  detonating  the  weapon. 


Figure  2:  Test  Structure  Floor  Plan 
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In  several  of  the  events,  the  test  structures 
were  instrumented  with  accelerometers  and  pressure 
transducers.  The  transducer  array  was  designed  to 
monitor  the  weapon  loading  on  the  structure, 
structural  response,  and  the  documentation  of  the 
Interior  environment.  In  addition  to  the 
instrumentation  array,  exterior  and  Interior 
cameras  documented  the  building  collapse/damage  and 
the  Interior  building  environment.  In  selected 
events,  the  test  structure  interior  included  office 
furniture,  computer  equipment,  and  test  manniquins. 


ANALYSIS  OVERVIEW 

A  parallel  effort  to  the  test  series  was  the 
finite  element  modeling  of  the  test  structures. 
The  Intent  of  the  study  was  to  identify  aspects  of 
the  finite  element  modeling  that  seemed  promising 
and,  more  Importantly,  limitations  and  other  areas 
requiring  further  development.  A  second  purpose  of 
the  analytical  study  was  to  support  the  test 
Instrumentation  setup  by  providing  pre-test 
predictions. 

The  analysis  effort  consisted  of  two  general 
phases:  pre-test  modeling  effort  and  the  post-test 
analysis.  The  pre-test  analysis  centered  on 
providing  the  pre-test  predicting  for  the 
Instrumentation  system  and  Included  the  development 
of  models  for  the  entire  structure,  substructure 
systems,  and  components.  The  pre-test  modeling 
effort  used  COSMOS?  a  general  purpose  finite 
element  code.  The  post-test  analysis  focused  on 
analytlng  selected  substructure  systems  and 
components  using  AOINA*  (Automatic  Dynamic 
Incremental  Nonlinear  Analysis)  finite  element 
code.  The  post-test  analysis  Included  nonlinear 
material  models  and  direct  comparisons  with  test 
results. 


PRE-TEST  analysis 

The  primary  task  in  the  pre-test  analysis  was 
the  development  of  finite  element  models  for  the 
entire  structure,  substructure  systems,  and 
components,  figure  3  shows  the  first  floor  segment 
of  the  entire  structure  finite  element 
representation.  The  models  were  developed  using 
COSMOS  Interactive  environment  and  were  constructed 
using  wall  and  floor  panels  as  the  basic  building 
blocks  for  the  mesh.  Connectivity  between  the  wall 
and  floor  panel  finite  element  blocks  replicated 
the  connection  between  the  wall  and  floor  panels  in 
the  pre-cast  test  structure,  for  the  casi-ln-place 
structure  model,  the  finite  element  mesh 
represented  a  monolithic  structural  system.  The 
substructure  models  represented  a  series  of  wall  or 
floor  panels  while  the  component  models  consisted 
of  single  wall  or  floor  panels.  In  each  case,  the 
natural  frequencies  and  mode  shapes  of  the  finite 
element  models  were  determined,  figure  a  shows  the 
first  three  mode  shapes  of  a  substructure  iysieiu 
consisting  of  three  wall  panels. 


Figure  3:  First  Floor  Segment  of  Test  Structure  Finite  Element 
Nooei 


Figure  ai  First  Three  Moat  Shepts  For  util  Panel  Subsystn 


During  the  development  of  the  finite  element 
models,  free-vibration  tests  were  conducted  on  the 
test  structures.  The  free-vibratlon  tests 
consisted  of  the  static  loading  of  the  test 
structure  at  a  selected  location  and  then  suddenly 
unloading  the  structure.  The  free-vibration 
response  of  the  structure  was  observed  by  an 
accelerometer  transducer  array.  The  pull -tests 
allowed  the  direct  observation  of  building 
substructure  natural  frequencies.  The  natural 
frequencies  observed  compared  well  with  the  natural 
frequencies  computed  by  the  finite  element  model. 

The  final  effort  of  the  pre-test  analysis  was 
to  predict  the  structural  response  for  the 
acceleration  transducer  array.  The  analysis 
consisted  of  estimating  blast  loading  on  the 
structure  caused  by  the  weapon  detonation  and  using 
the  modal  representation  of  the  structure  to 
predict  peak  structural  accelerations  at  the 
transducer  locations.  The  blast  loading  were 
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estimated  using  an  empirical  approach  which  using 
an  initial  impulse  loading  equal  to  the  first  three 
reflections  of  the  blast  and  a  time  varying 
pressure  load  equal  to  the  pseudo-static  confined 
blast  pressure. 


WEAPON  EFFECTS  TESTS 

The  weapon  effects  tests  conducted  on  the  test 
structures  provided  excellent  information  regarding 
the  general  vulnerabilities  of  the  two  construction 
types  tested,  collapse  mechanisms  for  pre-cast  and 
cast-in-place  structures.  blast  pressure 
propagation  through  unprotected  facilities, 
structural  response,  and  interior  environment 
definition.  Figure  5  shows  the  Interior  detonation 
of  the  test  munition  on  the  test  structure  while 
Figure  6  shows  a  typical  interior  pressure  time 
history  recorded  during  the  test  series.  The  test 
data  was  used  extensively  in  the  post-test  analysis 
which  is  discussed  in  the  next  section. 


POST-TEST  ANALYSIS 

The  post-test  analysis  centered  on  the 
analysis  of  selected  substructure  systems  and 
components  monitored  during  the  test  series  with 
active  instrumentation.  The  finite  element 
modeling  described  in  this  section  was  accomplished 
using  AOINA.  The  analysis  included  linear  and 
nonlinear  material  models  as  well  as  the 
elemlnation  of  structural  connections  and  elements 
during  the  loading  event.  The  analysis  used  the 
pressure  loading  measured  during  the  test  series 
(Figure  S)  and  standard  time  integration  methods  to 
predict  the  response  of  structural  subsystems, 
components,  and  connections.  Figure  6  shows  a 
typical  response  of  wall  panel  section  subjected  to 
an  airblast  loading. 

For  the  pre-cast  construction,  the  weak  points 
or  failure  modes  centered  on  the  panel  connections 
and  the  collapse  of  the  floor  panels.  For  the 
ca$t-1n-place  construction,  the  predicted  failure 
mode  followed  the  classic  shear  and  flexural 
failure.  Shear  failures  are  localiaed  near  the 
explosive  source  while  flexural  failure  is  more 
widespread.  In  the  vicinity  of  the  explosive 
source,  the  finite  element  representations  in  both 
cases  (pre-cast  and  cast-in-place)  quickly  exceed 
the  model  limitations  but  the  models  do  approximate 
when  the  structural  component  or  connection  fails 
(given  the  correct  failure  stresso';  or  forces). 
Future  analytical  efforts  should  focus  on  aoiha's 
capability  to  kill  elements  and  incrcmemaiiy 
progress  through  a  collapse  event.  Reasonable 
micro-  to  macro-  models  could  be  used  to  predict 
and  evaluate  conventional  weapon  effects  against 
structures  provided  that  key  empirical  models  are 
incorporated  into  the  analysis  process. 


Figure  5:  Interior  Detonation  on  Test  Structure 


Figure  6j  Typical  Interior  Pressure  Time  History 
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SUMMARY 

The  study  has  highlighted  the  capabilities  and 
limitations  of  the  finite  element  approach  In 
modeling  weapon  effects  against  structures.  The 
advantages  include  great  flexibility  in  model 
characterization,  generalized  loading  functions, 
nonlinear  analysis,  stability  analysis,  and 
elimination  of  destroyed  elements  (ADINA 
capabilities).  The  disadvantages  center  on  the 
lack  of  available  Information  on  the  material 
parameters  to  use,  numerical  Instabilities,  no 
generalized  approach  to  incorporate  empirical 
models,  and  the  appropriate  mix  of  macro-micro 
models  (reasonable  analysis  with  reasonable  cost). 
For  this  effort,  the  finite  element  results 
compared  favorably  with  the  test  data.  The 
modeling  exhibited  the  weak  elements  of  the 
structural  system  and  provided  better  pre-test 
estimates  for  the  instrumentation  system.  The 
post-test  analysis  demonstrated  failure  modes  which 
were  not  obvious  from  the  test  results.  It  is 
clear  that  an  appropriate  level  of  finite  element 
modeling  enhances  conventional  weapon  tests  against 
structures  particularly  in  the  pre-test  stages. 


ACKNOWLEDGEMENTS 

1 

COSMOS  Is  a  registered  trademark  of  the 
Structural  Research  Corporation. 

2 

AOINA  Is  a  registered  trademark  of  AOINA  RiO. 
Incorporated. 

The  authors  gratefully  acknowledge  the 
technical  support  of  Mr.  Marty  Fertal  who  conducted 
the  AOINA  simulations  described  in  this  paper. 


n 


D^UAVEO  VENTIIK:  of  ItrrEBHM.  BLASTS 


D.  U.  Kclaon.  J.  M.  UatC,  and  R.  L.  Holaaa 


U.  S,  Aiay  Enginaar  Vacarvaya  Expariaant  Sucion 
F.O.  Box  631 

Vtokaburs.  Mlaalaalppi  391S1-0631 


ABSTRACT:  VIhan  calculating  intaimal  blaac  loada 
for  chaabara  having  frangtbla  panala  covarlng  vanc 
opanlnga,  It  nay  ba  nacaaaary  to  conaldar  cha 
affact  oi  tha  panal  oonnaotlona  on  cha  vantlng 
tntClaCloft  tlM.  If  tha  panal  connacclona 
Caaporarlly  caatac  falluca,  tha  onaaC  of  vantlng 
can  ba  clgnlfleantly  dalayad.  Thla  dalay  vlll 
raaulc  In  a  gtaatar  Intacnal  loading  than  would  ba 
ealculatad  If  tha  panal  connacclona  fallad 
Inatantanaoualy.  Thla  papar  praaanta  a  aeehod, 
eouplad  with  axlaelng  analytical  cocpuccr  codoa, 
tixat  will  pradlcc  cha  claa  at  which  panal 
connacclona  fall  and  dia  caaulelng  louroal  hlaae 
load. 


URtOMICTXOH 

An  Incamal  axploalon  cyplcally  ptoducas  a 
eoaplax  load  eowpofod  of  hlgjfily  uanalant  aultlpla 
ahock  pcaaauraa  and  a  calaelvaly  long  duration  gaa 
praoiura.  In  a  eloaad  cHaabar,  thaaa  praaautaa 
can  ba  otimecurally  dattaging.  ttMtafora.  to  halp 
vaduco  Ineamal  praaauraa,  eponlngo  at#  loaaeUMO 
pKowidad  CO  ottiekly  vanc  tha  pcaaturat  out  of  tha 
chanbac.  Horaally.  thaaa  opanlnga  ata  eovacad 
with  a  frangibla  panal  chat  of  fora  nlniauil 
iNialatanea  to  blaac  praaausaa ,  Tha  atfaccivanaaa 
of  cha  vanc  opening  in  tiui  eaduccien  of  Ineanal 
praaautaa  la  dapaodant  on  cha  eaaponaa  of  Uu 
uauglbla  panal. 

Frangtbla  panala  ara  typically  llghcwalghc 
and  braak  avay  wlchout  aignlftcanely  confining  tha 
InCacnal  bloat.  Tha  panal  connscciona  aco  uaually 
daalgnad  to  fall.  KovoMnC  of  tha  panal  la 
aaavMMd  to  bagln  with  tha  arrival  of  Uta  Initial 
ahock  praaaufa.  Tho  panel  oowtwant  producaa  a 
varlabla  want  ataa  that  Incraataa  fro*  zaro  to  a 
aaxlaud  aa  tha  panal  oovaa  aWay  fron  cha  cpanlng. 
For  Um  parclcuUr  caaa  In  which  cha  panal  la 
racaaaad  within  tha  vanc  opening,  variahla  vantlng 
.doaa  nac  tnlciata  until  cha  panal  haa  novad  tha 
dlacanca  of  tha  racAea.  Tanctato  and  Htlaat* 
(Rafacanca  1)  hava  praaentad  a  Mthod  for 
calculating  cha  ineatnal  pcaaeuta  dacay  whan 
varlabla  vantlng  occura.  Including  tha  affact  of  a 
tocasaad  panal .  In  thla  awthod,  avallabla  fro* 
tha  iCaval  Civil  foglnaorlng  Laboratory  on  a 
alcroeoaipucar  progran  called  nUhC.  the  frangtbla 
panal  la  cevatad  aa  a  rigid  placa  that 
accalaraCaa  undar  tha  accloo  of  tha  Ineatnal  load. 


Panal  connacclona,  which  aro  uaually  daalgnad  to 
fall  quickly,  ara  ignorad.  FRANG  haa  boan  shown 
to  produca  good  aatiaatoa  of  inComal  blast  loads 
for  thaaa  condlclona. 

Thla  paper  addraaaaa  a  condition  In  which 
venting  la  dalayad  bacauaa  tha  frangtbla  panal 
com^actiona  aoaNncarily  raaiat  fsilura,  A  alniiar 
application,  not  sddraaaad  in  thla  papjtt,  is  one 
in  which  a  chai^ar  wall,  not  apaclfleally  daalgnad 
aa  a  vanc  area  but  known  to  ba  considerably  weaker 
chan  the  rate  of  cha  ohaabar.  nay  fall  aa  a 
frangible  panal.  In  both  caaoa.  If  an  aatliuta  of 
cha  taaecion  capacity  of  cha  panal  or  wall  can  ba 
obtaintd.  than  cha  aathod  pratanetd  in  this  papar 
would  ba  useful  for  aatlaatlng  the  intamal  load. 

A  ganaral  ouelina  of  the  Mthod  will  ba  praaancad 
first,  followed  by  an  llluaeraead  axaapla  of  tha 
Mthod  and  eeepatiaona  with  actual  case  saaulta. 

enter  or  oeureo  vektxho  om  XHTtfutAL  biast  loads 

vihan  panal  connacclona  do  not  fail 
laMdlacaly.  vaneiitg  la  dalayad.  This  delay  can 
causa  a  aigntftcane  inerasM  in  the  totairnAl 
loading  by  eo(<U!lning  the  blaac  and  allowing  cha 
gaa  ptasiurs  to  build  up  to  a  highar  laval  before 
venting  la  inleiatad.  (tote  that  this  la  a 
different  type  of  daisy  than  i^C  caused  by  a 
raceasad  panal  as  datcribad  above,  a  racaaaad 
panel  potsa  no  ealculatletial  problaai  for  the  FRaMG 
coda  bMauaa  cha  panal  ta  Still  aasuMd  to  begin 
Mving  undar  the  Lnfluanee  of  cha  Initial  ahock 
lapuUa.  For  the  case  whan  panel  connacclona 
causa  tha  dalay,  a  structural  oslcuUtlcn  la 
raqultad  to  datatwlna  the  tlM  at  which  tha 
cohnactlona  fall  and  patwl  MVaMne  begins. 
Calculating  Um  Ineatnal  load  In  Chaaa  caaaa 
Involvaa  aawaral  atapa  and  cannot  proaantly  ba 
done  using  e  elngla  coeipucer  coda.  For 
dlscusaiona  1«  this  papar,  ‘‘dalayad  venting*  will 
ba  used  to  Idantlfy  Caaaa  in  which  panel 
eonnacelons  do  not  fall  Inatancanaoualy.  Also, 
*tlM  of  fatluta*  will  ba  uaad  to  rafar  to  tha 
elM  at  whli^  thaaa  panal  conaactlons  do  fail. 

CAU»UTl«ie  tXKE  OF  FAILURE  AltS  THE  tUtOJM.  LOAD 

A  pclMry  point  in  ealcutacing  Intamal  blast 
loada  for  dalayad  venting  is  to  daearaino  tha  tlM 
of  connaoelon  fatluta.  Ones  tha  tlM  of  failuse 
la  known,  tha  praaaura  clM'hittory  can  ba  aaaily 
eonseruccad  using  axlaelng  coeputar  codas. 
Calculating  the  cLm  of  Callura  raqulras 
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consideration  of  the  alrblast  load  function  and 
tha  flexural  resistance  of  the  panel,  both  of 
Hhlch  ere  functions  of  tlae.  The  Internal  load  Is 
calculated  in  two  segments,  one  before  and  one 
after  the  panel  connection  failure.  The  two 
segments  are  then  Joined  together  to  form  the 
total  estimated  alfblarst  load  funtlon. 

The  method  presented  employs  four  computer 
programs.  In  addition  to  the  FRANC  program 
mentioned  above,  the  following  programs  are  used; 
6LAST1NW,  an  internal  blast  code  (Reference  2) , 
SDOF,  a  slngle'degree'of-freedom  dynamic 
structtiral  analysis  code  (Reference  3) ,  and 
FAILTIME,  a  code  which  Interfaces  with  SDOF  to 
calculate  the  time  of  panel  connection  failure 
(Reference  4). 

The  method  involves  five  major  steps  and 
requires  a  working  knowledge  of  single  degree  of 
freedom  equivalent  system  analytical  procedures. 

A  thorough  discussion  of  dynamic  structural 
response  and  SOOF  analysis  Is  found  In  References 
3  and  6.  In  the  steps  that  follow,  no  attempt  Is 
sade  to  list  all  of  the  assumptions  and 
limitations  of  the  four  codes  used.  This 
information  can  be  obtained  by  contacting  the 
parties  referenced  for  each  code. 

STEP  1:  CALOJUTE  INTERNAL  LOAD  FOR  CLOSED  CHAMBER 

If  It  Is  assumed  that  the  pressure  time- 
history  prior  to  the  time  of  failure  Is  Identical 
to  the  pressure  time -history  in  a  closed  chamber 
up  to  that  point  In  time,  then  the  BLASTING;  code 
can  be  used  to  quickly  construct  the  first  segment 
of  the  load  function.  The  first  step  Is  to  run 
BLASTINV  for  a  closed  chamber  Identical  to  the 
actual  chamber  with  Che  exception  that  venting  Is 
not  allowed.  This  calculation  produces  a  pressure 
time-history,  P(t),  for  a  fully  contained 
explosion  without  any  decay  associated  with 
venting.  At  this  point,  the  venting  initiation 
time  is  unknown. 

This  closed  chamber  pressure  time -history, 
F(t),  will  be  used  to  drive  the  calculations  In 
Steps  2  and  3  until  Che  time  of  failure  Is 
detemlned.  It  Is  necessary  to  Integrate  Che 
closed  chamber  pressure  time-history,  P(c),  to 
obtain  the  Impulse  clme-hlstory,  1(c),  for  use  In 
Step  4.  In  Step  4  Che  FRANC  code  will  be  used  to 
calculate  the  decaying  gas  pressure  time-history, 
Fg(C),  resulting  from  variable  venting  following 
the  time  of  failure.  Finally,  In  Step  3  the 
closed  chamber  pressure  clme-hlstory  will  be 
truncated  at  the  time  of  failure,  and  Che  decaying 
segment  of  the  pressure  clme-hlstory  from  Step  4 
will  be  appended  Co  construct  the  final  estimated 
pressure  clme-hlstory  within  Che  chamber. 

STEP  2:  CALCULATE  THE  DYNAMIC  RESPONSE  OF  THE 
PANEL 

The  panel  Is  analyzed  using  the  SDOF  code  to 
determine  Che  panel  center  displacement  Clme- 
hlstory,  D(C),  using  as  Input  the  closed  chamber 
pressure  clme-hlstory,  P(t),  of  Step  1;  the 
flexural  resistance  function  of  the  panel,  K(d); 
and  Che  mass  of  the  panel.  The  realstanco 
function,  K(d),  Is  Che  variation  In  panel  flexual 
resistance  with  displacement.  The  panel 
connections  are  assumed  infinitely  strong.  The 


SDOF  displacement  time-history,  D(C),  is 
calculated  well  past  the  time  of  expected  panel 
connection  failure.  This  displacement  time- 
history  is  used  in  Step  3  to  calculate  the 
resistance  time-history,  R(t),  of  the  panel. 

STEP  3:  CALCULATE  THE  TIME  AT  WHICH  PANEL 
CONNECTIONS  FAIL 

FAILTIME  first  constructs  a  resistance  time- 
history,  R(t),  using  Che  SDOF  displacement  time- 
history,  D(t),  from  Step  2  and  the  flexural 
repistance  function,  K(d),  of  the  panel.  Then 
FAILTIME  constructs  the  dynamic  reaction,  V{t),  by 
solving 

V(t)  -  a[P(t)]  +  b[R(C)]  (Equation  1) 

(Reference  6)  using  Che  resistance  time-history, 
R(c),  and  the  closed  chamber  pressure  time- 
history,  P(t),  from  Step  1.  In  this  equation,  "a" 
and  "b*  are  coefficients  that  vary  according  to 
whether  the  response  Is  elasti'c  or  plastic. 

Once  tha  dynamic  reaction,  V(t),  Is 
known,  Che  time  at  which  the  panel  connections 
fall  can  be  determined.  The  coiuiecclons  are 
assumed  to  fail  whan  the  dynamic  reaction,  V(t). 
at  the  ends  of  Che  panel  exceeds  the  connection 
capacity,  V(nax),  of  the  connections.  The 
connection  capacity  Is  Che  largest  dynamic 
reaction  that  can  be  sustained  by  the  panel 
connections  before  failure  occurs.  When  one  of 
several  modes  of  failure  are  possible  (for 
example,  bolt  failure  In  shear  or  tension,  anchor 
pullout,  concrete  breakout,  etc.),  the  connection 
capacity  Is  based  on  the  mode  of  failure  requiring 
the  least  reaction.  The  dynamic  reaction,  V(C), 
at  one  end  of  Che  panel  is  calculated  in  units  of 
force  (pounds) .  The  connection  capaulCy,  V(max) , 
must  also  be  expressed  in  units  of  force.  After 
entering  the  connection  capacity,  V(max) ,  FAILTIHE 
checks  to  see  if  Che  connection  capacity  was 
exceeded  by  the  dynamic  reaction,  V(c).  By 
plotting  the  dynamic  reaction,  V(t>,  with  the 
constant  connection  capacity,  V(max),  a 
decermination  can  be  made  as  to  when  and  if  the 
panel  connection  failed. 

FAILTIME  Is  written  for  simply  supported, 
uniformly  loaded,  one-way  panels.  Future 
revisions  of  the  FAILTIME  code  will  allow 
creatmenc  of  panels  having  all  four  sides 
supported. 

STEP  4;  VARIABU  VENTING  CALCULATION 

With  an  estimate  for  the  time  at  which  Che 
penal  connections  will  fall,  the  FRANC  code  can  be 
used  to  calculate  Che  gas  pressure  decay.  Fg(c), 

In  tV  chamber  as  the  panel  moves  away  from  Che 
open,.. .a  end  venting  occurs.  To  do  this  Che 
Impulse  acting  on  the  panel  at  Che  time  Che  panel 
connections  fail,  I(fail),  must  be  obtained  from 
Che  closed  chamber  impulse  clme-hlstory,  I(t), 
from  Step  1.  This  is  the  Impulse  Chat  is  used  by 
FRANC  CO  accelerate  the  panel,  (In  a  standard 
FRANC  calculation,  when  there  Is  no  delay  of 
venting  caused  by  the  panel  connections,  the 
Initial  shook  impulse  is  used  Co  acoslsrate  the 
panel).  After  additional  input.  Including  the 
equivalent  charge  weight,  the  panel  wsight,  end 
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dlaanaloncl  paramettrs  of  tha  chanber  and  vant 
araa,  FSANG  calculatas  tha  gas  prassura  claa> 
history,  Pg(t),  as  It  dacays  to  zaro. 

STEP  5;  CONSTRUCTION  OF  THE  TOTAL  PRESSURE  TIME- 
HISTORY 

As  stated  aarllar,  for  dalayad  venting  tha 
FRANC  calculation  rapraaants  only  that  sagaant  of 
tha  prassura  tiaa-hlstory  consisting  of  tha 
prassura  decay  following  the  time  of  failure.  Tha 
BLASTINtf  calculation,  F(c),  from  Stop  1,  based  on  a 
closed  ohaabar,  does  not  Include  tha  affects  of 
venting.  Tharafora,  the  final  step  Is  to  truncate 
the  closed  chaabar  prassura  tlaa-hlscory,  P(C),  at 
the  astlaatad  tlaa  of  failure  and  chan  append  tha 
FRANC  decaying  prasaura  clae-hlstory,  Pg(t),  at 
this  point  to  construct  tha  final  astlaata  of  the 
total  Incamal  blast  loading. 


COMPARISON  OF  GALCUUTIONS  VITH  TEST  RESULTS 

Tha  following  axaapla  compares  calculations 
from  the  flva-stap  aathod  of  analysis  and  actual 
test  results.  An  axploalva  charge  was  datonatad 
Inside  a  chaabar  having  a  voluae  of  approxlaataly 
2,600  cubic  feat  and  one  vent  araa  of  120  square 
feat.  Tha  opening  was  covered  with  a  panel  having 
an  approxlaata  weight  of  C.20  psl.  The  panel 
aaaaurad  10  feat  high  by  12  feet  wide  and  was 
connected  only  across  tha  cop  and  bottoa  adgaa. 

The  panel  was  considered  to  have  one-way  action  In 
flexure  and  slaple  supports.  An  18- Inch-wide 
section  of  the  panel  was  analyzed.  This  section 
had  an  astiaaCad  ccnneotlon  capacity,  V(Bax),  at 
aach  end  of  23,200  pounds  based  on  eoablnad  shear 
and  tensile  strasaes  In  the  penal  connections.  In 
the  following  discussion,  rafarsnee  to  the 
■baokvall*  will  Identify  tha  ohaabar  well  opposite 
the  vent  tree. 

Celculetlons  ware  aade,  following  tha  five 
scape  outlined  above.  First,  using  the  actual 
test  chamber  dimensions,  charge  date,  ate., 
EUStINH  oalouleted  the  evaraga  prassura  tlaa- 
history,  F(c),  on  the  baokvall  as  if  venting  had 
not  ooourrad;  that  Is,  for  a  closed  ohaabar 
(Figure  1) .  Using  this  oloaad  chaabar  pressure 
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time-history  and  the  structural  properties  of  Che 
panel,  SDOF  calculated  the  displacement  time- 
history,  D(t),  of  the  panel.  Tha  SDOF  Input  and 
resulting  displacement  time-history  are  shown  In 
Figures  2  and  3,  respectively.  FAILTIME  then 
constructed  tha  resistance  time-history,  R(t) , 
(Figure  4)  using  the  flexural  reslsranca  function, 
K(d) ,  and  the  displacement  time-history,  D(t) . 

With  tha  resistance  tlma-hlstory,  R(t),  and  tha 
closed  chamber  prassura  tlae-hlstory,  F(C), 
FAILTIME  Chan  solved  Equation  1  for  the  dynaalc 
reaction,  V(t).  In  Flgura  5,  tha  dynaalc 
reaction,  V(t),  Is  coaparad  to  tha  connection 
capacity,  V(aax),  of  the  panel.  The  calculated 
dynaalc  reaction  builds  and  eventually  exceeds  the 
connection  capacity  about  48  msec  after  detonation 
of  the  explosive.  Since  Che  connection  capacity 
was  barely  exceeded,  a  conservative  estlaate  for 
Che  Clae  of  failure  was  taken  as  35  msec. 


Load  Mais  Factor  *0.7 
Critical  Dampinq  Ratio  ■■  0.2 
Unload  Slop*  ■  30.19  psl/^nch 
Unit  maw  ■  0.000512  Ibesec^/lnch^ 
Figure  2.  Resistance  function  end  SDOF  Input 
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Figure  1.  BLASTIN  closed  chsaber  pressure  tine* 
history 


Flgura  1.  SDOF  dtsplaoeaent  tlaa-hlscory 
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Figure  4.  FAILIIME  resistance  cloe-hlseory 
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Figure  5.  Conpsrtson  of  FAILTIME  dynsnlc  reaction 
with  estlBacsd  rsooelon  oopoolty 


The  lapulee  at  the  tlsw  of  failure,  X(fatl), 
e<iual  to  1,750  pel  •usee,  wee  estluted  fro«  the 
closed  ehsuber  tupulse  eliae>hlseoty,  t(e).  In 
Figure  1.  lM)ulse  end  other  panel  end  ohsrge 

date  were  entered  Into  FRAMO  to  eelculete  the 
decay  of  the  gee  pressure,  Fg(c),  In  the  cheaber 
after  the  clue  of  failure.  This  segeenc  of  the 
pressure  elue-hlscoty  Is  shown  In  Figure  6. 
Finally,  the  closed  oheubet  pressure  clae^history, 
f<t),  woe  eruneated  ee  the  esclasted  eloe  of 
failure,  end  the  gee  preisuci!  ctits>hlscory  decay, 
rg(t),  wee  appended  (Figure  7)  to  obteln  the  final 
estlaete  of  the  total  Internal  loading. 

A  typical  recorded  pressure  and  lopulse  doe* 
history  on  the  bsckwell  froa  the  actual  test  Is 
shown  In  Figure  8.  The  gas  pressure  rise  within 
Che  ehsuber  Is  clearly  dlsoemable  as  a  long 
duration  swell  between  15  end  100  nseo  on  which 


the  highly  transient  shock  pressures  are 
super imposed.  Figure  9  is  a  comparison  of  the 
calculated  and  measured  average  Impulse  time- 
histories  on  the  backwall.  The  calculated  and 
measured  maximum  impulses  compare  well.  The 
calculated  maximum  impulse  (2,075  psi-msec) 
exceeds  the  measured  maximum  impulse  (1,830  psi- 
msec)  by  a  factor  of  only  1.13.  It  is  interesting 
to  note  that  the  difference  in  the  calculated 
maximum  Impulse  between  the  delayed  (2,075  psi- 
msec)  and  instantaneous  (750  psi-msec) 
connection  failure  is  a  factor  of  2.4.  .This 
observation  stresses  the  importance  of  carefully 
considering  the  time  of  connection  failure  when 
calculating  internal  loads. 


Figure  6.  F8AMG  gee  pressure  decay  cioa-history 
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Figure  7.  Total  esclaaced  pressure  time-hlscory 


36 


aooo  p 
1800  - 
1600  - 
1400  - 


K  ISOO 

M 

£ 


1000 
800  j- 

600  I- 

400 

aoo 
0 


TIME  -  MSEC 


Flgurt  8.  Typical  m««(ur«d  protsur*  and  lopulaa 
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CONCLUSIONS 

With  good  engineering  Judgment,  the  method 
presented  can  be  used  to  calculate  Internal  loads 
for  a  variety  of  venting  situations  not  treated 
before.  With  a  good  understanding  of  dynenlc 
structural  analyses  and  Internal  blast  phenomena, 
the  engineer  can  produce  good  approximations  for 
Internal  loads  when  the  Initiation  of  venting  Is 
dependent  on  the  delayed  failure  tine  of  the 
panel , 
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Fl|UEa  9.  Coaperleon  of  average  baokwall  Impulse 
ealoulaClon  with  test  dau 

In  addition,  high-speed  photography  shows 
that  panel  movssMnt  began  at  epproxlmatsly  30  msso, 
whUh  compares  well  with  the  oaloulated  S3  msec 
daisy  time. 

thsrsfora,  ths  flvs-etep  method  provided  a 
good  estimate  of  the  time  of  failure  end  tiw 
Internal  blast  loading  In  ths  oheabsr. 

Although  not  presented  hersln,  this  aethod 
wee  also  used  to  predict  the  time  of  failure  and 
Internal  load  for  anothar  teat  having  ths  same 
type  panel  but  a  lower  oonnsotlon  oapaolty.  The 
method  oorraotly  predicted  a  tlsM  of  failure  of  20 
maao.  To  help  put  the  effect  of  venting  delay  In 
perspective,  comparison  of  the  maximum  Impulss 
from  these  two  tests  reveals  an  approximate  40 
percent  Increase  by  delaying  the  time  of  failure 
from  20  to  SO  maeo. 
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ABSTRACT 

In  this  paper,  we  further  investigate  the  role 
of  Structure  Medium  Interaction  (SMI)  on  burled 
structure  response.  A  simple  SMI  model  can  be 
constructed  from  applying  the  boundary  condition 
between  the  soil  and  structure  that  requires  con¬ 
tinuity  of  both  stress  and  displacement  at  the 
soil  structure  interface.  Exploiting  these  simple 
boundary  conditions,  coupled  with  a  simple  model 
of  the  structural  deflection,  a  set  of  differen¬ 
tial  equations  that  represent  the  response  of  a 
burled  structure  to  explosively  produced  loading 
are  derived.  Limiting  cases  are  developed  for 
flexible,  perfectly  plastic  structures  to  fully 
demonstrate  the  interaction  between  the  incident 
ground  shock  loading  and  the  structural  defor¬ 
mation.  It  is  shown  that  deformation  is  con¬ 
trolled  by  applied  stress  and  displacement  of  the 
soil.  For  very  flexible  walls.  It  Is  shown  that 
the  wall  separates  from  the  soil  very  early  In 
time,  and  the  solution  degenerates  to  the  familiar 
Impulsive  load  response  case. 


INTRODUCTIOM 

Structure  medium  Interaction  (SMI)  effects 
play  an  Important  role  in  the  response  of  burled 
protective  structures  when  subjected  to  ground 
shock  loading  from  conventional  weapons.  The 
Influence  can  be  profound  for  highly  flexible 
structures,  with  the  effect  diminishing  as  the 
structure  becomes  rigid.  Several  papers  have 
appeared  In  the  recent  literature,  for  example, 
Drake,  Prank  and  Rochefort  (l),  Hlmiun  and 
Weldllngei*  (2)  and  (3),  that  have  Incorporated  SMI 
models  in  the  design  of  burled  structures. 

While  the  application  of  SMI  models  has  only 
recently  been  Introduced  Into  the  design  of  pro¬ 
tective  structures,  the  role  of  SMI  on  the 
response  of  embedded  structures  has  bean 
recognized  for  more  than  a  century.  For  example, 
Rayleigh  (4)  studied  the  radiation  field  produced 
by  plana  sound  waves  Impinging  on  rigid  and 
gaseous  spherical  Inclusions.  The  designs  of 
naval  structures  (l.e.,  submarines  and  ships)  have 
used  the  Taylor  plate  models  and  plane  wave 
interaction  models  sines  World  War  ll.  Perhaps 
the  most  widely  used  formulation  Is  the  OoubTy 
Asymptotic  Approximation  (OAA)  developed  by  Geers 
(S)  that  embodies  both  the  high  frequency  “plane 


wave  approximation"  and  the  low  frequency  effects 
where  the  motion  of  the  structure  and  fluid  are 
acting  Inphase. 

It  was  demonstrated  by  Drake,  et  a1.  (1),  that 
a  simple  RMI  model  can  be  constructed  from 
applying  the  boundary  condition  between  the  soil 
and  structure  that  requires  continuity  of  both 
stress  and  displacement  at  the  soil  structure 
Interface.  Exploiting  these  simple  boundary  con¬ 
ditions,  coupled  with  a  simple  model  of  the  struc¬ 
tural  deflection,  a  set  of  differential  equations 
that  represent  the  response  of  a  burled  structure 
to  explosively  produced  loading  were  developed. 

The  equations  have  direct  analogy  to  the  familiar 
single  degree-of-freedom  (SDOF)  system  with 
damping  and  are  easily  solved  by  analytical  and 
numerical  means.  A  comparison  of  mid-span  deflec¬ 
tions  meosured  on  burled  wall  explosive  tests  and 
computations  using  this  simple  model  was  within 
±30  percent  of  experiments. 

In  this  paper,  we  further  Investigate  the  role 
of  SMI  on  burled  structure  response.  Limiting 
cases  are  developed  for  flexible,  perfectly 
plastic  structures  to  fully  demonstrate  the 
Interaction  between  the  Incident  ground  shock 
loading  and  the  structural  deformation.  It  Is 
shown  that,  to  first  order,  deformation  is 
controlled  by  both  the  Incident  stress  and  the 
displacement  of  the  soil.  For  very  flexible 
walls.  It  is  shown  that  the  wall  separates  from 
the  soil  very  early  In  time,  and  the  solution 
degenerates  to  the  familiar  impulsive  load 
response  case. 

PROBLEM  DEFINITION 

A  formulation  of  the  SMI  model  Is  given  by 
Drake,  at  al.  (1)  by  a  combination  of  simple  wave 
propagation  theory  in  the  soil  and  rigid  body 
mechanisms  In  the  structuM.  Boundary  conditions 
of  continuity  of  both  stress  and  displacement  be¬ 
tween  the  soil  and  structure  are  Imposed  for  a  SMI 
model.  At  the  interface, 

01  -  Off  4-  Or  (1) 

whore  o^  1$  the  Interface  stress,  Off  is  the 
Incident  free-flald  stress  and  Or  Is  the 
reflected  stress  from  the  structure.  Also  at  the 
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boundary,  continuity  of  displacements  requires 
that 


Vff  -  Vp  -  u  (2) 

where  Vff  Is  the  particle  velocity  associated 
with  Off  ,  Vp  Is  the  reflected  particle  velocity 
and  u  Is  the  velocity  of  the  structure.  The 
equation  of  motion  for  the  structural  motion  is 


0^  ■  psLJl  +  R(u) 


(3) 


where  p*  Is  the  mass  density  of  the  structure,  L 
Is  the  structure  thickness,  and  R(u]  is  the 
resistance  per  unit  area.  For  an  elastlc- 
perfectly  plastic  structure 


R(u)  - 


u  <  Up 

U  >  Up 


(4) 


Thus  the  maximum  velocity  that  the  structure  can 
obtain  Is  twice  the  free-fleld  particle  velocity 
incident  to  the  structure.  It  can  also  be  seen 
from  £q.  3  that  the  Interface  stress  Is  reduced  to 
the  resistance  at  the  time  that  the  maximum  veloc¬ 
ity  Is  reached. 

The  role  of  the  properties  of  the  backfill 
material  can  be  Investigated  parametrically  by 
exploring  the  response  of  a  perfectly  plastic 
system.  This  solution  Is  easy  to  obtain  and  can 
be  used  to  Illustrate  the  salient  features  of  the 
SHI  problem.  Detailed  Integration  of  Eq.  3  with 
complex  resistance  functions  are  in  good  agreement 
with  the  simple  bounding  solutions. 

SOLUTIONS  FOR  PERFECTLY  PLASTIC  RESPONSE 

Consider  a  perfectly  plastic  structure  with  a 
constant  resistance  function. 


Using  the  relationship 

Op  ■  pcVp  (5) 

where  p  1s  the  mass  density  and  c  Is  the  prop¬ 
agation  velocity  of  the  soil  and  incorporating 
the  boundary  conditions,  Eqs.  1,  2,  and  4,  results 
In  the  expression 

PjLU  +  pcO  +  R(u)  •  Off  +  pcVff  (6) 

which  1$  the  equation  of  motion  for  the  structure 
that  includes  the  Interaction  effect.  Strictly 
speaking,  Eq.  5  Is  for  a  linear  soil,  but  for  real 
soils,  the  stress  and  velocity  time  histories  have 
different  waveforms  during  unloading.  Note  that 
the  SMI  effect  manifests  Itself  as  a  damping  term 
related  to  the  radiation  damping  provided  by  the 
reflected  wave  from  the  structure.  Also  note  that 

at  shock  Impingement  on  the  structure.  Thus  the 
Initial  reflection  factor  Is  two  times  the  Inci¬ 
dent  stress  which  arises  from  the  tacit  assumption 
that  the  structure  Is  moving  as  a  rigid  body.  As 
pointed  out  in  Reference  l,  substituting  a  reflec¬ 
tion  factor  based  on  acoustic  wavs  propagation 
theory  In  place  of  the  factor  of  two  will  not  con¬ 
serve  momentum  in  this  model.  Also,  the  reflec¬ 
tion  factor  produced  by  this  model  can  be  as  low 
as  one  —  depending  on  the  rise  time  of  the  Inci¬ 
dent  stress  pulse. 


^  "  l^x 


Eq.  6  Is  easily  integrated  to  give 


b 

/[Off  +  PcVff  -  Rflax]  dt  (9) 


where  n  ■  pc/m  and  m  >  pjL  . 


The  free-fleld  stress  and  velocity  time 
histories  can  be  estimated  as  (for  example,  see 
Reference  6) 


Off  •  pcVoe’<*^ 
Vff  .  Vos-Bt 


(10) 


where  a  ■  r/c  ,  J  -  1/2.5  a  and  r  is  the 
distance  to  the  structure  from  the  explosion. 

Note  that  a  simpler  form  of  the  velocity  pulse  Is 
assumed  here  to  facilitate  the  integration. 

Hence, 


v,  •  riii  *  rk 

-  Saa*  (i-e-nt)  (u) 

omax 


where  .  peVo  . 
Limiting  Case  n  »  a  >  B 


Koto  that  during  the  early  phase  of  the 
response  when  R(u)«^  ,  then  Eq.  3  is  •>  psLtl  , 
showing  the  Interface  stress  Is  simply  the  accel¬ 
eration  Of  the  structure  times  the  mass.  However, 
the  stress  quickly  decays  as  the  structure  Is 
rapidly  accelerated  to  the  particle  velocity  of 
the  soil  at  the  Interface.  The  maximum  structural 
velocity  occurs  whan  (1  >  0  , 

Oflax  •  2  Vff  -  R/pc  (8) 


For  most  flexible  structures,  n  »  a  >  8  • 
That  Is,  the  response  time  of  the  structure  mass 
Is  much  lass  than  the  duration  of  the  free-fleld 
stress  duration.  For  this  case,  as  t  >  l/n  , 
e*ht  0  ,  the  wall  velocity  becomes 


L  m  e-«t  +  (-8t  . 

Vo  %ax 
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From  this  expression  It  can  clearly  be  seen  that 
the  structure  velocity  simply  follows  the  velocity 
history  of  the  free-fleld  stress  and  particle 
velocity  components. 

The  displacement  of  the  structure  can  be 
easily  obtained  by  integrating  Eq.  12  as 


u  «  Vo  (l-e‘«t)  +  Vo  (i-e'S*’-)  -  ^max  t  (13) 
a  F"  pc 

Note  that  the  free-fleld  displacement  is 


Uff  •  Vq  (i,e‘6t) 


and  that  a  ~  2.50  .  so  that  ths  displacement  of 
the  structure  Is  following  the  displacement  of  the 
soil,  to  first  order. 

Thus  for  the  case  of  Rmax/pc  «  Vq  .  the  maxi¬ 
mum  displacement  of  the  structure  Is  clearly 
bounded  by 

Umax  s  Uff  (1+e)  (14) 

where  e  ■  0/o  <  1  .  Of  course  for  very  resistant 
structures  the  displacements  fall  below  this 
bound,  provided  that  the  resistance  Is  sufficient 
to  maintain  contact  between  the  soil  and  structure 
at  the  Interface. 

Drake,  et  al.  (l)  demonstrated  this  by  para¬ 
metric  calculations  and  plotted  the  results  as 
shown  In  figure  i  for  0  >  d  .  It  can  be  seen 
that  ths  maximum  deflection  is  proportional  to  the 
free-fleld  displacements  over  a  wide  range  of 
Input  conditions.  The  dashed  lines  indicate  where 
the  structure  separates  from  the  soil,  as  dis¬ 
cussed  In  the  next  section. 


FUE-rtCLO  STUSS/UStSTMia 

Figure  1.  Maximum  deflection  of  a  perfectly 
plastic  SDOF  with  SHI  model  from 
Drake,  et  al.  (l). 


Thus,  we  note  that  the  structural  sesponse  for 
flexible  structures  Is  dominated  by  ths  motion  of 
the  free-fleld  Incident  to  the  structure.  To 
first  order,  the  structure  velocity  quickly 
obtains  and  maintains  the  particle  velocity  of  the 
soil.  The  effect  of  the  structure  resistance  Is 
to  reduce  this  velocity  by  a  constant  amount. 

From  the  free-fleld  ground  shock  equations,  It 
Is  noted  that  at  a  constant  standoff  distance, 

Vq  »  constant  ^ 

Offlax  a  pc  >  (15) 

Uff  a  1/c  / 

Thus  a  tradeoff  In  the  design  and  backfill  proper¬ 
ties  Is  suggested.  High  quality  backfill 
materials  will  exhibit  higher  c  values,  resulting 
In  lower  free-fleld  displacuments,  but  higher 
Incident  stresses.  As  shown  In  figure  I,  varying 
the  resistance  provides  only  a  modest  decrease  In 
the  peak  deflection. 

Limiting  Case  for  Cavitation  at  Interface 

For  very  flexible  structures  (L/0  »  10),  the 
structure  mass  Is  very  quickly  accelr-ated  to 
velocities  approaching  twice  the  fre:  leld  par¬ 
ticle  velocity.  In  these  cases,  tension  can  devel¬ 
op  at  the  soil -structure  interface,  resulting  In 
the  structure  being  spalled  or  separated  from  the 
soil.  Rejoining  with  the  soil  may  occur  later  in 
time,  depending  ort  the  aecay  characteristics  of 
the  Incident  pu;se  and  the  structure  resistance. 

It  can  be  shown  from  the  solutions  obtained  In  Eq. 
9-U,  that 

Oi  ■  2  Off  e’Dt 

for  very  early  times.  Therefore,  from  Eq.  3, 

ml!  •»  2  pcVfle"*!^  -  Rd]jx 


and 


0  -  2  Vo  -  Rmax  t/m  (16) 

w'>ere  m  ■  Ps*-  •  displacement  occurs  In 

t>1$  case  ac  a  time,  t  ■  2  Vom/Ro)ax  •  S^ves 

Umax  “  ^  Vg  (17) 

which  Is  the  familiar  impulse  response  of  plastic 
systems. 

The  expresslin  in  Eq.  8  can  be  Improved  some¬ 
what  by  using  the  peak  velocity  of  the  wall  as 
given  by  Eq.  8,  and  estimating  the  time  of  peak 
'  velocity  as 


Wx  • 


ilfmaxiiSnajt 

(pc/m  *  a) 


(IB) 
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which  gives 


Ofliax  s  2-Rinax/gflvix 
1  +  aniypc 


(19) 


Thus  a  better  estimate  of  the  peak  displacement  is 

(20) 


2  Rmax  \  1  +  ow/pc  / 


Criteria  for  the  cavitation  at  the  interface 
can  be  established  from  Eq.  12.  The  Interface 
stress  can  be  calculated  from  Eq.  3  noting  that 


U  ■  -ae"®^  -0e“9^ 

^0 

resulting  In  the  expression 

•  d<  *  R  -  mVo  (ae"®t  +  Be'^t)  (21) 

In  order  to  remain  in  contact  with  the  soil, 

0^  >  0  ,  which  requires 


Ssas  i  Vo  0  (a-«t  +  e/a  e-fft)  (22) 

m 

For  purposes  of  developing  an  estimate,  note  that 
d  •  c/r  and  6  «  q/2.5  so  that 


Figure  2.  Contparlson  of  calculated  and  observed 
rotation  for  burled  wall  tests  from 
Drake,  et  a1,  (1). 


>  Vo£  or  Ssas  2  fei:  (23) 

w  f  ®max  P'' 

Thus  to  avoid  cavitation  and  minimize  the  struc¬ 
ture  deflection,  the  ratio  of  the  resistance  to 
the  incident  stress  mutt  be  greater  than  the  ratio 
of  mass  of  the  structure  to  the  mass  of  the  soil 
between  the  structure  and  the  explosion. 

Therefore,  In  the  case  of  very  flexible  struc¬ 
tures  where  cavitation  is  the  dominant  response 
mode,  the  maximum  deflection  Is  largely  Independent 
of  the  soil  medium,  increasing  somewhat  for  high 
velocity  soils  as  the  term  R/pc  becomes  small, 

COMPARISON  WITH  TESTS  AND  CALCUUITIONS 

A  number  of  burled  structure  tests  were  con¬ 
ducted  by  the  U,  S.  Army  Engineer  Waterways 
Experiment  station  in  support  of  the  U.  S.  Air 
Force  Engineering  and  Services  Center  (7),  The 
tests,  as  rsported  by  Baylot,  used  slabs  with 
length  to  thickness  ratios  of  5  and  10,  and  varied 
the  weapon  position  to  sustain  different  damage 
levels.  Calculated  wall  deflections  using  the  SNl 
model  described  in  this  paper  were  shown  to  be 
within  i30  percent  of  the  observed  deflections 
(see  figure  2).  All  of  the  salient  features 
observed  In  the  tests  were  successfully  predicted 


by  the  theory;  namely,  (a)  the  rapid  decay  of  the 
Interface  stress  following  the  Initial  peak,  (b) 
peak  Interface  stress  that  was  only  slightly 
higher  than  the  Incident  peak  stress  due  to  the 
finite  risetime  of  the  Input  pulse,  (c)  in  some 
cases,  separation  (l.a,,  ■  0)  at  the  soil 

structure  Interface,  (d)  late  time  interface 
stresses  that  approach  the  structural  resistance, 
(e)  and,  most  Importantly,  the  deflections  were 
accurately  calculated  by  a  simple  handbook  Ideali¬ 
zation  of  the  structure. 

In  these  tests,  the  slabs  with  length  to 
thickness  ratios  of  10  ware  observed  to  separate 
from  the  soil  shortly  after  the  shock  arrival. 

Peak  rotations  at  the  support  calculated  by  Eq.  20 
Is  shown  In  figure  3  compared  with  the  observed 
rotations.  It  can  be  seen  that  the  calculations 
compare  favorably  with  the  test  results  for  the 
charge  parallel  to  the  wall.  Departure  from  the 
theory  can  be  explained  on  the  variation  In  tho 
weapon  effects  due  to  weapon  orientation.  Lower 
defonnatlons  were  observed  for  weapons  orientated 
In  an  end-on  configuration  and  larger  for  a  ver¬ 
tical  orientation. 

Another  case  that  of  cavitation  dominating  the 
structural  response  is  clearly  shown  by  a  para¬ 
metric  study  by  Weldlinger  and  Hinman  (3).  Their 
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equation  20. 

example  is  that  of  a  S3-1nch-th1ck  slab  with  a 
span  of  42  feet  and  resistar.oe  of  Rn^v  >  45  psi 
subjected  to  a  partially  coupled  (f»0.4)  weapon  at 
the  ground  surface.  The  weapon  yield  was  1014  lbs 
of  TNT.  while  the  backfill  soil  was  sand,  with  a 
pc  ■  22  psi/fps,  c  ■  lOOQ  fps  and  n  ■  2.75  . 

For  the  problem  considered  the  incident  stress  was 
about  185  psi.  which  resulted  in  a  maximum  com> 
puted  displacement  of  2.9  Inches. 

Using  Eq.  20,  we  obtain  Umax  «  2.4  in.  which 
agrees  well  with  the  more  exact  solution. 

Both  SDOF  and  finite  element  (PE)  model  solu* 
tions  were  developed  In  Reference  3  for  this 
structure  In  several  soil  types,  as  shown  in 
figure  4.  From  the  bounds  developed  In  Eq.  23, 


!s§£  >  0.315 
Omax 


in  order  to  prevent  the  slab  from  separating  from 
the  soil.  In  the  baseline  case,  Rmax^^x  * 
45/185  ■  0.24  ,  so  this  clearly  falls  the  test. 

As  the  soil  Is  varied,  the  maximum  particle  ve1oc> 
ity  at  a  given  range  Is  nearly  constant  while 
the  maximum  stress  increases  proportional  to  pc  . 
Therefore  for  soils  with  c  >  750  fps  ,  the  Inci¬ 
dent  stress  can  always  cause  cavitation,  and  the 
bounding  solution  should  provide  good  results. 


Peak  displacements  estimated  from  5q.  20  are 
also  shown  In  comparison  to  the  SOOF  and  FE  com¬ 
putations  In  figure  2.  Since  we  are  not  certain 
about  the  detaiTs  of  how  the  material  was  modeled 
In  (3)  for  the  variation  In  c  ,  we  provide  two 
curves,  one  for  n  ■  2.75  and  n  ■  2.5,  which  bound 
the  computations.  The  simple,  impulse  load  esti¬ 
mate  Is  very  close  to  the  PE  calculation  for  the  n 

•  2.75  curve  up  to  c  *  3000  fps.  For  higher 
values  of  c  ,  it  appears  that  the  attenuation 
coefficient  was  varied  to  a  value  of  n  ■  2.5  at  c 

•  5000  fps. 


SUMMARY 

A  SMI  model  Is  proposed  that  accurately  models 
the  Interface  condition  between  a  struc>.ure  and 
soli.  During  the  loading  phase,  continuity  of 
both  stress  and  displacement  between  tha  soil  and 
structure  was  maintained.,  Resulting  equations  of 
motion  for  this  system  resenfble  those  of  a  SDOF 
system  with  damping.  However,  the  damping  term  Is 
the  result  of  satisfying  the  stated  boundary  con¬ 
ditions  and  Is  not  related  1n  any  way  to  viscous 
or  frictional  damping  effects. 

It  was  shown  that  Interface  stress  is  highly 
dependent  upon  the  Inertial  effects  and  the 
deflection  of  the  structural  section.  For  many 
Important  cases,  the  interface  stress  approaches 
zero  within  twenty-five  transit  times  through  the 
structural  section  and  then  slowly  approaches  the 
resistance  at  late  time.  It  was  shown  that  the 
maximum  structural  velocity  cannot  exceed  twice 
tne  Incident  free-fleld  particle  velocity  even  for 
very  low  resistance  structures.  In  contrast, 
current  design  methods  which  apply  the  free-fleld 
stresses  directly  as  the  interface  stress  violate 
the  basic  displacement  boundary  condition  and 
cause  the  structure  element  to  be  accelerated  to 
velocities  far  greater  than  physically  possible. 

A  comparison  of  theoretical  results  with 
burled  wall  experiments  Is  excellent.  Calculated 
wall  deflections  were  within  4-30  percent  of  those 
observed  In  tests.  Interface^stresses  were 
accurately  predicted  as  well.  The  theory  Is 
easily  solved  for  any  structure  and  Incident 
loading  by  analytical  methods  or  numerically  on 
desktop  or  programmable  calculators. 


Figure  4.  Maximum  deflection  computed  by  Eq.  20 
compared  with  SOOF  and  FE  analysis  of 
a  burled  wall,  from  Weldlinger  and 
Hlnman  (3). 
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ABSTRACT t 

Soil-Structura  Intacaction  ia  an  iaportant  conaidaration  in  calculating 
eha  raaponta  ot  buriad  atructuraa  to  blaat  af facta.  Raaulca  froa  tha  Shallow- 
Buriad  ScrucCuraa  (SBSJ  raaearch  prograa,  which  ia  aponaorad  by  tha  Oafanaa 
Nuclaar  Agancy,  hava  ahown  that  typical  buriad  cograand  and  control  typa 
atructuraa  can  aurviwa  aa  auch  aa  tanfold  graatar  paalc  ovarpraaaura  froa  a 
nuclaar  waapon  than  waa  thought  poaalbla.  Tha  unaxpactad  hardnaaa  of  tha 
buriad  atructuraa  raaultad  prlaarily  bacauaa  tha  affaeta  of  ioil-»tructura 
Intaraction  had  baan  undaraatiaatad.  Data  froa  tha  SBS  rtaaareh  hava  raaultad 
in  aignificant  ravialona  in  vulnarability  coaputational  aathoda  and  in  our 
aatiaataa  of  buriad  atructura  vulnarability  to  nuclaar  waapona. 

In  attaapta  to  axtrapolata  tha  analytical  aathoda  davelopad  in  tha  SBS 
rataaroh  progcaa  to  coaputa  buriad  atructural  raaponaa  to  oonvantional 
waapona,  eharo  ia  a  graat  daal  of  uncartainty  bacouaa  of  tha  localirad  loading 
and  raaponta  froa  conventional  waapona.  For  axaapla,  tha  cutvaa  ahown  in 
Pigura  1  indicata  about  tha  aaaa  pradiotad  atructural  raaponaa  in  a  backfill 
with  low  aaiaaic  velocity  (about  I, 030  fpa),  but  ahow  a  conaidartbla 
diffaranea  In  pradictad  atructural  raaponaa  at  highar  aaiaaic  valocitias. 

Thara  ara  aany  atructuraa  in  backfilla  with  highar  aaiaaic  valooltiaa, 
aapaoially  in  Europa.  Alao,  if  backfill  cyp«  doaa  not  affact  atructural 


•  Froa  lllnaan,  Eva  E.,  and  Uaidlingar,  P,  Procaadina  froa  Intarnationa! 
Syapoalua  on  Intaraction  of  Non-Nuclaar  Muaitiona  with  Structuraa,  Hannhoia, 
Watt  Oaraany,  1987. 
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r*spon*t,  a  eonaidarabla  coac  savings  could  be  realized  in  Cucure  conscruccion 
projects  by  not  requiring  a  select  backfill. 

Unfortunately,  there  are  very  little  data  fron  conventional  weapon  tests 
on  buried  structures  in  backfill  aaterials  with  seismic  velocities  outside  the 
range  of  about  800  to  1,300  fps.  The  diffarencas  indicated  by  the  curves  in 
Figure  1  have  recently  assumed  a  greater  importance,  since  the  analytical 
methods  fron  the  Army  Technical  Hanual  S-85S-1  predict  increasing  structural 
response  with  increasing  seismic  velocity,  and  analytical  methods  in  a  new 
draft  Air  Force  design  manual  predict  structural  response  to  be  very  nearly 
Independent  of  seismic  velocity. 

This  paper  will  review  the  applicable  theory,  evaluate  the  analytical 
methods  depicted  by  the  curves  in  Figure  1,  and  compare  predictiona  using 
these  methods  with  the  limited  data  base  available. 
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BEINfOaCEO  A2fl>  FSECAST  COSNECXIONS  UNDER  CIjOSE-IN  EXPLOSIONS 


T,  KrautluaiMr  and  M.  D«Suct«c 

D«p«xtMnt  of  Civil  «nd  Ninoral  Englnoorlng 
Ualvoraiey  of  Hinnoaoea.  U.S.A. 


Abatract 

Tha  caaponao  of  ratjif^Toad  and  piacaat 
OQisMcelona  undac  cha  application  of  clota'ln 
datonatlon  af facta  haa  baan  atudlad.  Typical 
connection  detail a  vata  analyzed  by  tha  finite 
alaaant  nachoda,  and  thalr  patforaanca  under 
aevara  ahotc--dutacloo  localised  loada  waa 
evaluated.  Ofaaarvad  iaflclanclaa  In  tha  behavior 
vara  axaalned.  and  aubaaquanc  daaign 
awdlfltatlona  etere  ahoun  to  correct  the  teaponae. 
Tha  procedure  for  avaluatln*  aeveral  knee-  and 
cao*Jolnta  la  peeaented  and  diaeuatad. 
teeoMMndatlona  are  provided  for  the  daalfo  of 
such  coonactlona. 


tfnually,  eonnactiona  ar«  aasuBcd  to  Imi 
rigid  which  •Hows  ths  adf  ining 

•l«Mnts  to  dAVelop  their  full  potential. 
ilow«v«r»  if  th«  oonnactlons  at*  not 
daaignad  adaquataly  that  will  not  ba  tha 
caaa*  and  tha  struoturna  will  fall 
ptaaaturaly.  £aanplaa  of  thia  type  of 
bahavidr  ava  tnafty  aa  dtaeuaaad  in 
(IctauthaiHaar.  19S7)>  whara  it  via  ahown 
that  floof'^arch  joints  axhibitad  aejor 
dahaga.  oaapita  tha  fact  that  thaaa 
•ttueturaa  wara  taatad  undar  a  aiaulatad 
nuclaar  aiwironaant,  aiailar  aodas  of 
t>ahavlor  could  ariat  undar  loealixad 
datonations.  Racant  atudiaa  on  tha 
aodalihg  Of  concrata  hohavior  undar  high 
rata  load  affacts  at  tha  Univaralty  of 
Ninnasdta  have  aada  It  pohhi^id  to 
accurataly  analyte  connacticna  aubjaotad  to 
this  type  of  loading. 


i<gfltt,DgflniU9n 

Tha  sagnituda  and  tlaa  pulse  history  of  tha 
blaat  load  ware  datarsinad  through 
interpolation  of  data  ahown  in  Coltharpt 
at.  al.(l985).  A  20"  wide  atrip  was 
suparisposed  on  tha  propoaad  prasuura  grid 
at  tha  worst  casa  location.  Tha  prassura 
waa  than  convartad  to  a  oonoantratad  load 
which  will  act  on  tha  and  of  a  oantilavar 
baaa.  Tha  resultant  peak  load  was  sst  at 
2451  kipa.  Basad  qn  tha  saaa  report,  the 
rise  tine  to  first  psak  waa  aat  at  0.04 
Mac.  Tl.a  duration  of  poaitiva  praaaura 
waa  sat  at  0.53  ayaao. 


Tan  connection  dataila  wsra  anaXyssd  for 
thair  structural  integrity  undar  blaat 
loading. 

Tha  analyais  was  coaplatad  with  a  finita 
aianant  coda  SAHSOlt  2  written  by  Schtayar, 
at.  al.(l984).  thia  coda  utilizaa  a 
conatiuativa  concrata  danaga  nodal  written 
by  Stsvana  and  Krauthaimar  (I9t6)  at  tha 
dnivarnity  of  Mlnnaaota.  in  tha 
conprasaion  dosaln  tha  nodal  utilisaa  a 
conbinatlon  of  plaatielty  and  contiuun 
danaga  nechanics  thaoriaa  with  e  nonlocal 
dafinition  Of  a  acalar  danaga  variable  to 
nodal  strain  softening.  In  tansila  regions 
a  nonaaaociatad  flow  rule  with  associated 
nodification  la  used.  Thia  procadura  will 
accurately  predict  strength  atatss  for 
tenaiofi  and  shear  while  holding  dllatlonal 
straina  to  a  nlninun  and  alininsting  tha 
posalbility  of  apontanaoua  anargy 
generation.  Tha  concrata  was  nodalad  with 
a  naab  coaposad  of  guadrstic  g  noda 
rectangular  and  6  noda  triangular  aleSants. 
Thu  sash  alsnant  aiza  waa  refined  until  a 
change  in  etrass  of  s%  or  less  waa  observed 
between  adjacent  alananta.  Thia  sash  was 
thsn  coiparad  to  a  further  ratinad  nash 
with  a  strass  daviation  of  5%  or  lass 
racordad  at  identical  points.  Tha  steal 
cainforcaaant  waa  a<ktalad  using  bar 
alananta  connactnd  to  tha  nldlina  nodas  of 
tha  concrata  alesanta. 
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The  adequacy  of  this  procedure  was  checked 
by  aodeling  reinforced  concrete  beans  and 
coaparing  the  results  to  those  measured 
experimentally  by  Feldnan  and  Siess  (1958) . 

BEAM 

MAXv  DISPLACEMENT 

(exper)  (samson2)  C-1 

3. O'*  3.2"  H-1 

8.4"  7.2"  These  results  are  very 

good  considering  the  high  degree  of 

nonlinearity  and  erratic  impulsive  loading. 


Below  is  a  brief  summary  of  each  connection 
and  it's  reinforcement  details.  All 

connections  represent  a  slice’  (or 
thickness)  of  20",  reinforcement  yield 
strength  of  60  ksi,  and  concrete  strength 
of  6  ksi  as  shown  in  Figure  1. 

Monolithic  Knee  Joint 

Five  different  geometries  with  varying 
reinforcement  were  tested. 


The  approach  used  to  model  beam/column 
connection  behavior  is  similar  to  that  used 
by  Nilsson  (1973) .  The  column  is  rigidly 
attached  at  its  base  (mid-height  between 
stories)  and  the  beam  is  cantilevered 
approximately  one  third  of  a  typical  span 
length.  This  method  was  chosen  in  order  to 
predict  a  worst  case  scenario  (i.e. 
beam/slab  discontinuity  due  to  openings  or 
localized  failure) .  A  concentrated  blast 
load  similar  to  that  outlined  above  was 
than  applied  to  the  end  of  the  cantilever. 
The  details  represent  a  beam/colinnn 
configuration.  However,  the  actual 
simulation  models  a  wall/slab  structure  by 
analyzing  a  slice  through  the  section 
subjected  to  worst  case  loadings.  This 
approach  does  not  consider  3-dimensional 
load  distribution  in  the  slab  for  two 
reasons.  The  blast  load  under 
consideration  is  somewhat  uniform  over  the 
width  of  the  stiructure  and  the 
susceptibility  to  repeated  loading  is 
likely. 

•Three  different  connection  types  with 
varying  details  were  chosen  by  their 
performance  in  previous  research.  These 
connections  were  then  tes';ed  by  the 
procedure  outlined  above  for  their 
structural  integrity.  Minor  alterations  in 
overall  dimensions,  area  and  location  of 
steel  reinforcing,  and  material  properties 
were  made  until  safe  and  efficient  behavior 
was  observed. 

The  parameters  used  to  determine  adequate 
connection  behavior  are  as  follows. 
Deflections  should  not  exceed  those 
specified  in  the  ACI  Building  Code  (318- 
33).  Direct  shear  should  not  surpass  the 
Hawkins  shear  limit  as  shown  in  Murtha  and 
Holland  (1982).  Flexural  shear  and 
compressive  stress  in  the  concrete  should 
not  exceed  ACI  specified  limits. 
Reinforcement  stresses  should  not  exceed 
1/3  yield  in  the  connection  region  and  2/3 
yield  at  the  hinge  location.  The  location 
of  the  hinge  should  occur  at  a  distance  d 
from  the  face  of  the  joint.  Following 
these  criteria  will  ensure  adequate 
behavior  under  repetitive  blast  loadings. 
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PRINARY  LONGITIOM 
REIHFCROENT 

STIRRUPS 

(CONFINEMENT  i  Sf£Afi) 

SECONOARY  REINFORCBeiT 
(SLA8-MI1 LQCATKJIS) 


TYPICAL  BEAM/'COLUHN  CROSS-SECTION 
Figure  l  ( SLAB-I’ALL ) 


The  first  connection  (BL-IA,  Figure  2) 
incorporates  an  18"  wide  by  20"  thick 
column  and  bean  detail.  The  beam  and 
column  are  reinforced  with  1%  steel  (3-19 
bars,  fy  »  60  ksi)  on  each  face.  Diagonal 
steel  (2  -  #8  bars)  is  placed  on  the  inside 
corner  of  the  beam-coluam  interface  at  a  45 
degree  angle  with  2"  of  cover.  The  #9  main 
flexural  bars  are  anchored  in  the  joint 
with  standard  90  degree  hooks.  Transverse 
stirrups  are  tied  around  the  #8  diagonal 
bars  and  the  outer  face  flexural  bars.  The 
stirrups  provide  confinement  for  the  joint 
core,  compression  resistance  for  the 
outside  flex)iral  bars,  crack  control  at  the 
interior  of  the  joint,  and  buckling 
resistance  for  the  diagonal  steel  under  a 
closing  moment.  The  size  and  spacing  of  the 
stirrups  can  bo  determined  by  the  procedure 
proposed  by  Park  and  Paulay  (1975). 

The  second  connection  (BSL2A,  Figure  3) 
incorporates  a  12"  wide  by  20"  thick  column 
and  beam  detail.  The  beam  and  column  are 
reinforced  with  1.5%  steel  (3  -  #9  bars,  fy 
■  60  ksi)  on  each  face.  Diagonal  stael  (3 
-  |8  bars)  is  placed  2"  from  the  exterior 
face  of  an  8"  diagonal  concrete  strut 
located  on  the  inside  corner  of  the  beam- 
column  interface.  The  #9  main  flexural 
bars  are  anchored  in  the  joint  with 
standard  90  degree  hooks.  Transverse 
stirrups  are  tied  around  the  18  diagonal 
bars  and  the  outer  face  flexural  bars. 
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Th«  third  connection  (Bh-fiA,  Figure  3) 
Incorporates  an  18"  wide  by  20"  thick 
colunn  and  beam  detail.  The  beam  and 
column  are  reinforced  with  0.6%  steel  (3  - 
#7  barSf  fy  *  60  kei)  on  each  face. 
Diagonal  steel  (2  >  *6  bars)  is  placed  2" 
from  the  exterior  face  of  an  8"  diagonal 
concrete  strut  located  on  the  inside  corner 
of  the  beam-column  interface.  The  |7  main 
flexural  bars  are  anchored  in  the  joint 
with  standard  90  degree  hooks.  Transverse 
stirrups  are  tied  around  the  #6  diagonal 
bars  and  the  outer  face  flexural  bars. 


Figure  3  connections  BSUA  and  BL-6A 


The  fourth  connection  (BL-6C)  is  identical 
to  BX1-6A  with  one  exception.  The  diagonal 
bar  steel  is  increased  to  3-#7s  which  is 
equal  to  the  amount  of  steel  used  for  the 
main  flexural  bars. 

The  fifth  connection  (BLS6A)  incorporates 
an  18"  wide  by  20"  thick  column  and  beam 
detail  with  sifcon  material  added  in  the 
connection  region  as  shown  in  Figure  4. 
The  material  properties  used  for  sifcon 
with  12%  by  volume  deformed  steel  wires  are 
as  follows:  modulus  of  elasticity  1200  ksi; 
compressive  strength  12  ksi;  and  a  tensile 
strength  of  1.8  ksi  (Homrich  and  Naaman, 
1988)  .  The  geometry  and  reinforcement  used 
in  this  detail  are  identical  to  the  BI1-6A 
connection. 
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Figure  4  Connection  BLSGA 


Preoaet  Knee  Joint 


TWO  types  of  precast  details  ware  tested. 
The  first  detail  (PL- ID,  Figure  S) 
inoorporatss  an  18"  by  20"  thick  beam  and 
column.  The  beam  is  plS'Sed  on  a  corbel 
which  protrudes  6"  from  the  face  of  the 
column.  A  1/2"  cotton  duck  pad  is  placed 
bstwesii  tha  bottom  of  the  bean  and  the  t^p 
of  the  corbel,  A  second  pad  in  placed 
between  tha  ond  of  the  beam  and  the  face  of 
the  column.  Tha  beam  is  ther  postensioned 
to  the  column  with  four  19  threaded  rebar 
as  shown  in  Figure  6.  The  <19  bars  are 
placed  in  3*  lor.g  eleeves  in  the  beam  and 
threaded  int'j  lenton  cov plots  located  In 
tha  colunn.  These  bars  are  than  stressed 
to  33  ksi  which  will  coriprass  tha  pad  and 
insure  a  uniform  seal.  Th^s  detail  (Figure 
6)  exhibited  the  bast  performance  in  a  teat 
of  seven  different  praoest  connect ions 
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subjected  to  severe  dynamic  loading 
(Jayashankar,  1987) .  In  addition  a 
continuous  steel  angle  with  a  3/8"  diagonal 
strut  located  at  20"  on  center  is  welded  to 
the  face  of  the  corbel  and  the  bottom  of 
the  beam.  The  design  of  the  corbel  should 
follow  standard  practice  as  set  fourth  by 
the  PCX  Design  Handbook  (1985). 


The  second  detail  (PL-23,  Figure  7) 
incorporates  an  18"  by  20"  thick  beam  and 
column.  The  beam  is  placed  on  a  corbel 
which  protrudes  6"  from  the  face  of  the 
column.  A  1/2"  cotton  duck  pad  Is  placed 
between  the  bottom  of  the  beam  and  the  top 
of  the  corbel.  A  4"  gap  between  the  end  of 
the  beam  and  the  face  of  the  column  will  be 
grouted  solid.  The  beam  is  than 
postensioncd  to  the  column  with  four  19 
threaded  rebar  similar  to  detail  PL-io. 


Monolithic  Tee  Joint 

Two  types  of  tea  joint  details  ware  tested. 
The  first  detail  (BT-8A,  Figure  8) 
incorporates  an  18"  wide  by  20"  thick 
column  and  beam.  The  beam  and  column  are 
reinforced  with  0.6%  steel  (3  -  #7  bars,  fy 
■  60  ksi)  on  each  face.  Diagonal  steel  (2 
•  #6  bars)  is  placed  2"  from  the  exterior 
face  of  a  4.5"  diagonal  concrete  strut 
located  on  the  inside  corner  of  the  beam- 
column  interface.  The  17  flexural  bars 
extending  from  the  beam  are  anchored  in  the 
joint  with  standard  90  degree  hooks.  The 
column  reinforcement  is  continuous  through 
the  joint. 
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Figure  8  Connection  BT-6A 


so 


The  second  detail  (BT-8A,  Figure  9) 
incorporates  the  same  beam  and  column 
dimensions  with  identical  steel 
reinforcement  as  for  BT-6A.  However,  the 
diagonal  strut  is  increased  to  9"  and  the 
diagonal  steel  located  2**  from  the  face  of 
the  strut  is  increased  to  3-#7  bars. 
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Figure  9  Connection  BT-8A 


RESULTS  AMD  DISCUSSION 

A  summary  of  results  for  the  five 
monolithic  knee  joint  connection  details 
are  shown  in  Table  1.  Peak  deflections  and 
stresses  are  shown  for  concrete  and  steel 
elements  at  critical  locations  in  tho 
detail.  These  values  are  compared  to  the 
allowable  and  a  determination  is  made  on 
the  reliability  of  the  detail. 


Table  1  Monolithic  Knee  Joints 
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The  detail  with  a  12*  beam  and  column 
(BSL2A}  shows  values'  of  deflections  and 


stresses  which  exceed  the  allowable  in  all 
respects.  The  size  of  the  members  in  this 
detail  are  not  sufficient  to  resist  the 
design  load.  Detail  BL-IA  required  a  large 
amount  of  steel  to  obtain  acceptable 
deflections  and  large  bar  stresses  were 
recorded  at  the  center  of  the  connection 
(element  #64)  .  For  these  reasons  a 
diagonal  strut  such  as  shown  in  details 
(BL-6A) ,  (BL-6C) ,  and  (BLS6A)  should  be 
added.  Detail  BL-6A  satisfied  all  the 
above  criteria  with  a  relatively  low  amount 
of  steel  reinforcement  (0.6%).  Increasing 
the  amount  of  diagonal  steel  from  half  of 
the  main  flexural  steel  has  little  effect 
on  the  connections  performance  (see  detail 
BL-6C) .  The  connection  which  performed  the 
best  out  of  the  five,  incorporated  "Sifcon" 
material  in  the  joint  region.  Both 
deflections  and  bar  stresses  in  the  joint 
region  decreased  dramatically. 

A  summary  of  the  precast  knee  joint 
connection  behavior  is  shown  in  Table  2. 
Detail  PL-ID  transfers  all  compressive 
loads  through  a  flexible,  reinforced  fiber 
pad.  Deflections  are  maintained  at  an 
acceptable  value,  however  compressive 
stress  in  the  top  reinforcement  is  well 
into  yield  which  may  cause  buckling.  The 
stress  in  the  tensile  reinforcement  does 
not  reach  yield,  however  the  stress  is  20% 
higher  than  allowable.  Reducing  the  amount 
of  prestress  to  10  ksi  should  alleviate 
this  problem. 


Table  2  Precast  Knee  Joints 
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Details  PL>2S  and  PL-4S  are  grouted  in 
place  with  a  6000  pal  concrete. 
Compressive  stress  in  the  top  reinforcement 
is  wall  within  allowable.  Tensile  stress 
in  the  bottom  bars  is  vary  near  yield. 
Reducing  the  amount  of  prastress  in  (PL*‘4S) 
to  10  ksi  has  lowered  the  stress  to  near 
acceptable  values.  Stiffening  the 
compressive  region  of  this  connection  with 
grout  has  greatly  increased  the  tension  on 
the  continuous  angle  assembly  located  at 
the  corbel.  The  spacing  of  the  stiffeners 
on  this  assembly  should  be  decreased  to  10* 
on  center. 
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The  results  for  the  monolithic  tee  joint 
are  shown  in  Table  3.  Detail  BT-6A  uses  a 
4  1/2"  diagonal  strut.  The  tensile 
reinforcement  in  the  beam  is  stressed  well 
into  the  yield  range.  Therefore  we 
increased  the  strut  size  to  9"  and  the  area 
of  the  diagonal  steel  to  (0.6%)  in  detail 
BT-rSA.  This  reduced  the  tensile  stresses 
in  the  joint  region  and  relocated  the  hinge 
to  an  acceptable  distance  from  the  joint 
interface.  Both  details  have  shown  high 
flexural  shear  stresses  at  the  interface 
and  interior  of  the  joint.  This  is  due  to 
a  stiff  column  flexible  beam  configuration. 
Correct  design  of  shear  stirrups  at  these 
regions  of  concentrated  moment  and  shear  is 
a  necessity.  The  reinforcement  in  the 
column  was  subject  to  very  light  stress  (20 
ksi  max. ) . 


Table  3  Monolithic  Tee  Joints 
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ABSTRACT:  In  this  papar,  a  rata-lndapondant 
constltutlva  aodal  for  plain  oonerata  la  propoaad 
for  application  to  the  analyala  of  lapulaa  loaded 
structural  aaobars.  The  nodal  eonblnss  a  conclnuun 
danaga  approach,  using  a  scalar  daaaga  varlabla, 
with  a  presaura  aanaitlva  plasticity  aodal.  Tha 
plasticity  nodal  Inoorporataa  a  nonasaoclatad  flov 
ruls  In  raglona  of  low  oonpraaslva  or  tanslla 
hydrostatic  prasauxas  and  an  aaaoclatad  flow  rula. 
Tha  eoncrata  nodal  la  ooablnad  with  a  uniaxial 
steal  nodal  and  a  layarad,  large  strain,  Tlnoshanko 
bean  alanant  to  parfom  tha  analysis  of  Inpulse 
loaded,  sinply  supported,  ralnforoad  oonorata 
baana. 


nmoDocnoii  amo  lACEORoras 

The  analysis  of  tha  rasponaa  of  Inpulsa  and 
blaic  loaded  itruoturas,  burltd  and  abova  ground, 
has  cocaivad  a  continuous  but  varying  lavol  of 
attention  over  the  past  40  yaara.  Due  to  tha 
alnost  Infinite  nunbar  of  pamutationa  of  a  given 
structure's  paranatara  (l.a,  gaonatry,  natsrlal 
propartlaa,  depth  of  burial,  ate.)  and  antlolpatad 
thraata,  and,  due  to  tha  costa  of  parfomlng  full* 
and  snall*aoala  testa  on  such  struoturaa,  tha 
aswunt  of  available  oxparlnantal  dau,  while  broad, 
la  alto  scant  rolatlva  to  any  particular 
oonhlnatlon  of  attuoturo  and  Inpulaa  load.  Thus, 
thn  davalopnant  of  analytloal/oonputatlonal  tools 
are  required  In  order  to:  1)  understand  tha 
oonplax  nonlinear  behavior  of  tha  structure  (and 
toll.  If  burled),  2)  parfom  paraiMtclo  studies, 
and  3)  develop  design  guldallnaa. 

One  analytical  tool  la  the  Finite  Blamnc 
Method,  which  has  bean  applied  with  great  tucoasa 
to  gaoaatrlcally  and  utarlally  nonlinear  oontlnuut 
and  ttruotural  problana,  Tha  gaoMtrio 
nenltnearttlaa  can  bo  uodellad  with  tha  wall  known 
klnSMtle  fomulatlont  of  large  dlsplactMnt 
analyaat  (Lagranglan,  Cularlan,  Updatad 
Ugranglan);  the  Mtarlal  nonllnearltlas  ara  not  as 
taslly  accountod  for.  'Uhlla  adequate  *00010  exist 
for  tha  nonlinear  bahavlor  of  steal,  the 
davalopMnt  of  an  acourata,  self*contttcanc,  and 
unified  aodal  for  tha  nonlinear  rasponaa  of  plain 
eoncrata  la  still  an  ectlva  area  of  intaraat. 

In  order  to  achUva  a  high*qualley  nuaerleal 
solution  CO  a  given  boundary  value  problaa,  an 
iegirovad.  rate  independent,  strain  softening. 


Theodor  Krauthsimar 
Assoc.  Prof,  of  Civil  Engrg. 
University  of  Mlnnasota 
Minneapolis,  MN  55453 


nonlocal  Continuum  Danaga/Flastlcity  (CDF)  model 
for  plain  concrete  is  presented  in  this  paper  and 
certain  attributes  of  this  model  are  developed  to 
facilitate  its  application  in  structural  elamenta, 
such  as  beans,  plates,  and  shells. 

Given  the  complex  bahavlor  of  concrete,  it  is 
not  ourprlslng  that  a  large  number  of  distinctly 
different  constitutive  models  for  concrete  have 
been  proposed  over  tha  years.  Hare,  only  the 
Continuum  Damage  Meohanlos/Flaatlalty  models  and 
models  for  strain  softening  are  discussed. 

The  nonlinear  response  of  concrete  is  created 
through  tha  combination  of  mlerooraok  growth  and 
frictional  slip.  Tha  theory  of  Continuum  Damage 
Mechanics  accounts  for  tha  phenomena  of  strength 
and  stiffness  degradation  due  to  mlorooraoklng;  the 
plastic  flow  and  prapaak  nonlinearity  of  conorato 
can  be  modelled  with  the  theory  of  Plaatlolty. 
Recently,  the  two  theories  have  boon  combined  to 
form  a  nunbar  of  successful  approaches  for 
modelling  plain  concrete.  One  of  the  first 
acb.inpts  at  combining  the  damage  and  plasticity 
phenonona  was  made  by  Baxant  and  Kin  '79,  who 
merged  a  oonvantlonul  plaatlolty  approach  with  a 
strain  basad  fracturing  theory.  Utar,  Han  and 
Chen  *67  combined  a  stress  bated  hardening 
plasticity  theory  for  the  pca*fallure  response  of 
concrete  with  a  strain  based  plastlo/fraoturlng 
thaory  for  the  response  after  failure.  Using  tha 
Internal  variable  theory  of  themodynanlca ,  Yacdanl 
and  Sehreyer  '87  developed  a  Von  Hlses  plasticity 
surface  and  a  mean  presaura  saniltlve  damage 
surface  that  aoconmodetes  the  two  nodes  (shear  and 
tensile)  of  cracking  by  caking  Into  account  tha 
current  stress  stats  at  the  crack  surface  as 
detamlned  by  the  orthogonal  projections  of  the 
stress  censer.  Also,  Slmo  at  al.  '87  coupled  a 
conventional  stress  based  Cap  modal  with  a  scalar 
damage  variable  that  evolves  using  tha  rule 
suggested  by  Nazars  '82;  this  approach  yielded 
excellent  results,  lastly,  Franttlskonls  and  Oesal 
'87  combined  a  scalar  damage  approach  with 
plesclclty  theory  and  found  good  egreeaenc  between 
predictions  and  test  results. 

The  application  of  softening,  race  independent 
constitutive  medals  In  the  creataanc  of  initial 
boundary  value  problams  leads  to  e  loss  of 
hypsrboliclty  In  the  governing  equations  of  motion 
and  to  mesh  dependency,  local Ltet ion  of  strains, 
end  erroneous  predictions  of  energy  dissipation  In 
numerical  calculations;  thus,  the  potential  for 
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such  negative  side  effects  exists  for  each  of  the 
models  mentioned  in  the  previous  paragraph.  As 
discussed  by  ^ead  and  Hegemier  ‘84,  strain 
softening  in  concrete  is  the  direct  result  of  the 
formation  of  discrete  internal  structure  due  to 
microcracking.  Thus,  in  the  softening  range, 
concrete  can  no  longer  be  represented  as  a 
continuum  and  some  form  of  representation  of  the 
internal  structure  must  be  Included  in  the 
constitutive  model.  In  order  to  overcome  this 
difflcuty  in  modelling,  various  methods  for 
enforcing  finite  energy  dissipation  over  a  discrete 
region  of  the  body  and  for  removing  the  mesh 
dependency  have  been  proposed;  these  include:  a 
composite  damage  formulation  which  incorporates  a 
"daMge  volume”  and  distinguishes  between  tensile 
and  shear  cracking  (Wiliam,  Bicanic,  and  Sture 
'  84) ;  the  introduction  of  higher  order  spatial 
derivatives  into  the  strain  displacement  relations 
(Lasry  and  Belytschko  '87);  and,  the  introduction 
of  the  strain  gradient  into  the  definition  of  the 
strength  or  yield  function  (Schreyer  and  Chen  *86). 

One  method  of  particular  interest  involves  the 
use  of  "nonlocal”  constitutive  laws,  in  which  the 
dependant  variable  at  a  point  is  not  solely  a 
function  of  the  state  variables  at  that  point  but 
rather  depends  on  what  occurs  in  the  neighborhood 
of  chat  point.  Such  an  approach  was  formulated  for 
concrete  by  Bezant  ec  al.  ‘84,  who  averaged  both 
strasa  and  strain  over  the  neighborhood  of  the 
point  and  used  the  nonlocal  definitions  of  the 
stress  and  strain  directly  in  the  governing 
equations.  Recently,  Pijaudler* Cabot  and  Bezant 
'87  applied  an  averaging  approach  to  a  nonlocal 
definition  of  the  daatage  variable  in  a  Continuum 
Damage  Hechenies  model.  With  this  approach,  Che 
field  equations  have  the  standard  form  and  no  extra 
boundary  conditions  are  needed.  Their  work  ia 
particularly  interaating  in  chat  they  suoeaaafuUy 
combined  their  nonloccl  spproach  with  a  layered 
beam  element  to  analyse  the  response  of  a 
statically  loaded  beam;  they  showed  chat,  due  to 
Euler's  hypothesis  of  plane  cross  sections,  the 
length  for  damage  averaging  can  not  be  less  chan 
the  beam  depth.  It  should  be  noted  that  Che  use  of 
the  beam  depth  for  averaging  the  damage  can  be 
validated  physically  by  exatilning  the  hinge 
fomaCiona  in  beams  that  have  bean  loaded  to 
failure  (Corley  '66);  in  most  oases,  the  hinge 
length  is  on  the  order  of  8S  to  93  peccant  of  the 
total  depth. 

As  is  well  known,  rotatory  inertia  and  shear 
affects  significantly  alter  the  response  of 
structural  elements  that  are  loaded  over  a  vary 
short  time  period  relative  to  the  element's  natural 
period,  and,  chetefora,  a  layered,  large 
displacement,  large  strain,  Timoshenko  beam  element 
is  developed.  The  proposed  oonstitutive  model  ia 
combined  with  this  Timoshenko  beam  element  end  a 
■train  hardening  plasticity  model  for  the 
reinforcing  steal  to  succaeafully  analyze  two ' 
simply  supported,  impulse  loaded,  reinforced 
concrete  beams  taiiced  by  Feldman  and  Sless  'SB. 

A  IIOHLOCAL  COKTKnnm  DAMACE/rUStXGXTT  HODEI.  FOR 
rUIM  CbilCUTg 

The  concrete  model  developed  herein  is  a 
combination  of  Continuum  Damage  Hechanies  and 
Rleaticity  theories  and  it  uses  a  nonlocal 
definition  of  a  scalar  damage  variable  to  modal 


strain  softening.  Various  researchers  have 
combined  plasticity  and  continuum  damage  approaches 
to  successfully  reproduce  test  results  of  material 
elements;  however,  without  a  consistent  technique 
for  stabilizing  the  strain  softening  calculations, 
they  appear  doomed  to  the  sane  shortcoming 
discussed  above  (mesh  dependency,  unrealistic 
energy  dissipation,  etc.)  when  applied  to  realistic 
boundary  value  problems .  The  nonlocal  approach 
described  in  the  following  successfully 
incorporates  strain  softening  effects. 

The  Continuum  Damage  Mechanics  portion  of  the 
approach  developed  herein-  and  discussed  by  Stevens 
and  Krauthammer  '88  (to  appear) ,  is  based  loosely 
on  the  approach  proposed  by  Frantzlskonis  and  Desai 
'87,  in  which  tha  strength  and  stiffness 
degradation  are  modeled  with  a  scalar  damage 
parameter.  It  can  rightfully  be  argued  that  a 
complete  three  dimensional  treatment  of 
ffllcrocracklng  effects  on  the  response  of  concrete 
requires  a  higher  order  definition  of  the  damage 
variable,  and,  first,  second,  fourth,  and  even 
eighth  order  tensors  have  been  proposed.  However, 
in  cases  where  damage  directionality  is  not  a 
dominant  feature,  such  as  in  the  plane  strain  or 
plane  stress  conditions  of  beams,  frames,  end 
slabs,  a  scalar  parameter  is  sufficient  (Resande 
'87). 

The  determination  of  how  the  averaging  of 
damage  is  to  be  performed  (in  order  to  develope  a 
nonlocal  damage  parameter)  and  tha  value  of  the 
characteristic  length,  i,  (over  which  the  averaging 
is  to  be  performed)  are  based  on  considerations  of 
cha  proposed  application  of  tha  oonatitutlva  modal. 
As  PiJaudlsfCabot  and  Bezant  *87  show,  the  uss  of 
the  Eular  hypothesis  in  a  beam  analysla  raquiras 
that  tha  averaging  length  be  greater  or  equal  to 
Che  beam  depth.  Tharafort,  alnca  this  constitutive 
model  will  be  Implemented  Into  a  Timoshenko  beam 
tlsmsnc,  which  Is  based  on  tha  Eular  prlnolpla  and. 
In  which,  the  shear  strain  la  rtpctssntsd 
klnsmstloslly  by  a  rotation  of  tha  orots  section 
similar  to  the  flsxural  rotation  (as  dlacusssd 
latsr) .  tha  damags  averaging  Is  performed 
unlaxlally  at  each  Isysr  of  cha  beam  alamsnt,  In 
Che  direction  perellel  to  the  long  exle  of  the 
beem,  end  the  cheteoterletlo  length,  i  la  taken  aa 
the  beem  depth.  Thla  one  dlmenaional  averaging 
scheme  rectclcta  the  pooslbllley  of  itreln 
softening  through  the  depth  of  tha  beam;  thue, 
while  the  nonlocal  lealar  deaege  parameter  doea 
degrade  both  ahear  and  normal  atreBees.  a  monoconlc 
loading  path  of  pure  aheer  at  a  orose>e«ocion  of 
the  beam  will  result  In  a  constitutive  reeponse 
Chet  Is  elsstle/etraln  hardening  plaeclc  In  sheer, 
providing  damaga  hae  not  occurred  at  Cha 
Integration  polnca  adjacent  to  the  crose  saoclon. 
Thle  shortcoming  ml^t  be  oveceoma  through  use  of 
two  nonlocal  acaler  damage  parsmecsca,  rsprsssnclng 
csnslls  end  shear  deimaga,  teparacaly.  However,  ee 
the  reeulta  presented  later  ehcw,  ehla  uniaxial 
damage  averaging  ceohnlqus  la  aul cable  for  the 
bsams  that  ware  analyzed. 

The  next  leaue  la  the  fora  of  the  evolution 
aquatlona  for  cha  local  damage  parameter:  In  this 
approach,  the  evolution  equations  will  be 
paramecerlzod  using  the  concept  of  equivalent 
tensile  Strain  (Mazars  '82;  Ortiz  '83;  Slmo  and  Ju 
•87). 

Tho  reaponsa  of  cha  topical  (undamaged) 


nuit«rt«l  Is  controlled  by  the  plasticity  portion  of. 
the  CDF  model,  which  employs  a  continuous,  smooth 
yield  surface  that  combines  a  strain  hardening 
modified  DtMcker*Prager  failure  surface  with  a 
curved  cap;  this  model  was  developed  by  Schreyer 
and  Bean  '87,  who  refer  to  It  as  a  "Frager-Orucker* 
cap  model  to  differentiate  It  from  the  Drucker- 
Frager  cap  model  commonly  used  for  modelling  soils. 

A  typical  modified  Frager-Drucker  yield  surface 
In  the  J{ ,  F  plane  Is  plotted  In  Figure  1 . 

Over  the  majority  of  the  yield  surface,  the 
plastic  strain  tensor  evolves  through  use  of  an 
associated  flow  rule;  however,  when  the  model  Is 
specialized  from  the  three  dimensional  state  to  a 
plane  strain  state  (as  Is  typically  assumed  for 
analyses  of  structures  that  are  "long"  In  the 
direction  perpendicular  to  the  plane  of  loading) ,  a 
nonastoclated  flow  rule  provides  more  reasonable 
results  In  regions  of  tensile  hydrostatic  pressures 
and  low  compressive  hydrostatic  pressures.  Under 
plane  strain  conditions,  the  associated  flow  rule 
generates  dilation,  which  creates.  In  turn,  an 
unrealistic  confinement  stress  In  the  out>of-plane 
direction;  since  concrete  Is  a  pressure  sensitive 
material,  this  out*of*plane  compressive  stress  can 
lead  to  erroneous  predictions  of  shear  and  tensile 
strengths . 

In  order  to  obviate  this  shortcoming,  the 
concrete  model  proposed  herein  uses  a  nonassoolated 
flow  rule  In  regions  of  tensile  or  low  compressive 
(lets  than  0.1  fj)  hydrostatic  pressure  to  reduce 
the  dllatlonal  strain  predictions  while  at  the  same 
time  allowing  adequate  pradletlons  of  strength. 

The  nonassoclated  flow  rule  Is  Implemented  with 
a  plastic  potential  that  Is  chosen  as  another 
modified  Druoker-Pragsr  surface  as  shown  In  Figure 
2. 

Thera  are  two  shortcomings  to  the  use  of  a 
nonassooiatod  flow  rule.  First,  nonstsoclatsd  flow 
rules  lead  to  asymmetric  tangential  elastlfplastlc 
compllanos  tensors,  which.  In  cum,  are  difficult 
to  accommodate  oomputatlonally  In  the  typical 
solution  schemas  used  In  Finite  Element 
applications.  However,  since  the  cantral 
difference  tsehniqua  la  used  to  Integrate  the 
squatloni  of  motion  In  the  Flnlta  Blamant  approach 
discussed  later,  the  tangential  oompllanoa  la  never 
explicitly  caloulatsd  or  Inverted  and  this 
difficulty  Is  avoided  on  the  computational  level. 
Second,  and  most  Importantly,  Che  main  criticism 
lavalad  against  the  use  of  nonassoclated  flow  ruloa 
Is  the  potential  for  spontaneous  energy  generation. 
This  unpleasant  fsaturn  was  first  discussed  by 
Xl'lushln  *61  (a  summary  of  11‘lushln's  proof  la 
prasencad  by  Sandler  and  Rubin  *87). 

The  modified  Fragar*Dcucker  plasticity  modal 
has  the  following  advaneagas;  1)  strain  hardening 
Is  predicted  for  all  paths,  2)  more  ductility  la 
pradlcced  for  paths  assoolatsd  with  largo  maan 
prosauros,  3)  the  yield  surfaoa  Intersects  the 
hydroscatlo  axis  at  right  angles  for  both  positive 
aitd  nsgatlva  maan  pressuras  so  no  special  algorithm 
Is  nasdad  at  chase  Intersection  points,  and  t)  the 
flow  surface  la  continuous  and  has  a  continuous 
derivative  everywhere  so  a  comer  algorithm  is  not 
requltad  (Schreyer  and  Bean  *87).  Although  the 
number  of  required  material  constants  la  large,  the 
values  may  be  deduced  directly  from  conventional 
uniaxial  and  criaxlal  test  results;  also,  Schreyer 
and  Bean  *87  and  Stevens  and  Krauthammer  ‘88 


present  representative  values  for  low,  medium,  and 
high  strength  concretes. 

The  proposed  concrete  model  (using  the  local 
definition  of  the  damage  variable)  was  evaluated 
successfully,  on  a  material  level,  through 
comparisons  with  three  sets  of  concrete  test  data; 
for  a  full  presentation  of  the  results,  the  reader 
Is  directed  to  Stevens  and  Krauthammer  '88. 

TIMOSHENKO  BEAM  ELEMENT 

The  kinematic  relations  for  the  layered 
Timoshenko  beam  element  are  based  on  the  usual 
assumptions;  1)  plane  sections  originally  normal 
to  the  n'.<utral  axis  remain  plane  after  deformation 
but  not  necessarily  normal  to  the  neutral  axis  and 
2)  the  shear  strain  Is  constant  over  the  beam's 
cross  section  and  may  be  defined  as  a  rotation  of 
the  cross  section.  Assuming  small  to  moderate 
rotations  and  small  axial  strains  and,  using  the 
definition  of  the  Green- Lagrange  strain  tensor,  the 
nonzero,  nonlinear,  strain  displacement  relations 
may  be  written  as 
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where  U  m  U(f^)  la  the  axial  displacement  of  the 
neutral  axis  along  the  beam,  7  •  7(fi)  la  the  shear 
dafonaatlon,  fi  la  tha  total  rotation  of  cha  croaa 
sactlon,  V  •  la  tha  transvezae  displacement 

of  the  neutral  axis  along  the  beam,  are  the 
local  beam  coordinates  (1  ••1,2),  i  •  V‘,  and  (  )*  ■ 
derivative  w.r.t.  the  Independent  variable.  In 
this  formulation,  the  axial  dlaplaoament  and  tha 
shaar  rotation  are  Interpolated  with  linear  shape 
functions  acting  on  the  nodal  varlablaa,  and  the 
transverse  dlsplscement  Is  Intarpolatsd  using  Chs 
standard  cubic  Htrmetlen  polynomials  acting  on  the 
total  rotations,  or  fit  ~  whare  and 

art  tha  flaxural  and  shear  rotations,  raspeotlvely, 
sc  nods  1  of  the  beam.  The  sheer  rotations  are 
retained  us  Independent  degrees  of  frsadom  end 
continuity  of  ehear  scrota  alementa  la  pretarvad. 

In  the  Implementation  of  the  Flnlta  Element 
method  using  central  difference  rfme  Integration, 
tha  oelouleclon  of  the  Intamel  force  vector  must 
be  performed.  The  Internal  .^oroe  Victor  is 
equivalent.  In  a  virtual  work  cans#,  to  tha 
Internal  alement  etreiaea  created  during  tha 
deformation  of  the  element  and  It  Is  daflnad  as 


OlM®!  <lVol  . (3) 
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whars  IFih^I  is  the  Internal  force  vector 
equlvalsnc  to  the  alsmant's  stress  state,  [B|*  la 
the  three  dlmenalonal  strain  displacement  matrix, 
and  [ff|  la  the  Cauchy  stress  tensor  (matrix).  For 
Incagratlon  through  the  depth,  the  element  Is 
subdivided  Into  10  layers  and  the  Integration  la 
performed  by  s  summation  of  the  stresses  at  each 
layer;  the  streeees  In  the  eteel  end  concrete  are 
calculated  with  tha  modale  discussed  above  and  e 
perfset  bond  la  assumed  between  the  reinforcement 
end  conoreca.  Tha  Incegretloa  over  the  length  of 
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tha  «l«iMnt  Is  parfornad  using  thraa  point  Gaussian 
Intagratlon.  Also,  tha  movaaant  of  tha  nautral 
axis  dua  to  tha  natarlal  nonllnearltles  Is  trackad 
with  an  itaratlve  technlqua  davalopad  by  Pugllsl 
and  Krauthumor  '87. 


BEAM  AMALTSIS  Affi)  BISGUSSION 

The  Tlaoshanko  beaa  alament  was  combined  with 
tha  proposed  CDF  concrete  modal  and  a  uniaxial 
alastlc/straln  hardening  plastic  model  for  the 
reinforcement  steel  (sea  Park  and  Paulay  '75);  the 
shear  contribution  of  tha  steal  rainforcemant  is 
assumed  small  compared  to  the  shear  force  resisted 
by  the  concrete  and,  thus,  tha  steel  response  Is 
modelled  as  uniaxial.  The  resulting  approach  was 
applied  to  the  analysis  of  simply  supported 
reinforced  conoreca  beams,  which  were  subjected  to 
large  aaplltuda,  short  duration  loads  at  midspan 
(Feldman  and  Sloss  '58).  Figure  3  presents  the 
dimensions,  cross  sectional  geometry  and  some  of 
the  Inscrumontatlon  for  beams  HI  and  Cl;  Figures  4 
and  5  present  the  linear  approximations  of  the 
applied  loadings  of  beams  HI  and  Cl,  respectively. 
During  casting,  both  beams  experienced  considerable 
damage  and  large  permanent  deflection.  The 
residual  mldspan  displacements  were  approximately  8 
Inches  and  2.5  inches  In  beams  HI  and  Cl, 
respectively. 

Due  to  the  symmetry  of  the  beasts,  loads,  and 
boundary  conditions,  only  one  half  of  the  beam  was 
modeled,  using  the  coarse  and  refined  meshes  shown 
in  Figure  6.  At  the  right  hand  node,  the 
horlsoneal  displacement,  shear  rotation,  end 
bending  rotation  ere  rascreinad. 

As  is  typical  with  sioet  test  reports,  only  the 
oompressive  strength  of  tha  ooncrete  and  the  yield 
atrain  end  stress  of  Che  steel  were  reported  by 
Feldman  and  Sitae  '58.  Tha  ooncrete  compreesive 
strengths  of  beasts  HI  and  Cl  vers  3960.0  pal  and 
5830.0  psi,  reepectively,  so  “typical''  values  for 
stadium  strength  concrete  are  used  in  Che  analysis, 
as  listed  in  Table  1.  Also,  typical  values  for  the 
reinforcing  steel  were  assumed,  as  listed  in  Table 
2.  the  dsneitles  of  the  concrete  and  steel  were 
taken  as  ISO  pof  end  490  pof,  raapactively. 

The  equations  of  motion  were  integrated 
temporally  with  the  conditionally  stable  cantrel 
difference  technique;  Che  integraciofl  time  seep  was 
chosen  as  half  the  orteioal  time  step  (which  was 
governed  by  the  flexural  frequency). 

Comparisons  bacweeo  the  experimental  and 
oaloulated  displacement  histories  of  Cages  3,  4, 
and  5  of  beam  HI  are  shown  in  Figure  7  for  coarse 
mash.  Meeh-eensltivity  was  checked  by  using  the 
model  with  refined  mesh  end  it  was  obrorved  that 
this  algorithm  is  not  sensitive  to  the  fineness  of 
the  model.  As  this  figure  shove,  tnere  is  good 
agreement  b.ttween  peek  displacement,  residual 
displacement  and  the  time  when  peak  displscsaencs 
occur. 

Thi  rssults  in  Figurs  7  when  considered 
elongwith  the  meeh- independence  of  thle  elgorithm 
end  the  feet  thet  e  significant  smount  of  strain 
softsnlng  has  occurrsd  et  midspen,  show  that  ths 
campuescionsl  difflcultlea  normally  peessne  with 
strain  softening  models  have  been  ellevleted  with 
the  nonlocel  CDF  model.  Xn  order  to  verify  the 
benafite  of  the  nonlocal  approach,  beam  HI  waa 


reanalyzed  using  the  coarse  and  refined  meshes  and 
Che  CDF  model  without  the  damage  averaging;  the 
computed  histories  of  the  midspan  displacement  for 
the  local  and  nonlocal,  coarse  and  refined  mesh 
analyses  are  plotted  In  Figure  8.  As  shown  in  this 
figure,  the  local  strain  softening  model 
overestimates  the  displacement  with  the  coarse 
mesh,  and,  with  the  refined  mash,  the  displacement 
appears  to  continue  to  Infinity.  An  examination  of 
the  output  for  these  local  analyses  showed  that  tha 
damage  localized  Into  the  element  directly  next  to 
the  load;  for  the  refined  mesh,  this  element  became 
a  kinematic  hinge  and  hence  the'  Increasing 
displacement. 

Comparisons  between  the  measured  reaction  at 
the  left  end  of  the  beam  and  the  calculated 
reaction  using  the  coarse  mesh  for  beam  HI  are 
shown  in  Figure  9.  Again,  there  is  very  good 
agreement  with  the  peak  reaction  and  the  time  of 
peak  reaction,  as  well  as  the  time  when  the 
reaction  force  returns  to  zero.  However,  after 
approxlmecely  60  msec,  the  calculated  reaction 
osclllstss  with  s  higher  frequency  chan  the 
measured  reaction;  this  difference  is  due  Co  Che 
fact  that  Che  unloading  stiffness  predicted  with 
Che  continuum  damage/plasciclcy  model  Is  coo  large 
(as  discussed  in  Stevens  and  Krauthammer  '88), 
resulting  in  a  stlffer  beam  with  s  higher  nscursl 
frequency.  A  correction  to  Che  constitutive  model 
would  Improve  the  results  for  the  response  after  60 
msec;  however,  Che  key  point  is  that  the  peak  vslua 
is  corrsctly  prsdictsd  as  i.s  ths  time  et  which  the 
peak  occurs  and  ths  elms  wbsn  ths  rssotion  cscucns 
to  zero. 


(MMGIUSIOHS 

A  nonlocal  Continuum  Dsasge/Flaseioicy  modal 
has  bsan  dsvslopsd  and  implamntsd  for  ths  snslyasa 
of  rsinforcsd  oonersts  bssms.  Xn  ths  dsvslopmenc 
of  tha  modal,  prsviously  proposed  focmulstlons  wees 
modifisd  and  combined  into  s  comprehensive 
constitutive  relationship.  This  modal  has  bean 
vsriflsd  for  various  strsss  paths  and  by  ths 
snslys.ls  of  rtinforosd  oonersts  bssms  under  impulss 
loading. 

Ths  rssults  of  ths  snslysss  of  tns  two  impulse 
loaded  beams  agrss  quits  wall  with  the  experimental 
data.  Tha  snslysss  prtdicts  ths  dlsplscsmsnt, 
reaction,  end  strain  historlss  very  olossly.  The 
observed  mesh  indspendsnes  la  particularly  plsssing 
considering  the  shortcomings  usually  associated 
with  strain  loftsning  constituclvs  models. 

Xn  oonelusion,  it  appears  that  the  oomblnselon 
of  cha  Timoshenko  beam  element  with  the  nonlocal 
COP  model  and  ths  unlnxlsl  steal  modal  is  capable 
of  rapraaantlng  the  dynamic  bshsvlsr  of  reinforced 
ooncrete  beams  quite  scourstsly. 
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Abthncb 

This  pq>er  sonunttizes  the  result  of  a  Phase  I  SBIR  re¬ 
search  efiigtt  widi  an  overall  objective  of  reviewing  and  analyzing 
stochastic  process  methodologies  for  their  appUcability  to  protec¬ 
tive  stnicttire  design  [1],  The  reseaich  focused  on  the  following 
tasks:  (1)  develop  an  overview  of  the  phases  of  protective  struc¬ 
ture  design  and  identify  speciGc  areu  where  probability-based 
methods  are  applicable;  (2)  identify  the  steps  required  for  the  de¬ 
velopment  of  reliability-based  dt^gn  methods  for  protective 
structures;  and  (3)  illustrate  the  application  of  reliability  methods 
to  the  anatysis  of  an  aboveground  structure  and  a  buried  sttuctum 
The  IluM  I  tesulta  indicate  that  Aiodamental  improvements  to  the 
daetp  of  protective  structum  can  be  achieved  through  research 
amiqipUcatiaa  of  modem  coocepu  of  structural  teliablUfy^  These 
improvemeius  include;  (1)  systematic  jdantiflcatinw  of  uncer- 
laMae  and  failure  modm;  (2)  development  of  leUabllity-based 
dc:dgn  factors;  ind  (3)  development  of  in  analysis  cool  Cor  suivlv- 
ible  structures  tad  faeilitifis. 


USES  OFPROBABlUSnC  METHODS  IN  PROTECTIVE 
DESIGN 

The  protective  dasiipt  process  involves  a  number  of  activl- 
liec  or  phiiM  that  tnteraa  in  vaitoui  ways  over  the  design  life 
cycle  of  •  facility.  Rve  basic  phases  wete  tdcatifled:  (Othreu 
amdustioiu  (2)  tnission  snd  perfotiAsnci  criteria  developm^ 
design  and  analysis;  (4)  cotuKrucUon;  and  <3)  teat  and  operation. 
Thaie  phases  cover  ^  entire  range  of  opeiatiotui  of  protective 
design,  PoUowtegisabtiefsummaty  CromRerereoce  1  ofeach 
phase  with  tespem  w  key  uncertsittiki  and  the  potential  applica¬ 
tions  of  probabiUitic  mmhods. 


Threat  Evaluation  Phase 

The  objective  of  threat  evaluation  is  to  dmermine  the 
wartime  attack  strategy  and  types  of  weapons  to  which  the  facil- 
ity/structute  would  likely  be  subjected.  Three  general  fypes  of 
uncertainties  are  associated  with  threat  evaluation:  (1)  attack 
straiegiea  (decision  vitiables  such  u  aimpoints,  number  of  sorties, 
and  we^xHt  types)  that  the  enemy  force  might  employ  on  a  facil¬ 
ity:  (2)  the  modeling  or  prediction  etiois  asaociated  with  a  model 
of  attack  and  enemy  weapon  characteristics;  and  (3)  die  random 
uncertainties  (such  at  C^,  fuae  performance,  bomb  autvivabti- 
Ity)  about  the  mean  values  of  precUctedpcrfonnance.  Probabilis¬ 
tic  methodologlea  with  direct  ajqiUcation  to  treating  and  analyzing 
these  uncrrtairitiea  include  (1)  decision  analytia  arid  game  ih^; 
(2)  Mome  Carlo  Simulation;  and  (3)  exoeme  value  theory. 


MisahMt/Paifomiuca  Critaria  Davalopnaat  Phaaa 

In  the  misiion  and  petfotminca  critaria  developmeiu 
phasa,thesiaicdmiitiooistran3laiedlniotBaio(deiigRi«quira> 
mous  for  each  facility,  which  can  be  used  to  develop  deatgA  crita¬ 
ria  for  individual  structure.  General  design  criteria  specify  the 
icceptable  structural  perforaunce  in  terms  of  Umiting  damage 
relative  to  internal  equipment  or  pataormel  In  addition  to  the 
general  design  criteria,  there  are  often  requirements  on  derign 
pnrfomuncc.  These  requirements  may  be  treated  as  constraints 
that  the  design  must  sailify.  0>at.sorv{vibiUfy,ntaititiittability. 
rellabiltty.  and  cooMtuctability  are  riements  ihM  may  be  expreiaed 
u  general  criteria  or  cornets  on  the  dmign  performance. 
These  pararacttn  are  essentially  decirion  variabiea  that  are  eitab- 
Ushed  during  the  criteria  development  phase.  These  vatiibles 
specify  the  requirements  of  the  facilify,  as  such,  there  is  gen¬ 
try  no  need  to  treat  uncertainties.  However,  un'»tta*flties  that 
result  Crom  the  predicted  peifotmince  of  s  derign  to  meet  titeia 
criterit  nuy  be  signUlcant  and  should  be  treated  in  the  analyiis 

rib«M 


Design  and  Analysis  Phase 

Hie  design  and  analysis  phase  is  the  central  activity  of 
protective  design  and  the  one  with  the  greatest  focus  of  researoh 
activity.  It  can  be  subdivided  into  three  elements:  design  synthe¬ 
sis,  system  failure  mode  analysis,  and  design  analysis.  Design 
synthesis  begins  with  the  development  of  a  protective  design 
concqit  based  tm  the  ^wcified  threats  and  the  general  design 
requitements.  Stnictural  form,  material,  construction  method, 
lotion,  and  extenial  protection  are  attributes  that  the  designer 
can  select  in  the  develt^ent  of  the  design  conc^L  Rerearch  in 
this  area  deals  with  basic  topics,  such  as  improved  materials,  as 
well  as  innovative  concepts  for  protective  constructitm  and  basing 
altemadvos.  In  addition,  research  tools  from  optimal  design  the¬ 
ory  and  artificial  intelligence  are  relevant 

System  failure  mode  analysis  involves  (he  identification  of 
potential  failure  modes  of  each  dement  and  devdoping  the  struc¬ 
ture  and  facility  failure  paths.  For  quantitative  analysis  and  com¬ 
plex  systems,  event  tree  and  fault  tree  methodohi^es  are  useful 
tools  to  perform  system  modeling,  identify  vuinerabilities.  and 
assess  consequences  of  damage  and  loss  of  function.  The  integra¬ 
tion  of  event  and  fault  trees  provides  an  analytic  approach  for 
systematic  identification  and  modeling  of  the  sequence  of  failure 
and  the  detomituiion  of  critical  dementt  of  a  complex  sys¬ 
tem.  They  provide  the  bes:  available  means  for  identifying  and 
understanding  ftciii^  dedgn  and  operation  far  a  manner  that  leads 
to  a  quantification  of  system  nliability. 

The  design  inatyais  phase  follows  the  development  of  the 
prdiminaiy  design  concept  and  the  identification  of  all  relevant 
biluremodes.  The  four  basic  ;Aeps  are:  (I)  obudn  htformation  on 
weapon  characteristics;  (2)  asresa  weapon  eifecu;  (3)  (le^op 
stnicural  loads  from  speeiSed  weapon  direeu.  chanttusistics  aa^ 
effects:  and  (4)  tnaly»  stniciural  response.  Thii  phare  is  the 
central  aotivify  of  many  itructurd  anafysisfdeiignos  for  cooven- 
tionai  weapon  effects. 

Reliabiliity-based  design  (RBD)  methods  provide  a  ra¬ 
tional  method  to  treat  uncertaimiea  in  tho  design  and  tnaiysls 
phase.  Cnee  these  methods  are  devetop&l  and  documented  in  a 
research  phase,  the  designer  would  work  with  deslgr.  fwttors  (such 
as  load  and  retisunce  facton)  and 's  not  tequhed  to  perform 
**probabtlistkt‘'  analytU.  As  summtrired  by  EUingwood  ei  ai  (2), 
the  principal  sdvanuges  of  probabUlitiu  Umit  state  design  tr«  (1) 
more  consistent  reliibiUty  (survivability}  U  attained  for  dinenau 
design  situations  becauM  the  dlfiutnt  variabilities  of  the  various 
swtctunl  strengihs  and  toads  are  treated;  (2)  the  reliability  level 
cm  be  chosen  to  lefioct  the  consequences  damage  on  structural 
collapse;  (3)  the  designer  aidilevts  a  better  undemanding  of  the 
fiindametual  sintctural  requliemenu  and  of  the  behavior  of  the 
sinicture  in  meeting  three  requirements;  (4)  the  design  process  is 
simpliSed  by  encouraging  th  >  same  design  philosophy  ind  proce- 
durre  to  be  adopted  for  all  materials  of  coosouction;  (5)  better 
Judgment  can  be  applied  to  nonrouiine  situations;  and  (6)  design 
manuals  can  be  updated  in  a  muio  rational  manner. 


Many  subdisciplines  of  structural  rellabilify  analysis  pro¬ 
vide  the  methods  for  developing  and  implementing  reliability- 
based  design.  These  include  sectmd  moment  methods,  stochastic 
load  combination,  stochastic  finite  elements,  and  design  code 
calibration.  Common  mathematical  toob  include  applied  statis¬ 
tics,  Monte  Carlo  simulation,  enor  propagation,  and  experimental 
design.  Optimixation  tlieory  has  also  been  applied  to  structural 
design  to  idenr'fy  designs  Cat  cptimlze  performance  (minimum 
cost,  maximum  reliability)  for  specified  design  aiteria. 

Constructioa 

Permanent  protected  structures  are  oftot  constructed  in- 
place  and  are  therefore  subject  to  construction  errors  and  toler¬ 
ances.  While  research  and  some  data  on  nonprotected  construc¬ 
tion  errors  and  uncertaintire  exist,  research  has  not  been  po- 
fonned  for  protective  constniction.  A  general  belief  b  that  these 
uncertaintire  are  small  compared  to  those  resulting  from  the  re- 
jqxMise  analysis.  It  b  noted  that  constructitm  enres  associated 
with  field  fortification  and  nonenghteered  sfructural  iqrgndes 
should  be  treated  separately. 


Test  add  Operation 

The  tMt  and  opmiion  phase  tepioenis  the  final  activity  of 
the  integrated  design  pnK«n  ami  sauetural  life  eyeb.  It  can  be 
divided  into  two  time  periods:  peacetime  and  wartime,  frt  peace¬ 
time,  controlled  expedmenti  and  tests  are  perfonned  to  assess 
smKXunI  peifonnnace.  ExperirneudetiipcaaoepunBdmeihods 
of  statistical  inference  am  useful  prolNMlbtio  rttcthm^ 
appUcationa.  In  wintoe.  the  sfruettma  wiU  be  subjected  to  ef¬ 
fects  from  enemy  weapons.  Mwy  souciures  will  be  damaged  to 
different  degrees.  Damage  aisresmena  and  analyses  of  safety 
will  need  to  be  performed,  both  it  •  field  level,  and  later,  in  a 
more  deuUed  cn^neering  analysis.  Useful  methodologire  are 
systm  identificatitm,  cngirreetiiig  databases,  attd  expett  systems 
(schnology. 


BURIED  mUCTtlRE  ANALYSIS  EXAMPLE 

Thb  example  competre  determittbtic  and  probabiUstic 
analyses  of  a  box-sh^wd  structure  buried  in  sandy  clay.  The 
threat  is  givre  and  Urn  damsge  piobeHlity  u  to  be  deteiinbed  as  a 
h»iait«(>ficalcd«andoIf(X)for(>.5<A<3.0ft/lb“.  The 
i'ructure  contains  critical  shock  sensitive  commutucations  equlp- 
meau,  which  has  a  frequency  of  10  Mx  in  the  domhuat  reiqwnse 
mode  and  b  hard  mounted  to  the  floor.  The  siruciuro  is  asnuned 
to  fail  to  perfotm  Us  mission  if  the  equipment  b  damaged. 

Detarmielatk  Asratysb  Meihodotety 

The  scope  of  tire  tictennuustis  analytis  metiads  (im<  thb 
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eMoqde  is  limited  to  existing  ilesignmaiiutlqrpe  methods.  Two 
CUhim  modes  ire  evaluited:  will  flexuil  Mm  isd  excessive 
floor  sccelentions. 

WaUFtuaulPattun.  The  deterministic  method  for  com* 
puting  die  flexuil  response  of  die  buried  will  is  the  equivslent 
singleHl^rBe«of>firMdom  sGractute>medium  interictjon  (SMI) 
model(3].  Hiesolutioaoftheeqittttioaofmodoagivesthedis- 
plscement  of  the  center  of  the  will.  In  Shis  ipproech,  die  peik  flee 
field  stress  is  multiplied  by  in  equivslent  uniform  loed  futor 
vveo  by  Kiger  [4].  Utis  bctor  sccounts  for  the  nonunifotm  ns- 
tnre  of  the  loads  due  to  convcadoasl  explosives  tnd  is  s  fttaction 
of  die  wan  a^wctislio  wl  the  distance  between  the  bomb  lad  the 
wall 


Hguie  1  shows  the  center  deflecdon  of  the  saucnue  wall 
as  a  Affledoa  of  scaled  stindoff  (H)  predicted  by  the  SMI  model 
Also  shown  in  the  figure  is  the  empiricil  NatUnal  Defease  Re- 
seardi  Committee  (NDRQ  damage  curve  for  buried  reinfoieed 
concrete  walls.  Note  that  the  NDRC  curve  is  for  walU  with  lengdi 
todilc)a»ssistiosGAi)iniherinteof5:l  iol?;l  indficedimen- 
iIdos  ratio  (a/b)  of  about  3:5.  The  exsntpte  problem  wall,  how¬ 
ever,  has  an  lAt  ratio  of  apfminutely  2J:l  and  a  face  dtmtasiona 
tadoofqiiinutimattfyoj^i.  Htnce,  the  flexual  model  uaed  in  the 
SMI  snaly^  it  cxpe^  u  be  conservative  fo  predkdnf  damage 
for  dtis  veiy  thick  wall  Rased  tm  damage  level  as  a  (Unction  of 
f^ppoK  rtxation  [52,  tto  SMI  modal  pm^  that  wiU  damage  is 
laes  than  or  equal  fosH^  for  X>2S.  moderate  for  L7<kai2J, 
heavy  for  1.4  <  X  £  1.7,  m4  is  St  itKipum  Mure  for  X  S 
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HGUREl.  DetermteiMic  Ptedictfoo  of  Wail  Dcflectloit 
tyfi* 


FIfm  1 4howa  a  ahaip  inmsliioobst^^ 
ageMteefor  snsUchangeefaiX,  Thisladttaioihs«eepiittiof 
d»  daCciUitw  vetws  raaga  curt*,  wMch 
fug  lUenutlkMi  of  petit  Cm  Said  straMWifo  tinge.  ThesMpneH 
of  fosottve  preset  t  ptoblstt  fotlm  an^ 
datctminiitfo  evaluation,  itts^tgastt  thatha  hi  able  fo  ptedimthat 
amall,  almoM  unooiiotHfot^  eha^  in  large 

chasiei  ht  predicted  damage.  Since  all  dta  hi  tha 

aaal)^  have  bean  n^tiectad.  diaie  is  iw  quMiiiutive  pto^^ 
smamsthetaargfabatweentheaemsa^ 

In-SMirAtiwShedi:  Tbefo-s&Pcturaabo^cnvinMmeiu 
hi  dkinttttnaed  by  aminutting  dtt  rei^^ 
the  strucnire  floor.  Ifltfo  tespamm  tputtn  fo  obtained  from 
mater  of  the  maxhtnun  floor  aceefatation,  velocUy.and  dieplace. 
mem.  u^  the  ptucedtm  hi  Refoteocea  6  and  7,  this  ex- 

taqile,  only  die  hcrittmul  amion  cautiiAf  from  dm  ipaciiied 

The  leapoiwi  Spectra  for  equipment  hard  tnouaM  to 
floor  fo  eadmaiad  by  scaling  dm  peak  floor  eeceteradoii,  velocity, 
md  ^rptoMwiwt  bf  2.0. 1.5,  and  1.2«  terpectively,  u  recom- 
meodadbyKifirlT).  theaefimomtivcdmamptiliedfflotioain 
datccifonaiomvefociiyitflddisptioemeMriitioMofthespecaa. 
R^iectivciy.  Tben«iliiii|aaidyshindicatctthmfarequ^^ 
c^pacityo{10g*s(foriahondttiatiotipulaetypicaioracot!tven- 
tinMlWMpc<(XmN»gdfaihOT  would  bapwdkte^ 


RsUal)iIii^4iaitd  inai^  bvolvw  dm  performs 
teqiHiPCf  of  eepf,  ftltwe  mode  oMlyi^f 

tkm  of  UMWtahitiae.  prior  tt  foe  pmM^^  Tbese 

iispaaredmicrtbedfodeMfoRefoito^  Oneoffoskeyotti- 
of  foif  pvfefow  ir  *hf*  q|'****'*'V*fd*w  of  patamjWr 

appmifhai  Fm  txfwplt.  Table  1  foe  resutis  tXf  dm 

mtoeKahuy  analyafo  for'foe  flaxmid  Mum  mode,  wbkh  am  ba^ 
on  befo  dkam  dam  «a4ytfo  andctmsparisansof  modal  pmii^^ 
and  etqmrimeeiA'Vfosfe  afpimpi^ 

MoemX^  fonuimion  was  uaed  to  im^pagaie  dieae  us* 
ceitafoilm  forcu^  dm  dmanmiiidiifo  smaeoM 
tbe  foUowini  patagraplm  sitttmirixe  dm  results  of  foe  saalyds 
for  esdi  ftfoM  mode  (wall  flexuie  and  fo-situemm  foodfl.  foe 
combM  fotiim  ptobiiiiUt> .  MdssntiiiviQi  analysis. 

WmB  Fktutv  Oeemgy.  tabfoZcoaqiaiteefoedeteaninis- 
tfo  (ttcniaafl  and  foe  imiabUity-beicd  piedictkiM  of  foe  deflectm 
orfoecaaieroffoewall  TheditforeeoMinfoeineanandMmfofl 
deflectient  art  due  to  foe  IS  percent  under  bias  (tee  t)hi 
foe  taoifiiaal  iBKidal  and  dm  fam  foat  dm  wad  tesponae  b  nonlinear 
Another  impanaht  obeervaUon  in  Table  2  is  dm  targe  atsndiid 
dmdationa  of  foe  reliability-based  aoalysia  rtnul^ 

Hfum  2  slmws  pccdicied  damege  as  a  fonmioo  of  scaled 
range  for  foe  deiMBHiiditic  and  ptobtWilinicaotlysii.  Notefoeso- 


61 


Tabte  1.  UncotiinQ'  Analysis  for  Buried  Structure  Example 


l&ceitiintyType. 

Model  and  Pvametos 

Uncertainty  Cliaracterization 

Mean/ 

Nominal 

Coefficient 
of  Variation 

1.  Prediction  Error 

Wait  Flexure  Model 

Max  Deflecdcn 

1.15 

0.25 

L/d  Outside  Dau 

1.0 

0.15 

Limited  Dau 

1.0 

0.05 

Groundshock  Model 

Peak  Free-Field  Stress 

1.0 

0.65 

2.  Random  Uncertainties 
Structure 

Concrete  Density 

1.0 

0.05 

Concrete,  E, 

1.0 

0.10 

Concrete,/^' 

1.4 

0.20 

SteeL/ 

1.5 

0.15 

Son 

IM  Weight 

1.0 

0.05 

Seismic  Velocity 

1.0 

0.30 

Free-Field  Stress 

Stress-Decay,  a 

1,0 

0.10 

Attenuation  Coefficient.  N 

Treated  in  Model  Pred.  Error 

Yield.  W 

Deterministic 

Range,  R 

Deterministis  (Par.  Variation) 

Damage  Criteria 

1.0 

0.20 

(Support  Rotation) 

Table  2.  Buried  Wall  Center  Deflection  Statistics 


Deflection(in) 

Range 

Deterministic 

Probabilistic 

(it/lo 

standara 

(Nominal) 

Mean 

Deviation 

1.0 

-• 

.« 

• 

1.5 

7.9 

12.6 

20.7 

2,0 

3.6 

5.0 

7.6 

2.5 

1.4 

2.2 

3.9 

3.0 

0.4 

1.0 

2,3 

*WtUCoUi|)ied 


called  "cookie-cutter”  shape  of  the  determinisdc  qiproach.  This 
sensitivity  presents  a  problem  to  a  designer  in  light  of  the  obvious 
uncertainties  associated  with  the  predictions  of  structural  re¬ 
sponse.  The  probabilistic  ^tproach  presents  a  more  realistic  pic¬ 
ture  in  terms  of  a  continuous,  monotonically  decreasing  curve  of 
damage  probability  versus  scaled  range. 


HOURE2.  Comparison  of  Mkted  Wall  Response 


IH-Struc&trt  Shock,  llie  probabilistic  analysis  of  in¬ 
structure  shock  omsidered  uncenainty  in  the  firee  Add  ground 
shock,  motion  of  the  structure  floor,  dynamic  properties  of  the 
equipment,  and  response  of  the  equipment  [1],  The  results  are 
compared  in  Table  3  to  the  deterministic  analysis.  In  contrast  to 
the  wall  response  failure  mode,  there  is  voy  little  difference  in 
themnninal  and  mean  values. 

The  predicted  in-structure  shock  accelerations  can  be  used 
to  estimate  probability  of  equipment  failure.  For  this  example, 
the  uncertainty  in  the  capacity  of  the  equipment  (fragility)  Is  mod¬ 
eled  using  a  lognormal  distribution  with  a  mean  capacity  of  lOg's 
and  a  coefflclent  of  variation  cf  25  percenL  Figure  3  shows  the 
predicted  probability  of  equipment  failure  as  a  function  of  scaled 
range. 

Combiiud  Faifun  Modts.  The  structure  mission  is  com¬ 
promised  if  either  the  waU  sustains  heavy  damage  or  a  critical 
piece  of  equipment  fails  by  In-structure  shock.  These  damage 
modes  have  b^  combined  exactly  within  the  Monte  Carlo  simu¬ 
lation.  Figure  3  shows  the  results  of  this  combined  failure  mode 
analysis.  The  probability  of  failure  for  the  combined  modes  is 
equal  to  or  sligiuly  higher  than  the  maximum  of  either  failure 
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Table  3.  Equipment  Reapmue  Statistics  Due  to  bi-Structure  Shock 


Defiection(m) 

Range 

Probabilistic 

(ft/lb 

Detenuinistic 

(Nomiiul) 

Mean 

Standard 

Deviation 

0,5 

14,0 

11.9 

1.0 

9.6 

9.7 

8.1 

1.5 

7.0 

.,1 

5.9 

2.0 

5.3 

5.3 

4.5 

15 

4.1 

4.2 

3.5 

3.0 

3.3 

3.3 

18 

mechanism.  IIib  difference  between  the  maximum  of  the  individ¬ 
ual  modes  and  the  combined  curve  is  not  large  because  only  two 
failure  modes  have  been  considered  and  these  are  correlated. 
Also,  the  wall  probability  of  failure  is  much  larger  than  the  equip¬ 
ment  probability  of  failure  (excqtt  at  the  higher  scale  ranges)  so 
that  wall  failure  dominates.  At  a  scaled  range  of  3.0,  the  two 
models  have  almost  equal  failure  probabilides  and  the  combined 
failuie  probability  is  approximately  double. 


FIOURE3.  bulividual  and  Combined  Failure  Probabilities 

VHuriaiHt^  Rtmkint.  Die  reliability-based  analysis  indi¬ 
cates  large  uncertaintii  -  about  the  pitdicied  mean  values.  Uncer¬ 
tainty  ranking  and  sensnivlty  analysis  are  useftil  byproducts  of  a 
pn>bid)ilisile  analysis.  example,  the  source  of  the  large  uncer¬ 
tainty  in  the  predicted  v,  i|  deflection  can  be  quantified  by  rank 
evaiuatiaii,basedaa|iartitir‘(0$i*<;i}.  Ihe  rank  analysis  was 


performed  by  )q)plying  stepwise  r^ression  analysis  procedures  to 
the  Monte  Carlo  simulation  ouqnus.  Ibis  analysis  indicates  that 
uncertainties  in  predicting  structural  wall  re^nse  are  dominated 
by  the  ground  ^ock  model  uncertainties  in  prediction  of  flee- 
field  stress  (t*  =  0.77).  Soil  seismic  velocity  (f  =  0.03)  and  SMI 
flexure  reqwnse  model  (i*  =  0.02)  uncertainties  contributed  to  the 
uncertainty,  but  were  much  less  important  Uncertainties  in  struc¬ 
tural  material  properties,  strengths,  etc.,  did  not  significantly  con¬ 
tribute  to  the  predicted  response  uncertainty. 


RELIABILITY-BASED  DESIGN  FACTORS 


Probability-based  analyses  are  not  always  expedient  for 
design,  particularly  when  the  structural  analyst  does  not  have  a 
background  in  probabilistic  methods.  For  this  reason,  nliability 
based  design  factors  (RBDF)  are  developed  as  a  research  product 
for  application  by  the  designer,  much  u  safety  factois  axe  used  in 
traditioail  design. 

The  RBDF  foimat  is  based  on  the  eoneqx  of  a  design 
factes',  which  is  the  ratio  of  the  required  nominal  capacity  to  the 
aomlNo/ structural  response.  The  use  of  the  tcnn/icmiHa/ implies 
that  the  capacity  and  ffi^xmse  are  computed  using  otdinaiy  deter- 
minisdeamiyais  procedures.  Hence,  the  designer  cootpuies  struc¬ 
tural  respense  (say,  for  example,  wall  ductilU:^  or  equipment  ac- 
ceteradon),  using  standard  methods,  compares  this  to  the  (Opacity 
(ductility  capacity  or  acceleration  capacity).  If  the  ratio  of  ctqMc- 
ity  to  lesponse  dme  not  ineet  or  excels  the  required  design  factor, 
then  the  design  is  modified  to  increase  the  ca^ty  or  reduce  the 
response.  The  requited  design  factor  is  selected  a  ubledtat 
fives  design  faettms  as  a  {ymctlon  of  the  ‘"'^al  pn^ility.  For 
example,  Table  4  fummatiz»  preliminary  design  fseton  devd- 
oped  in  the  Phase  I  research  effoit. 


TiR>le4.  DcsipFicton  for  In-Strucoae  Shock  and 
Buried  Wall  Deflection 


Desired  Reliabliity 
Goal  (Probability 
of  Survival) 

Design  Factor 

tn-Smic.  Shock' 

g 

1 

5 

0.75 

1.3 

4.9 

0.85 

1.8 

3.0 

0.90 

11 

4.0 

0.95 

18 

6.2 

*  IRequitadequiimesaGctkiaiiattfaptcltyiadeii^mf^^ 


**RequiKd  wtU  deficeiicn  captdty  xktifn  facto' I 
pt«ihae4  wall  ddUctlni 
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These  preliminaiy  factors  for  design  reliability  (sundvabiliQ^) 
goals  of  0.75, 0.85, 0.90,  and  0.95  range  from  about  1.3  to  3.0  for 
in-stnicture  shock  and  from  2  to  6  for  buried  wall  response.  Note 
that  the  factors  are  considerably  higher  for  the  buried  wall  design, 
quandtaiively  reflecting  ths  mudi  greater  uncertainty  in  predicting 
wall  deflection  response.  It  should  be  pointed  out  that,  in  addition 
to  design  application,  the  RBDF  tables  can  be  used  as  an  analysis 
tool.  For  a  given  structure  and  a  given  threat,  the  analyst  need 
only  perfoim  a  conventional  determinisdc  analysis  and  compute 
the  inherent  design  factor  (ratio  of  cqtacity  to  response).  By  al¬ 
tering  the  design  factor  table  with  this  computed  ratio  the  proba¬ 
bility  of  survival  is  obtained. 


CONCLUSIONS 

The  following  conclusiotu  are  based  on  the  results  sum* 
marized  herein  and  presented  in  detail  in  Reference  1: 

1.  The  potential  uses  of  probability-based  approaches  to 
the  analysis  and  design  of  protective  structures  have  been 
identified  and  the  advanuges  illustrated. 

2.  There  is  a  significant  amount  of  data  in  many  of  the 

key  areas  that  has  not  been  systematically  analyzed  Cor  un¬ 
certainties  and  model  etion.  Preliminary  analysis  of  some 
of  these  dsia  and  comparisons  to  design  methods  indicates 
that  nominal  design  are  not  always  conservative;  (P.  < 

0  JO)  for  the  oues  cotuddeted. 

3.  Fundamental  improvementa  to  the  analysis  and  design 
of  Air  Force  facilities  can  be  achieved  through  research 
and  application  of  modem  concqtts  of  structural  reliabll- 
hy.  These  iinprovemenij  would  result  from: 

(a)  development  of  a  reliabllity-bised  design  method¬ 
ology  that  vvould  put  all  protective  structure  desip 
on  a  consistent,  cost-effective,  and  a  balanced 
basis,  ind  would  biing  protective  structure  desip 
up-to-date  with  conventional  desip; 

(b)  identification  of  key  uncertainties  in  protective 
sttuoiute  analysis  and  desip  that  would  allow  a 
prioritization  of  important  research  areas: 

(c)  development  of  a  systems  analysis  tool  for  struc¬ 
tures  and  facilities  to  develop  improved  concepts 
of  survivabte  smicturcs  for  easting  threats. 
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EXTEDIENT  METHODS  OF  PROTECTION  TO  MITIGATE  STRUCTURAL  DAMAGE  AND  SPALL 
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ABSTRACT 

A  series  of  full-scsle  tests  using  generel 
purpose  bombs  wes  conducted  on  e  herdened  reinforced 
oonorete  sboveground  structure.  Seversl  expedient 
sMthods  of  nltlgsting  structursl  dsnsge  due  to 
eirblsst  end  fregment  Inpsots  were  tested,  and  the 
response  of  the  wells  of  the  structure  was 
■oni Cored. 

Results  from  these  tests  Indlcete  that  bare 
reinforced  concrete  walls  will  be  significantly 
damaged  by  the  given  weapon.  On  the  exterior 
surfeoa,  several  inches  of  concrete  cover  will  be 
destroyed  by  fragment  Impaocs,  and  on  the  Interior 
surface,  hlgh’Velocity  spall  will  occur.  This 
damage  can  be  pravanted  by  the  use  of  any  number  of 
expedient  aethoda  outlined  in  this  paper.  The  most 
eesC*effeoclve  method  appears  to  be  a  soil  berm. 

BACXSRiOUND 

The  U.  S.  Air  Force,  Europe  (USAFE)  is 
xaiponaibla  for  Che  design  and  construotlon  of 
miliCacy  facilities  which  meat  the  NATO  samihardened 
stnustura  re<|uireaencs  for  protection  sgeinst 
conventional  weapons .  NATO  re^iraments  and  current 
design  procedures  have  led  to  the  use  of  heavily 
reinforced  29.5 -inch  thick  concrata  walla  to  resist 
these  leads.  However,  teats  conduotsd  in  1980-1962 
by  the  German  Bundaswehr  Infrastructure  Staff 
(Rsfersncea  1  and  2)  Indicatsd  that  these  walls  were 
typically  over-reinforced  and  chat  more  eoonomloal 
designs  with  less  steal  reinforcement  would  provide 
the  required  protection. 

A  series  of  NATO  Samiherdaned  Design  Criteria 
Full  'Scala  tests  ware  conducted  by  the  Air  Force 
Snsinasring  and  Services  Center  (Af^C)  and  the  U.S. 
Army  Usterways  Experiment  Station  (UES)  during  the 
auamar  of  1987  (Reference  5).  The  objeotive  of  thia 
teet  aariee  wae  to  determine  the  effacta  of 
revetments  and  eand  beteui  on  the  raeponae  of 
reinforced  oonereee  wells  eubjeoted  to  the  nearby 
detonation  cC  e  cased  munition. 

TEST  DfiSCRiniON 

The  fuil-seale  sortas  wae  composed  of  sevaral 
teste  with  e  nearby  detonstion  of  a  genarsl  puepoae 
bomb  located  at  varLoua  positiona  around  the  test 
structure.  The  bomba  v<era  placed  at  the  seme  range 
from  the  teat  structure  for  each  of  the  teat  eventa 
.discussed  in  this  paper.  The  teat  structure  wae  of 


reinforced  concrete  coitscructlon,  approxlnstely  SO 
feet  wide  by  60  feet  long.  All  of  the  tests  were 
Instrumented  to  record  slrblasc  pressure  loading  on 
Che  wells,  In-structuce  accsleretlon,  and  wall 
dsfiseclon  versus  cine. 

Each  of  Che  wells  discussed  in  this  paper  were 
2S.3-lnch  thick  sections  with  0.25  percent  stasl 
relnforcsmsnt,  One  bars  vail  was  tasted  «s  s 
bassllns;  additional  pcotacclve  m<»ssurss  testad 
included  s  sand  bans,  s  sacrificial  prs-oast  panel, 
ssnd-grld  ravecmsnts,  and  s  portable  Blcburg 
cevsement. 

Each  of  Che  2S,S-lnch  walls  had  0.25  paresnt 
vscttosl  atsal  rsinforosnsnc  in  each  face  oonslsclng 
of  Numbsc  5  bars  apscad  at  5  inches.  Horlsontal 
ratnforcanant  oonsistsd  of  Numbar  5  bars  sptosd  at 
4.3  inches.  Stitgla-lag  stirrups  fabrlestad  from 
Number  3  bars  provided  ahasr  rslnforcsmsnc,  The 
stirrups  wart  apscad  at  IS  Inches  on  each  vsrcloai 
bar,  and  acaggerad  at  9  Inchss  off-center  every 
other  bar.  Each  of  the  walls  had  a  olssr  span  of 
about  13  fsst-2  inohas.  and  e  eonocacs  cover  of 
about  3  l(,chsa  on  ths  sxcsrtor  faos  and  3/4  Inch  on 
the  interior  face. 

Several  concepts  dsslgnsd  to  attenuact  Che 
elrblesc  end  ulntmlse  spalling  were  considered  In 
the  full-scale  teet:  (1)  e  4-foot-7-lnoh  tend  berm 
with  a  2:3  slope;  (2)  s  S-lnch  chick  llghrly 
reinforced  pre-cast  panel ,  anchored  to  the  structure 
with  removable  bolts  at  ths  cop  and  with  a  small  gap 
bsewssn  panel  and  structure ;  (3)  ssnd-gtid 
ravetmsnts;  (4)  porcsbla  Bliburg  reveemsnea. 

1.  Sand  Berm 

Previous  half- tesla  casts  (Rsfsrsnee  4)  have 
Ikultcstsd  Chat  the  addition  of  a  barm  to  the 
excsrlor  of  ths  structurs  will  reduce  midepsn 
deflections  and  virtually  sliminata  spalling.  A 
bans  of  silty  sand  (taken  from  the  test  site)  with 
e  slope  of  2'3  was  plaead  against  one  25.5  inch 
wall.  Ths  bstm  sxtsndad  only  part  way  up  the  wall, 
since  previous  esses  had  Indicatsd  ehsc  the  most 
severs  blast  loads  and  fragment  damage  occurred  neer 
the  ground  eurface.  .  The  berm  was  about  4  fast- 7 
tnchea  high. 

2.  Pre-Ceet  Panels 

The  Carman  Bundaswshr  Infrastructure  Staff 
had  previously  conductsd  half-seals  tests  of  4-  and 
6- inch  prs-esst  panels  with  favorable  results.  In 
the  German  tssts,  the  pre-cast  panels  wars 
complatsly  destroyed,  but  they  did  ptotset  ths  main 
attuccura  walla  from  fragment  damage  and  alimlnsted 
spall.  The  panel  used  in  the  full-aeala  eertes  wes 
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6  inehas  thick. 

3 .  Sand  Grids 

Sand-grids  vara  devalopad  at  WES  as  a  rapid 
■sans  of  constructing  roadways  on  baachas  for 
■llltary  vehlclas  (Rafaranca  S).  Thesa  sand-grids 
hava  slnca  basn  taatad  as  ravatnants  for  artillery 
anplacaaants  with  good  rasults  (Reference  6) . 
Sand-grids  consist  of  8-lnch  high  hlgh-danslty 
polyathylana  strips  connected  by  ultrasonic  or  heat 
welds .  They  are  manufactured  and  shipped  In 
collapsed  4- Inch  thick  sections  that  expand  to  20 
fast  during  construction.  Each  expanded  grid 
section  is  either  4  or  8  feat  wide  by  20  feat  long 
and  contains  a  honeyeonb  arrangaoant  of  calls. 

4.  Bltburg  Revatmants 

The  Bltburg  portable  revetments  are  modular 
reinforced  concrete  sections  Intended  to  protect 
Items  lass  than  6  feet  tall  (seo  Figure  1).  Each 
section  was  about  6  feat-7  Inches  wide  and  6  feec-7 
Inches  high. 

RESULTS 

Although  each  of  the  test  events  was  unique,  a 
few  general  obsarvetlons  can  be  aade.  For  all 
aboveground  tests  against  unprotected  walls,  damage 
to  the  exterior  of  the  structure  was  similar.  The 
cover  concrete  on  the  exterior  face  of  the  lower 
tialf  of  the  wall  was  eroded  away  by  the  fragments, 
exposing  the  reinforcing  steel.  Spalling  on  the 
Interior  face  of  the  unprotected  walls  was  limit. d 
to  the  depth  of  the  relnforoamant,  and  was  greatest 
In  the  lower  one- third  of  the  wall  span.  Ko  spall 
was  aasoolated  with  the  pre-cast  panels  or  sand 
bem,  Mid  there  was  very  little  sr .  d^ge  with 
the  Bltburg  revetments  or  sand-grid  revecmanti. 
Direct  comparisons  of  peak  pressure  measurements 
made  on  the  exterior  face  of  ears  wall  are  given  In 
Table  1,  and  wall  response  f  each  event  Is  given 
In  Table  2.  A  detailed  (bsorlptlon  of  all  test 
results  follows. 

1.  Pre-Cast  Panels 

The  lower  half  of  each  of  the  penele  wee 
completely  destroyed  during  the  test  svent,  while 
the  upper  halves  wsrs  left  Intact,  but  fall  to  ths 
ground.  Tha  panda  pravantad  all  but  a  faw 
fragisantB  from  panetreclng  the  wall  of  tha  atrueeura 
<Sae  Flgurt  2).  tlie  avaraga  paok  praeaura  behind 
tha  panela)  wee  about  330  pat,  oompared  to  peak 
praasuras  of  2000-3000  psl  at  elmllar  rangat  on  an 
unprotacttd  wall.  The  peek  eoceletetlon  of  the 
wall,  meaeurad  3  faet  from  the  floor  surface  and 
oppoerte  the  weapon,  was  sbout  8900  g's,  end  the 
peek  meaaured  deflection  wee  about  1.8  Inches.  The 
interior  of  the  test  well  suscained  minor  oraoklng, 
but  had  no  spall  damage. 

2.  Bltburg  Revetments 

The  Bltburg  revetsMnts  wets  originally  designed 
to  protect  fuel  storage  tanks  snd  ochsr  Itsms  lass 
then  6  feat  tell  from  fregmentatlon  end  alrblast 
effects.  Tha  Bltburg  design  wee  eeleeted  for 
testing  because  of  Its  reedy  evellablltcy  on  US  Air 
Force  besae  and  because  of  Its  portability.  The 
revetments  were  placed  directly  sgelnat  the  vail  of 
tha  structure,  with  about  8  Inchei  clearance  between 
tha  well  and  the  vertical  section  of  the  revetments. 
Ths  rsvetmsnts  directly  In  front  of  the  threat 
wespon  were  almost  oomplstely  dsstroyed  by  s 


combination  of  alrblast  and  fragmentation  effects 
(Figure  3),  but  the  wall  behind  the  --evetments  was 
largely  undamaged.  As  In  previous  test  events, 
there  was  extensive  concrete  cratering  on  the 
unprotected  portion  of  the  wall.  The  Interior  of 
Che  test  wall  sustained  minor  cracking  In  addition 
to  slight  spall  damage  (Figure  4) . 

The  peak  alrblast  measurements  from  this  event 
Illustrate  the  benefit  of  uslr?  revetments  Co  shield 
a  structure  from  alrblast:  the  peak  pressure 
measured  3  feet,  4  Inches  from  Che  ground  surface 
and  directly  opposite  the  weapon  was  about  232  psl, 
while  the  peak  pressure  o  feet,  7  Inches  above  Che 
ground  surface  (Just  above  Che  cop  of  Che 
revetments)  was  about  2270  psl.  The  peak 
acceleration  of  the  wall,  measured  5  feet  from  the 
floor  surface  and  ooposlce  the  weapon,  was  about  900 
g's,  and  Che  peak  measured  deflection  was  about  1.4 
Inches . 

3 .  Sand-Grid  Revetments 

The  sand- grid  configuration  tested  consisted  of 
twelve  layers  of  4-fooc  wide  sand-grid  sections 
(total  height  about  8  feet),  placed  with  a  small 
clearance  (about  4  inches)  between  the  revetment  and 
the  structure.  Construction  of  the  revetment  Is 
shown  In  Figure  5.  The  sand-grid,  as  expected,  was 
badly  domagad  by  the  test  event  (Figure  6) ,  but 
allowed  no  new  fragment  panecracions  Into  tha 
structura  balow  about  10  faet,  and  significantly 
attanuatad  the  peak  pressures  msssurad  on  tha 
exterior  of  tha  atruetura  wall.  The  peak  alrblaat 
maasuramanta  on  tha  structura  and  behind  the 
revetment  ranged  from  about  109  pal  to  297  pal.  Tha 
peak  alrblaat  maaauramant  on  tha  wall  was  about  2300 
Dsl,  occurring  Juac  balow  the  roof,  directly 
opposite  tha  weapon.  Tha  peak  acoalaratlon  of  the 
vail,  maasursd  5  feat  from  the  floor  surface  and 
oppoalte  tha  weapon,  was  about  170  g’t,  and  tha  peak 
■aaturad  daflacclon  was  sbout  0.7  Inchas.  The 
interior  of  tha  tost  wall  sustained  minor  cracking 
In  addition  to  slight  spall  damage, 

4 ,  Sand  Barm 

This  tset  event  wee  designed  to  evaluate  the 
performance  of  a  send  berm  In  eetenuaelng  the  damage 
to  the  exterior  of  the  structure.  Conaerucead  from 
a  silty  sand  taken  from  the  teat  alta,  tha  barm  was 
about  4  ftac-7  Ir.ehea  high  with  a  2:3  slope.  The 
berm  Is  shown  pretest  In  Figure  7. 

A  postaat  photograph  of  tha  berm  is  shown  in 
Figure  8.  Tha  peak  Interface  preesurea  measured  at 
the  interface  of  the  berm  and  structure  ranged  from 
74  pal  CO  312  pal  dlruccty  opposite  the  weapon.  The 
peak  alrblesc  meesuccmants  on  the  wall  end  above  the 
berm  ranged  from  326  psl  near  the  roof  of  the 
structura  to  3380  pal  directly  above  the  cop  of  the 
berm.  The  berm  proved  to  be  one  of  the  more 
affective  maaauree  evaluated  In  mlclgecing 
structural  dsmega.  Cratering  of  the  exterior  face 
of  Che  atrueeura  was  limited  to  thee  portion  of  the 
wall  above  Che  berm.  This  arse  could  easily  be 
atnlmtzad  by  Increasing  the  height  of  the  barm.  Aa 
shown  in  Figure  8,  an  added  advancags  to  the  usa  of 
soil  barms  is  thslr  inhsranc  multiple  strike 
capability;  most  of  the  msterlal  of  ths  berm  wee 
left  In  place.  Tha  peak  maeaured  aocslaraelon  of 
Che  wall  about  3  feat  from  tha  floor  surface  Was 
about  1730  g's.  thara  waa  no  spall  and  only  minor 
flexural  cracking  aaaocUcad  with  thla  east  event. 
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5.  B«s«Iln«  25. S  inch  U«ll 

This  test  event  was  designed  to  test  a  baseline 
unprotected  25.5  inch  wall.  For  this  test  event, 
the  weapon  was  kept  at  the  same  standoff  distance 
fron  the  structure  as  in  previous  events. 

As  expected,  this  test  resulted  in  nuch  more 
damage  to  ^e  exterior  of  the  stnteture  than  any  of 
the  previously  mentioned  tests  (Figure  9) .  Unlike 
any  of  the  previous  test  events,  this  event  resulted 
in  a  large  quantity  of  hlgh>spead  spall  fragments 
(Figure  10} .  The  interior  of  the  structure  was 
littered  with  spall  fragments  located  anywhere  from 
innodiately  inside  the  test  wall  to  the  other  side 
of  the  structure,  about  60  feet  distant.  The  peak 
pressures  measured  outside  the  structure  and 
immediately  opposite  the  weapon  ranged  from  400  psl 
to  over  2800  psi.  The  peak  acceleration  of  the 
wall,  measured  5  feet  from  the  floor  surface  and 
opposite  the  weapon,  was  about  21000  g's,  and  the 
P'tak  measured  deflection  was  about  2.5  inches. 

SUMHAB.Y 

Loading  on  the  exterior  face  of  the  structure 
from  the  detonation  of  the  criteria  threat  is  due  to 
a  combination  of  alrblast  and  fragmant  Impact.  For 
walls  unprotected  by  bams  or  revetments,  the  peak 
presaures  are  highest  near  the  ground  surface  and 
directly  opposite  the  weapon,  and  lower  on  ocher 
portions  of  the  wall.  The  airblasc  distribution  is 
highly  cranalent,  concentratad  near  the  bottom  of 
Che  wall  at  early  cine  after  detonation  and  decaying 
at  a  fast  race  to  a  lower  magnitude  later  in  time, 
The  greatest  density  of  fragmnt  liapacca  occurs  on 
the  lower  portion  of  the  wall.  These  impacts 
destroy  the  cover  of  concrots  In  this  area,  exposing 
the  exterior  reinforcing  steal  and  weakening  the 
wall. 

For  walla  with  berms  end  revetmanta,  the  peak 
pressures  are  highest  on  the  upper  portion  of  the 
wall  (just  above  the  barm  or  revetment),  and  are 
siplficantly  reduced  on  the  lower  portion  of  the 
wall  under  the  cover  of  the  berm  or  revetment.  Both 
the  average  blase  pressure  and  impulse  are 
significantly  lowar  for  ths  protected  valla.  The 
berms  and  ravatmante  are  also  very  affective  in 
stopping  fragmanta. 

All  axpadienc  methods  used  In  the  full*eoala 
test  to  reduce  spslletlon  were  effective.  The  use 
of  harms  altminsees  spsll  tneirsly  snd  is  the  moet 
cost  effective  solution.  The  portable  Bitburg 
revenancs  work  well  (for  a  single  hit  only),  and 
are  portable  and  easily  lepleced.  The  send  grid  is 
an  expedient  meant  of  protection  for  single  htce. 
The  polyethylene  grid  can  be  compressed  for  easy 
storage  prior  to  it's  use. 

Since  ths  primary  demsge  meeheniem  for  the  bare 
walla  was  apsllation,  Incraasing  the  atsel 
reinforcarnsne  from  0.25  percent  Is  Ineffective, 
unlaee  additional  shaar  rsinforcsmsnt  is  added  end 
tebar  spacing  is  minlmlzad.  Walls  protseced  by  s 
ban  or  raveemsne  psrformsd  ths  bast  due  mostly  to 
the  fact  that  ths  barm  or  revetment  reducae  the 
load.  They  respond  in  a  flexural  mods  with  very 
little  daf lection  end  only  minor  scniotural  derege. 


Th(j  re,^earch  was  sponsored  by  Che  U.S.  Air 
Force  I.ngineaiing  and  Services  Center,  Tyndall  AFB, 
FL,  under  the  direction  of  LTC  Bob  Majka. 
Perjilsslon  to  publish  this  paper  was  granted  by  Che 
Office,  Chief  of  Engineers. 
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Table  1.  Peak  Pressure  Comparisons 


Distance  up  Wall 

Peak  Pressure  Measurement,  psl 

Baseline* 

Pre-Cast 

Panels 

Bltburg 

Revetment 

Sand-Grid 

Revetment 

Sand  Berm 

0 

7035. 

385. 

714. 

125. 

512. 

2'-3" 

482. 

3' -4“ 

1720. 

370. 

232. 

150. 

6'-7" 

2565. 

359. 

2270. 

297. 

>3300. 

13'-2- 

1214. 

203. 

497. 

2500. 

806. 

*  Composite  from  several  test  events. 


Table  2.  Wall  Response  Comparisons 


Figure  4.  Interior  oi  structure  following  Bltburg 
revetnent  test 


Figure  6,  Damage  to  Sand-Grid  revetnent 


Figure  8.  Poscest  view  of  sand  berm 


Figure  10.  Interior  damage  to  baseline 


Figure  Damage  to  baeellue  test  wall 
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PaOTECTIVE  CONSTBUCTION  DESIGN  VALIDATION 

LT  COL  aOBEBT  J.  MAJKA  AMD  MR  WALTER  C.  BUCHKOLTZ 

HQ  XXR  FORCE  ENGIMBBRING  AMD  SERVICES  CENTER 


Ale  base  facilities  that  will  cesist 
the  effects  of  a  conventional 
weapons  system  ace  expensive. 
Elements  include  tons  of  conccete. 
tightly  spaced  heavy  celnCocclng 
bars,  and  costly  blast-cesistant  aic 
valves  which  pcever4t  explosive 
ovecpcesBUces  fcom  enteeing  the 
sheltee  through  heating  or  aic 
conditioning  systems. 

Recent  validation  of  the  North 
Atlantic  Treaty  Organisation  (NATO) 
aircraft  shelter  design,  conducted 
by  the  Ait  Force  at  Tyndall  afb. 
Pla.«  proved  that  less  costly 
sheltee  design  could  be  Just  as 
.effective. 

The  test  goal  was  to  determine  the 
optimum  combination  of  struotucal 
features  needed  to  protect  personnel 
and  equipment  while  reducing 
consteuctlon  costa. 


Background 

Readquactecs.  U.3.  Aic  Pocess. 
Bucops  (USAPE)  is  responsible  for 
the  design  and  construction  of 
military  CaciUties  that  comply  with 
NATO  protective  structure 

requirements. 

Tests  by  the  German  Oundesweht 
infrastructure  Staff  showad  that  tho 
typical  semlhacdoned  facility  was 
ovetceinforoed  and  that  a  more 
economical  design  was  appropriate. 
Shelters  are  designed  to  resist 
large  localised  loads  fcom 
conventional  weapon  detonations  on 
axtecior  vails. 

Two  considerations  not  satisfied 
in  the  German  concept  were  Interior 
wall  spall  and  defleerlon  of  walls. 
Spall  occurs  when  blast  pressure 
loading  or  the  inpacr  of  a  high 
energy  object  (such  as  a  large 
fragment)  cause  a  stress  wave  to 
travel  chCougU  a  wall  and  reflect 


off  the  intecloc  wall  surface.  The 
conccete  between  the  celnfotclng 
bars  and  the  surface  is  spalled  at 
velocities  high  enough  to  damage 
equipment  or  critically  injure 
personnel  inside  the  structure. 

Spalling  is  related  to  the  cate  of 
deflection.  Deflection  can  also 
displace  the  primary  shelter 
structure  fcom  its  foundations. 

Test  methods  were  significant 
because  they  were  done  in  two 
phases.  The  first  utilized  scale 
model  BtcuoturoB.  After  confidence 
limits  were  established  on  concepts 
being  costed,  a  full-scale  test  was 
used.  Host  earlier  shelter  testing 
used  costly  full-scale  shelters 
subjected  to  explosive  events.  This 
cestcicCsd  a  multiple  trial  approach. 

Values  galnsd  fcom  the  scale-model 
effort  were  then  used  Co  dsetgn  a 
full-scala  facility.  Research 
dollars  were  saved  by  incorporating 
these  results  into  one  full-scale 
test. 

The  Cull-soaie  test  offered  an 
opportunity  to  evaluato  several 
protective  blast  valves  and  doors  at 
Che  same  time.  These  were 
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previously  evaluated  only  against 
blast  pressures  generated  in  a 
laboratory  and  not  tested  against 
fragment  and  blast  loadings 
generated  by  conventional  weapons. 

All  constructign  was  accomplished 
in-house  by  the  HQ  AFBSC  Operations 
Support  Branch.  further  saving 
money.  support  in  design, 
instrumentation  and  photography  was 
received  from  the  corps  of 
Engineers.  Waterways  Experiment 
Station  (WES).  Vicksburg.  Miss. 

The  first  scaled  test  series 
addressed  the  deflection  problem. 
Walls  were  constructed  with  various 
percents  of  principal  steel 
reinforcement.  Also  evaluated  to 
determine  the  bast  shear 
reinforcement  and  resistance  to 
deflection  were  the  performance  of 
stirrups  (concrete  shoes  which  hold 
walls  in  place)  and  dowels  (large 
pins  which  anchor  overstructures  to 
their  foundations). 

The  second  scaled  test  series 
evaluated  various  ways  to  prevent 
spall. 

All  Costs  were  conducted  at  the  HQ 
AFBSC  Facilities  and  Pavements  Test 
Sits. 

Subscale  Testing 

Half-scale,  rsinfocoed  concrete 
box  steuctures  were  used  in  both 
scaled  test  series.  Bach  was 
designed  to  model  vail,  coof  and 
floor  slab  eeocione  of  a  typical 
aamihacdened  facility.  Scaled, 
caecd  cheeges  at  a  specified  chreet 
distance  weee  ueed  to  simulate  a 
oonventionai  weapon  detonation. 

Fourteen  ecaied  teste  were 
conducted  in  two  teeies  (six  in 
Secies  1  and  eight  in  Seciae  tl). 
In  Secies  1.  vacious  peecsncagss  of 
Che  Bccuctueal  volume  occupied  by 
principal  eceel  including 

high-denstcy.  medium  low  density, 
and  two  levels  of  low  density  were 
ueed  in  each  face.  (The  cucrent 
design  tsguirss  ons  of  the 
high-denslty  standecds).  The  floor 
and  toot  designs  were  consistent 
with  current  acmihatdened  deeign 
eeguireaence.  TOec  eectione  were 
placed  against  an  L-shapsd  reaction 
Bcruocuee  so  that  woveiasnc  was 
mloimlssd  (Figure  1).  Paec  testing 
indicated  that  shear  failure  in  the 
wall  was  possible,  single- leg  shear 
etlccupt  were  ueed  in  tour  tests, 
while  dowels  were  used  in  two  tests. 

The  Setioa  ll  ccets  evaluated 
three  methode  for  controlling  tpail 
(Figure  2):  earth  berms,  intscioc 
spall  placet,  and  thicker  walla. 


Berms  were  constructed  from  silty 
sand  with  slopes  of  1:1.5  extending 
partially  up  the  walls. 

Series  I  tests  pcovlded  blast  and 
fragment  patterns  and  established 
critical  heights  for  the  protective 
devices.  steel  spall  membranes  were 
Installed  halfway  up  on  the  interior 
low-thickness  wall.  Thicker  walls 
were  also  tested.  Because  of 
reinforcement,  the  thicket  wall  had 
the  same  flexural  capacity  as  one 
with  a  high-percentage  of  steel 
celnf occement  and  low  thickness. 

Full-Scale  Facility  Tasting 

A  full-scale  facility  (15ta  wide  x 
lam  long,  x  7.5m  high)  was  completed 
in  April,  1987.  A  partial  penthouse 
along  the  east  side  (2.7m  high  by 
3.7m  wide  by  18M  long)  contained 
various  blast  valves. 

One  set  of  blast  valves  was 
connected  to  an  operational  test 
ventilation  system. 


The  blast  valvee  were  placed  at 
three  locations  to  test  their 
pectocBsnce  under  difCerent 

conditions  (see  Figure  3  on  next 
page)  along  the  east  wall  of  the 
penthouse,  in  a  protective  structure 
at  Che  north  end  of  the  penthouse, 
and  behind  a  blast  wail  along  the 
west  wall  of  the  facility.  This  was 
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thtt  lowest  location  and  subjected 
the  valves  to  the  most  ccitlcal 
enviconmenti  Including  bomb 
fcagments,  dust,  soil,  and  conccete 
fragments  from  the  facility 
components. 


The  south,  east  and  west  walls 
were  65  cm  thlcK  with  high  principal 
steel  reinforcement  at  both  faces. 
Because  the  north  wall  was  80  cm 
thick  with  medium  principal  steel  at 
both  faces.  It  had  the  same  flexural 
capacity  as  the  65  cm  wall. 


The  five  blast  doors  tested 
Included  two  from  WES,  and  one  each 
from  Luwa  Ltd.  (Switzerland),  Temet 
USA  INC. ,  and  the  United  Kingdom 
Public  Services  Agency,  manufactured 
by  Energy  Equipment  Co.  Four  blast 
valve  designs  were  tested  from  Luwa 
Corp.,  J.P.  Sheltec  (Sweden),  Bately 
Valve  Co.  (England),  and  Temet  USA 
Inc. 

The  north  half  of  the  facility 
also  had  a  basement  constructed 
underneath  It  to  be  tested  against 
burled  weapons.  Motion  data,  using 
standard  office  furniture  and 
Instrumented  mannequins,  were 
recorded  principally  in  the  basement 
and  on  the  first  floor. 
In-structure  shock  or  acceleration 
of  the  floor  was  recorded  to 
determine  equipment  response  and 
operability.  A  backup  power 
generator  bolted  directly  to  the 
floor,  with  no  shock  isloation 

equipment,  was  operational  during 
Che  fleet  two  tests.  Spall  plates 
inside  of  Che  eaet  wall  were  welded 
Co  a  steel  frame  and  anchored  along 
the  calling  and  floor  by  angle 

iron.  Anchors  extended  from  the 
places  and  wart  cast  in  place  when 
Che  concrete  walls  ware  pourad. 

Several  other  protective  syetens 
were  tested  to  determine  their 

ecfecclveneee  in  preventing  spall. 
Full-height  precast  concrete  panels 
were  bolted  to  the  south  wall,  with 
3-lnah  air  gape  between  the  wall  and 
tha  panel.  Btcburg  csvetmente 
approximately  3  m  high  were  also 
tested  (see  riguce  4).  A  sand  berm 
(see  Figure  $  on  next  page) 

eonstrucced  along  the  west  wall  wae 
teeced  for  epall  protection,  ae  wae 
a  send  grid  eyetem  along  the  south 
wall.  Beseembllng  a  corrugated 
pepet  packing  carton  bottle 
protector,  when  pulled  apart,  the 
grid  material  forme  a  honeycomb 
partem.  sand  is  poured  Into  the 
honeycomb  openings  until  they  ace 
filled,  after  which  new  sand  and 
honeycomb  layers  ace  completed  until 
the  eand  grid  is  approximately  2  m 
high. 

Chamber  cevetmente  were  also 
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t«st«d  against  a  differanc  weapons 
Cheeat>  because  they  wece  located 
facthee  £com  the  detonation.  These 
cevetments  can  be  tasked  to  3  m 
height  and  connected  in4»ide  with 
steel  stcaps  and  filled  with  sand. 
They  ate  made  Ccoo  standard  concrete 
and  arc  celnfocced  with  steel  mesh, 
steel  fiber  or  polypropylene  fiber. 
The  mesh  and  fibers  wece  also 
combined  for  reinforcement  in 
several  chambers  (Figure  6). 

Bguipment  Test  Results 

The  half-scale  tests  provided 
valuable  installation  data.  '.'he 
anchors,  steel  frame  and  angl<$  iron 
installation  wece  the  keys,  to  the 
successful  performance, 

All  blast  doors  received  damage 
in  different  amounts  from  the  bomb 
fragments.  Locking  devices,  hinges 
and  door  frames  were  weak.  Several 
doors  were  subjected  to  nulciple 
detonations.  Outer  steel  surfaces 
wece  damaged  but  no  doors  were 
penetrated  completely  bv  the  bomb 
fcagmente.  All  dooic  provided 
adeguace  blast  protection,  although 
one  door  fcame,  damaged  by 
Ccagments,  was  breached  and  intncUc 
facility  damage  occuced. 

All  blast  valvts  functioned 
adequately  when  first  tested, 
Location  of  ths  valvss  was  critical 
bseause  dust,  soil  and  c-onceets 
fcagmsnts  affectsd  thsir  opetation. 
After  tepsatsd  testing.  several 
valves  failed  to  open.  The  valves 
acs  normally  Instalisd  in  banks  of 
ssvecai  valves  each.  Overall 
peefotmanee  of  the  valves  is  still 
being  snaiysod. 

the  ventiUtlon  system  and  power 
genetacoc  functioned  properly  during 
the  testa,  with  no  damage  to  systesi 
components  Ths  spall  plates 
functioned  extremely  well.  with 
minimum  wall  defleccioo  and  nc 
spalling. 

Structural  tsmt  Results 

the  facility  systems  parforitnd 
extcsmoly  well  against  <ip>iclCi»d 
wsapens  threats.  theca  ace  some 
additional  conclusions:  Rlast  dcots 
must  not  bo  In  the  weapon 
line^ot-sight  or  eevace  fragment 
damage  may  result.  Blast  valves 
must  be  pUcnd  above  the  facility. 
In  a  protective  enclosure.  to 
ainimlte  or  eliminate  bla«it 
pressure,  bomb  fragment  daa^ge,  and 
damage  from  soil  or  concrete 


fragments. 

Thick  walls  with  sheer  stirrups 
and  a  reduced  percent  of  reln'tocclng 
steel  worked  extcaraely  well,  with 
deflection  low.  The  unpcocsctsd 
wall  experienced  spall  problems,  but 
when  protected  by  precast  panels, 
Bitbucg  cevetments,  sand  grid  oc 
sand  berm,  no  spall  was  observed. 
The  sand  berm  provided  ths  best 
overall  protection. 

The  chamber  revetments  combining 
steel  mesh  and  fiber  reinforcement 
(steel  oc  polypropylene)  wece 
effective.  Although  sections  wece 
subjected  to  six  detonations,  no 
Ccagments  penetrated  the  rear 
surfaces,  demonstrating  excellent 
protection  for  equipment  or 
facilities  against  fragment  damage. 
The  -  thickest  wall  provided  no 
additional  protection  over  the  less 
thick  wall.  therefore  additional 
consctuction  costs  ace  not 
justified.  Deflection  was  minimal 
but  spall  still  occuced. 
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SHOCK  ISOlATXOfl  FG.I  COm’EHTIONAL  WEAPONS 


S.  A.  Kigar,  J.  H.  Weachersby,  D.  W.  Hyde 


U.S.  Army  Engineer  Waterways  Experiment  Station 
Vicksburg,  Mississippi 


ABSTRACT 

A  series  of  full>scale  tests  using  general 
purpose  bombs  was  conducted  on  a  hardened  reinforced 
concrete  aboveground  structure  with  a  partial 
basement.  Weapons  were  placed  in  both  above-  and 
below-ground  configurations,  and  the  response  of 
various  types  of  equipment,  computers,  and 
instrumented  mannequins  were  monitored. 

Results  from  these  tests  indicate  that  any 
equipment  not  restrained  against  movement,  i.e. 
.<)imply  sitting  on  Che  floo.;  or  isolation  pads,  will 
undergo  slgnlficart  rigid  body  displacements.  In 
general  the  best  isolator  performance  was  obtained 
from  the  cupnounc  series  neoprene  isolators.  Results 
also  indicate  that  the  only  real  danger  of  injury 
(from  shock)  to  personnel  is  from  falling  eq’  Ipment, 
e.g.  bookshelves,  cabinets,  or  light  fixtures. 


BACKGROUND 

Shock  lsolaci''n  for  equipment  in  blast-resistant 
structures  can  be  a  costly  and  uncertain  procedure 
for  conventional  weapon  threats,  especially  for 
aboveground  structures  with  combined  airblast  and 
fragment  loading.  The  in- structure  shock  environment 
is  known  only  approximately  at  best.  Methods  for 
calculate...''  in-structurw  shock  for  aboveground 
structures  generally  assume  plane  wave  loading  on  the 
struct  re,  instead  of  the  highly  localized  transient 
loads  gat, 'Mated  by  conventional  weapons,  and  they 
typica'ly  Ignore  fragment  loads.  Also,  there  are 
almost  no  data  on  modem  communication  and  computer 
ecuipmenr  fragility.  This  usually  leads  to  an  over- 
designed  and  expensive  shock- isolation  system  whose 
performance  Is  uncertain. 

A  series  of  NATO  Semihardened  Design  Criteria 
Full-Scale  tests  was  conductsd  by  the  Air  Force 
Engineering  end  Services  Center  (AFESC)  and  the  U.S. 
Army  Waterways  Experiment  Station  (WES)  during  the 
summer  of  1987  (Reference  1).  During  August  1988, 
the  WES  conducted  a  second  sarlss  of  tests  on  the 
seme  full-scale  structure  (Reference  8).  The  1988 
tests  wars  sponsored  by  the  Defense  Nuclear  Agency 
(DKA)  xn  cooperation  with  AFESC  to  further  evaluate 
in-strueture  shock  and  shock  isolation  mathodc. 
Equipment  tested  in  this  second  series  included  a 
oospater,  several  pieces  of  slectronio  equipment,  a 
larga  air-handling  unit,  and  a  l,000*kw  generator. 


Shock- isolation  methods  were  evaluated  by  repeating 
tests  under  several  different  isolation  conditions, 
including  a  baseline  test  with  the  equipment  hard- 
mounted. 

In  the  1987  test  series  all  equipment  was  bolted 
directly  to  the  floor  in  a  hard-mounted 
configuration,  or  simply  unattached,  as  with  some 
desks  and  bookshelves.  Data  from  these  tests 
indicate  relatively  high  acceleration  levels  within 
the  structure,  but  the  generator  and  air-handling 
equipment  remained  operational  before,  during,  and 
after  the  tests.  In  Figure  1  the  shock  spectra 
generated  from  a  typical  acceleration  data  record 
measured  near  the  base  of  the  air-handling  unit  is 
compared  to  a  fragility  curve  for  alr-handllng 
equipment  from  7MS>8SS'l  (Reference  2).  The  measured 
acceleration  levels  clearly  exceed  the  allowable 
limits  indtestad  by  the  fragility  curve.  Since  the 
slr-hendling  unit  continued  to  operate,  the  fragility 
curve  (at  least  in  this  case)  is  coniervative .  Peak 
in*structure  spectra  accelerations  'rom  1000  g's  for 
burled  shots  to  10,000  g's  for  aboveground  shots  were 
calculated  from  the  acceleration  data;  however, 
relative  displacements  associated  with  these 
accelerations  were  very  small  (0,01  Inch  to  0.001 
inch) .  Figure  2  shews  a  postest  view  of  pr-^- 
posltlonad  mannequins  in  sitting  and  standing 
positions.  Data  from  these  tests  indicate  the  only 
injury  to  personnel  would  have  been  from  falling 
objects,  such  as  the  bookshelf  in  Figure  2.  Data 
from  the  test  s*rles  conducted  in  1988  are  summarized 
below. 

TEST  DESCRIPTION 

A  series  of  three  experiments  using  buried 
general  purpose  bombs  were  conducted  on  Che  full 
scale  semihardened  test  structure  at  Tyndall  AFB,  FL. 
The  objective  of  these  experiments  was  to  detsrmlna 
the  shock  attsnuaclng  capabilities  of  four  simple 
shock  isolation  systems.  Several  pieces  of  equipment 
ware  placed  inside  the  building,  and  acceleration 
measurements  made  with  Che  equipment  both  hard- 
mounted  end  shock  isolated  using  various  shock 
isolation  methods.  The  equipment  varied  in  weight 
from  SO  lbs.  Co  2000  lbs.  The  equipment  layout  for 
Che  three  experiments  is  shown  In  Figure  3. 

The  general  purpose  bombs  viere  placed  in  Che 
same  ouried  location  for  each  test.  Each  item  of 
equipment  was  hard-mounted,  i.e.  bolted  directly  Co 
the  floor,  in  at  luart  one  test,  and  isolated  from 
Che  floor  in  various  configurations  for  the  remaining 
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tests  as  shown  In  Table  1.  Horizontal  and  vertical 
acceleration  was  measured  on  each  item  of  equipment 
and  on  the  floor  near  the  item  in  every  test.  Tne 
Isolation  systems  consisted  of  Cupmoiint  Series 
Isolators,  500  Series  Isolators,  and  neoprene  pads, 
as  shown  in  Figure  4.  The  neoprene  pads  were  in  some 
instances  simply  placed  beneath  the  equipment  while 
in  other  instances  the  pads  were  secured  to  the 
equipment  and  floor  slab  as  shown  in  Figure  5. 
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RESULTS 

Peak  values  of  acceleration  for  each  item  of 
equipment  for  the  various  mounting  configurations  are 
tabulated  in  Table  2.  Comparing  hard-mounted  and 
isolated  acceleration  values  from  Table  2  indicates 
that  the  Cupmount  Series  mounts  reduced  the  peak 
acceleration  in  the  X  (horizontal)  direction  by  about 
50  percent,  and  reduced  the  peak  acceleration  in  the 
Y  (vertical)  direction  by  about  30  percent.  The  500 
Series  mounts  Increased  the  average  peak  acceleration 
in  the  X  direction  by  about  20  percent  while 
decreasing  the  peak  acceleration  in  the  Y  direction 
by  about  15  percent.  The  neoprene  pads,  when  not 
attached  to  the  floor  slab,  decreased  the 
acceleration  in  the  X  direction  by  about  15  percent, 
but  increased  the  acceleration  in  the  Y  direction  by 
about  5  percent.  Also,  larger  displacements  were 
seen  with  the  unattached  pads  than  with  any  other 
isolators.  For  example.  Figure  6  shows  an 
approximate  11  inch  rigid  body  displacement  of  the 
1000-kw  generator  when  it  was  placed  on  two  layers 
(about  1/4  inch  thick)  of  neoprene  pads.  The 
neoprene  pads,  when  attached  to  the  equipment  and 
floor  as  shown  in  Figure  5,  performed  very  well. 
They  reduced  the  peak  acceleration  in  the  X  direction 
by  about  43  percent  and  reduced  the  peak  acceleration 
in  the  Y  direction  by  about  34  percent.  These 
numbers  were  very  similar  to  those  obtained  using  the 
cupmount  series  mounts.  All  of  the  data  from  the 
1988  tests  are  given  in  Reference  3. 

SUMMARY 

Relatively  high  peak  in-structure  acceleration 
levels  (in  excess  of  1,000  g's)  occurred  in  these 
tests.  however,  relative  displacements  were  very 
small  and  no  equipment  damage  or  personnel  injuries 
were  observed,  except  due  to  falling  objects.  More 
data  are  needed  on  equipment  fragility  in  this  very 
high- frequency,  highr^ shook  environment.  Rigid  body 
motions  of  unattached  equipment  and  furniture  can  be 
a  problem.  Positive  connections  to  the  floor  slab 
should  always  be  required,  tnd  caller  items  should 
be  prevented  from  overturning.  The  commonly  used 
Cupmount  Series  isolator  performed  very  well,  and  can 
be  expected  to  reduce  peak  acceleration  values  by 
about  50%  compared  to  a  hard-mounted  configuration. 
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Figure  1.  Shock  spectra  for  aboveground  detonation 


Figure  2.  Posteat  view  of  mannequins 
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1000-KW  GENERATOR 


Table  1 »  Shock  isolatioi  oonf iguraticns 


Test  1 

Test  2 

Test  3 

1,000  kW  generator 

Hard-mounted 

Cupmount  series 

Neoprene  pads 
two  layers 

Air-handling  unit 

Hard-mounted 

500  series 

Cupmount  series 

T1  Silent  700 

Mo.  1 

Neoprene  pads 

500  series 

Cupmount  series 

TI  Silent  700 

Mo.  2 

Hard-mounted 

Hard-mounted 

Hard-mounted 

Oscillosoope  rack 

Mo.  1 

500  series 

Cupmount  series 

Neoprene  pads 

Oscillosoope  rack 

Mo.  2 

Hard-mounted 

Hard-mounted 

Hard-mounted 

Tektronix  4081 

Mo.  1 

Neoprene  pads 

Hard-mounted 

Cupmount  series 

Tektronix  4081 

Mo.  2 

Neoprene  pads 

Cupnomt  series 

Hard-mounted 

Tektronix  4081 

No.  3 

Neoprene  pads 

Cupmount  series 

Hard-mounted 

Tektronix  4081 

Mo,  4 

Neoprene  pads 

Cupmount  series 

Hard-mounted 

Table  2.  Sumaxy  of  peak  acceleration  data  (in  g's) 


Item 

Hard-mounted 

500  series 

Mocsgcno-Bada 

1000  kH  Generator 

X-32.3 

X-16.1 

X-4.19 

Y-40.4 

Y-30.2 

- - 

Y-14.1 

Air-handling  vxiit 

X-9.81 

X-2.01 

X«2.35 

Y-13.2 

Y-3.16 

Y«2.82 

TI  Silent  700 

X-26.7* 

X-17.0 

X«21.1 

X-22.2 

Y-36.8* 

Y-24.3 

Y-24v3 

Y-79.6 

1  tack 

X-5.49* 

XX8.32 

X-17.2 

X-2.44** 

Y-31.8* 

Y-21.1 

X«65v6 

Y-14.9** 

Tektronix  4081 

X>9.53 

X-2.98 

— 

X-12.8 

No.  1 

Y-2S.9 

Y-14.1 

— 

Y-30.6 

'Mctronix  4081 

X-7.60 

X-4.98 

— 

X-43.0 

No.  2 

Y-24.1 

Y-15.8 

— — 

Y-46.9 

Tektronix  4081 

X-9.87 

X-3.0S 

X-7.86 

No.  3 

Y-24.5 

Y-17.7 

— — 

Y-32.6 

Tektronix  4081 

X-5.52 

X-7.48 

X-7.77 

No.  4 

Y-22.6 

Y**21 .8 

— 

Y-7.49 

Denotee  average  of  several  readings. 

Neoiixene  pads  were  attadwd  to  both  the  floor  and  the  eguipnent  rack  as  shown  in 
Figure  3. 

This  type  of  mount  not  used  for  the  item. 
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A  series  of  model  experiments  have  been  conducted 
to  investigate  the  dynamic  response  of  pile 
foundations  in  a  blast  and  shock  environment.  The 
model  experiments  were  conducted  at  the  Cambridge 
Geotechnical  Centrifuge  facility,  Cambridge  UK. 
Two  experiments  were  carried  out  on  a  pair  of 
single  piles  and  one  on  a  group  of  six  piles  at 
60g,  and  a  fourth  experiment  on  a  pair  of  single 
piles  was  conducted  at  Ig  for  comparison.  The 
model  piles  were  hollow  aluminium  alloy  tubes, 
Instrusiented  with  strain  gauges  in  the  horizontal 
and  vertical  planes  to  measure  bending  and  axial 
strain.  The  test  bed  was  a  fine  grained  sand, 
saturated  with  water.  The  blast  load  on  the 
foundation  was  provided  by  a  2  gm  charge  of 
pentaerythritol  tetranltrate  (PBTN)  placed  below 
the  surface  at  a  depth  of  one-half  of  the  pile 
length  and  detonated  at  a  distance  of 
approximately  one  crater  radius.  The  paper 
describes  the  experimental  techniques  employed  and 
presents  results  from  the  model  tests.  The 
results  highlight  the  differences  in  pile 
response,  particularly  with  depth,  between  the  Ig 
and  the  80g  experiments  and  confirm  the  importance 
of  correctly  scaling  geostatlc  stresses. 


1.0  INTRODUCTION 

Evaluation  of  the  response  of  foundations  to  blast 
loads  is  complicated  by  the  fact  that  soil 
properties  are  sensitive  functions  of  overburden 
pressure.  in  order  to  properly  reproduce  this 
effect  in  a  model  one  must  either  use  a  full  scale 
model  or  increase  the  rate  at  which  the  overburden 
stress  increases  with  depth.  Field  testing  on  a 
full  scale  prototype  pile  foundation  is  diffioult, 
generally  expensive  and  in  some  cases  not 
feasible.  Therefore,  a  database  from  which  to 
develop  a  rational  design  method  cannot  readily  be 
developed  using  field  testing  alone.  The  series 
of  centrifuge  experiments  repotted  in  this  paper 
were  undertaken  to  determine  the  feasibility  of 
using  subscale  models  on  a  centrifuge  to  collect 
valid  and  reliable  data  on  the  blast  response  of 
pile  foundations.  Discussion  of  the  use  of  the 
geotechnical  centrifuge  for  dynamic  modelling  and 
on  the  necessary  scaling  relations  for  Ig  and  for 
centrifuge  modelling  may  be  found  in  the 
literature,  Schofield  (1981),  Schofield  and 
Steadman  (1988) . 

The  objective  of  these  experiments  was  to  show 
that  the  centrifuge  is  an  appropriate  means  to 


gain  an  understanding  of  pile  foundation  response 
under  blast  loading.  The  design  parameters  were 
selected  to  be  representative  of  a  typical  hollow 
reinforced  concrete  pile  in  a  sand  foundation  but 
not  a  precise  model  of  a  specific  prototype.  This 
approach  was  selected  because  the  authors  believed 
that  there  was  more  to  be  learned  at  this  stage 
from  a  study  of  a  generic  problem  which  identified 
phenomena,  rather  than  from  models  which  attempted 
to  reproduce  details  of  a  sice  specific  problem. 
Hence  Che  results  are  more  readily  extended  to 
other  systems. 

The  paper  describes  the  design  and  instrumentation 
of  the  modal  piles,  and  aspects  of  model 
construction.  Results  are  presented  in  terms  of 
deformations  and  bending  momenta  as  a  function  of 
length  for  two  Isolated  single  piles  fixed  at 
ground  level.  The  cratering  and  ground  motion 
data  collected  are  presented  in  a  companion  paper, 
Gaffney  et  al.  (1989).  Results  from  the  initial 
experiments  were  reported  by  Felice  at  al.  (1988) . 


2.C  EXFERIMENTAI,  PROGRAMME 

2 . 1  Model  charge 

Four  experiments  were  conducted  on  single  piles 
and  a  line  of  piles  forming  a  pile  group  as 
described  above.  The  charge  was  selected  to 
simulate  the  detonation  of  a  1000  lb  conventional 
munition.  At  l/60th  scale  ,  2  gm  of  explosive  was 
required  (scaling  relations  ace  shown  in  Table  1) . 
This  was  accompished  by  packing  1.8  gm  PETN  inside 
a  plastic  sphere.  The  remaining  0.2  gm  was  made 
up  by  the  detonator.  The  average  density  of  the 
PETN  powder  was  0.9  gm/cm^.  The  detonator,  a 
Reynolds  RP-80  exploding  bridge  wire  was  placed  in 
the  sphoce  to  initiate  the  explosion  from  the 
centre.  The  charge  was  coated  with  a  thin  film  of 
epoxy  to  protect  it  against  moisture. 

2.2  Model  piles 

The  model  piles  were  constructed  from  8. 35  mm 
outer  diameter  aluminium  alloy  (ducal)  tube  with 
an  inner  diameter  of  5.0  mm.  Prior  to 
instrumenting,  the  piles  wore  turned  down  to  an 
outer  diameter  of  5.8  mm.  Each  pile  extended  147 
mm  below  its  top  fixity.  The  elastic  Young's 
modulus  and  density  of  dural  are  69  x  10^  MPa  and 
2.83  gm/cm^,  respectively.  The  manufacturer's 
specification  set  the  elastic  limit  as  255  MPa  and 
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th*  mlnimtun  taaslle  strength  as  310  MPa.  These 
data  give  a  model  bonding  stiffness  of  El  -  1.72 
Nm^  and  a  fully  plastic  moment  capacity  of  Mp  - 
3.625  Itm.  From  the  scaling  relations,  the 
equivalent  prototype  pile  is  seen  to  be  8.S2  m  in 
length,  with  a  bending  stiffness  of  1.72  x  60^  • 
22.3  MKm^,  and  a  plastic  moment  capacity  of  3.625 
X  60^  •  0.783  KNm.  The  head  of  each  pile  was 
fixed  into  a  steel  gantry  at  a  range  of  either  125 
mm  (7.5  m  prototype)  or  158  mm  (9.5  m  prototype), 
as  shown  in  Pig.  1. 

Selected  piles  and  the  ground  beam  were 
Instrumented  with  full  bridge  strain  gauge 
circuits  to  record  bonding  or  axial  strain.  To 
protect  the  gauges  and  wiring  from  moisture,  they 
were  coated  with  a  polyurethane  varnish  and 
covered  with  a  heat  shrin)c  plastic  tubing, 


2.3  Model  construction  and  layout 

Pig.  1  shows  the  containment  system  used  for  the 
latter  two  of  the  four  experiments.  Two  circular 
tubs  were  used,  one  sitting  on  a  rubber  mat  inside 
the  other,  with  an  air  gap  of  approximately  24  mm 
separating  them  around  the  perimeter.  This  was  a 
development  from  the  containment  used  for  the 
earlier  experiments  (RSS.130  and  131,  which  used 
only  a  single  tub)  following  concerns  raised  over 
the  level  of  safoty  that  a  single  tub  provided. 

Breexe  bloe)(  (a  porous  concrete  patio  bloc)c)  was 
placed  at  the  base  of  the  model  to  simulate  an 
underlying  bedcoc)c.  The  model  was  then 
constructed  by  pouring  a  uniform  dry  sand  layer  to 
a  depth  of  ISO  mm  (9  m  at  prototype  scale) .  The 
aand  used  in  the  model  was  a  Leighton  Buxxard 
100/200  sand  with  a  nominal  grain  sise  of  0.12  nrn 
and  specific  gravity  of  2.65.  100/200  denote  the 
British  Standard  sieve  sixes  through  which  the 
sand  should  pass/b«  retained. 

Prior  to  the  sand-pouring  the  piles,  which  had 
been  clamped  at  the  top  into  a  heavy  steel  gantry, 
were  fixed  in  position  by  locating  the  gantry 
onto  the  model  chamber.  In  each  case  the  gantry 
spanned  across  the  tub,  and  was  securely  bolted  to 
the  stiff  rim  of  the  outer  tub.  For  the  piir  of 
piles  experiments,  one  pile  was  positioned  to  be 
just  inside  the  crater  with  the  second  near  the 
crater  lip,  about  90*  further  round  the  crater 
perimeter.  A  small  clearance  existed  between  the 
piles  and  the  breexe  bloc)c  beneath  to  avoid 
unpredictable  axial  bearing  forces  in  the  piles. 

Sand  was  rained  from  a  hopper  suspended  above  the 
tub  using  a  constant  height  of  drop.  Pouring  was 
interrupted  to  allow  the  placing  of  transducers  in 
the  free  field  (see  companion  paper  for  a 
description  and  discussion  of  the  free  field 
instruntentatlon)  .  The  sand  was  levelled  and  the 
model  saturated  by  sealing  the  chamber  with  a 
heavy  lid  and  drawing  in .  a  calculated  volume  of 
water  under  vacuum.  in  the  second  series  of 
experiments,  RSS.140  and  141,  CO2  was  flushed 
through  the  chamber  pricr  to  the  introduction  of 
water . 


For  all  four  experiments,  the  charge  was  placed  at 
a  depth  of  74  mm  (half  the  pile  length) 
immediately  prior  to  mounting  the  pac)cage  on  the 
centrifuge  swinging  platform.  A  short  thin  walled 
brass  tube  which  had  been  placed  in  the  sand  bed 
during  sand  pouring  to  mar)c  the  location  of  the 
charge  was  excavated  to  the  correct  depth  where 
the  charge  was  placed  at  the  bottom  of  the  hole. 
The  brass  tube  was  then  bac)cfllled  with  the 
excavated  sand  and  then  extracted  by  gently 
vibrating  the  tube  and  pulling  it  upwards.  The 
water  level  inside  the  tube  was  maintained  at  a 
constant  level  during  the  process. 

The  completed  modal  was  then  mounted  on  the 
centrifuge  and  accelerated  to  60  g.  The  FS-10 
firing  control  unit,  which  had  also  been  mounted 
on  the  centrifuge,  was  triggered  remote.'.y  from  the 
control  room. 

The  procedure  for  the  1  g  experiment  was  identical 
in  all  respects  except  that  it  was  detonated  on 
the  laboratory  floor  instead  of  on  the  centrifuge. 
Care  was  ta)cen  to  ensure  that  the  model  was  level 
and  that  the  ground  water  level  was  exactly  at  the 
ground  surface. 

Table  2  summarises  the  model  parameters  for  each 
of  the  model  tests. 


3.0  RESULTS  AMD  DISCUSSIOM 
3 . 1  Bending  moments 

All  results  are  plotted  in  terms  of  prototype 
dimensions.  Fig.  2  shows  time  histories  of  the 
development  of  bending  moment  from  a  typical 
isolated  single  pile  in  a  60g  experiment.  A  rapid 
build-up  of  bending  moment  is  followed  by  a  slower 
decay,  with  the  duration  of  plastic  straining 
being  about  0.6  seconds  (10  msec  'real  time'  in 
the  model)  .  The  records,  which  are  from  strain 
gauge  bridges  at  different  depths  on  the  same 
pile,  show  strongly  consistent  data,  and  this 
enables  detailed  consideration  to  be  given  to  the 
profile  of  bending  moment  with  time  along  the 
pile. 

In  contrast.  Fig.  3  shows  the  equivalent  data  from 
a  Ig  experiment.  There  is  a  longer  period  at  or 
near  the  pea)(  bending  strain  at  each  depth,  and  a 
less  rapid  build-up  to  the  pea)c  strain.  The 
magnitudes  of  pea)c  strain  are  comparable  between 
the  Ig  and  the  60g  model  as  would  be  expected 
since  in  both  cases  pea)(  strains  are  limited  by 
the  plastic  moment  capacity  of  the  pile. 

It  is  clear  from  the  plots  of  displacement  profile 
shown  below  that  all  piles  developed  plastic 
hinges  and  failed,  with  large  lateral  permanent 
displacements.  The  pea)c  bending  strains  recorded 
by  the  strain  gauges  are  large  and  in  excess  of 
the  manufacturer's  stated  elastic  range.  On  a  few 
CQcasions  the  local  bending  strains  were  sc  large 
(2-3%)  that  the  gauges  brojre  down  and  the  signal 
from  that  circuit  was  lost.  Clearly  as  strains  in 
the  outer  fibres  of  the  dural  piles  exceed  a  yield 
point  the  relationship  between  moment  and  measured 
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bending  strain  is  no  longer  linear.  Thus  although 
peak  moments  are  of  great  Interest,  the 
development  of  moment  as  a  function  of  time  and 
depth  may  provide  more  significant  clues  to  the 
ceponse  of  the  pile  from  an  analytical  view-point. 
Furthermore,  In  the  design  of  the  model  piles 
correct  scaling  of  the  bending  stiffness  was 
chosen  at  the  expense  of  the  correct  yield 
characteristics . 

Fig.  4  con^ares  the  development  of  bending  moment 
with  time  from  an  Isolated  single  pile  In  a  60g 
experiment  to  the  development  of  moment  In  a  Ig 
experiment . 

Time  Intervals  have  been  chosen  to  be  shortly 
after  the  first  significant  bulld-up  of  strain  and 
then  at  equal  time  Intervals  of  46.9  msec  until 
the  peak  strain  was  reached,  or  shortly 
thereafter.  A  'best  fit'  fourth  order  polynomial 
has  been  drawn  through  the  peak  data  of  both  Ig 
and  60g  experiments. 

One  clear  feature  of  the  data  Is  the  marked 
similarity  between  the  strain  distributions  at  the 
early  (tl)  stage,  In  contrast  to  the  later 
bulld-up  towards  a  peak.  At  the  t2  or  t3  stage 
the  peak  moment  Is  clearly  at  a  shallower  depth  In 
the  60g  experiment  than  In  the  ig  experiment 
<whlch  consistently  shows  a  peak  near  shot  depth) . 
The  60g  model,  however,  has  a  distribution  biased 
upwards,  towards  the  ground  surface. 

This  is  likely  to  be  due  to  the  increased 
stiffness  of  the  sand  bed  In  the  60g  experiment  In 
comparison  to  the  Ig  model  (a  factor  of  8  under 
static  lead  conditions) .  However,  this  factor  may 
have  been  considerably  reduced  around  the  time  of 
the  maximum  moments  <t2  or  t3  In  rig.  4)  because 
of  the  large  negative  pore  pressures  In  the  Ig 
experiment  which  followed  the  peak  pore  pressure 
'wave*.  In  the  60g  experiments,  the  pore  pressure 
wave  lasted  longer  and  decayed  considerably  more 
slowly. 

The  early  time  data  of  bending  strains  Is  enlarged 
In  Fig.  S  together  with  a  polynomial  which  follows 
the  general  trend  of  the  data.  At  this  time  In 
both  the  Ig  and  the  60g  experiments,  the  pore 
pressure  uround  the  pile  was  approaching  a  peak, 
with  a  corresponding  reduction  In  stiffness. 

The  polynomial  In  Fig.  S  was  used  to  deduce  the 
general  trend  of  the  pressure  distribution,  by 
differentiation.  A  check  was  made,  by  Integrating 
the  polynomial,  that  the  slope  at  the  pile  top  and 
the  displacement  at  the  pile  bottom  were  minimised 
or  zero.  It  Is  clear  that  the  trend  of  the 
pressure  distribution  at  this  time  la  as  shown  In 
Fig.  6,  In  which  a  high  load  at  the  p'lo  cop 
decays  to  near  zero  around  the  central  portion  of 
Che  pile  together  with  some  restraint  near  the 
pile  bottom. 

Preliminary  conclusions  are  that  the  central 
portion  of  the  pile  moves  with  the  liquefied  soli, 
without  load  or  resistance.  Near  the  top  of  the 
pile,  the  soil  la  being  driven  past  the  stationary 
pile  head,  Imposing  very  substantial  loading  and 
this  Is  the  most  likely  location  of  the  first 


plastic  hinge. 

3.2  Displacements 

Fig.  7  shows  the  profile  of  displacement  with 
depth  for  piles  (1)  and  (2)  (range  7.Sra  and  9.Sm) 
In  both  a  Ig  (141)  and  60g  (140)  experiment. 
Measurements  were  made  using  an  image  analyser  at 
AFWL. 

The  initial  slopes  of  the  piles  are  very 
consistent,  but  the  Increased  restraint  caused  by 
the  higher  overburden  stresses  at  the  base  of  the 
60g  model  limits  Che  outward  movement  of  these 
piles.  Clearly  these  displacements  are  much  more 
closely  linked  to  the  distribution  of  soil  strain 
in  the  sand  beds,  which  Is  discussed  In  more 
detail  by  Gaffney  et  al.  (1989). 

The  depth  of  the  central  plastic  hinge  is  less  in 
Che  60g  than  In  the  Ig  experiments,  which  Is  in 
agreement  with  the  observations  of  bending  strain 
in  Fig.  4. 

A  critical  element  In  the  prediction  of  the  range 
of  damage  to  piles  in  the  field  is  Che  nature  and 
magnitude  of  the  pore  pressure  wave.  The  fast 
decay  of  the  pore  pressure  wave  In  the  Ig  model  Is 
likely  to  bo  due  to  the  strong  dilation  front 
which  followed  Immediately  behind.  Although  at 
high  g  a  dilation  front  will  recover  the  full 
field  (or  prototype)  strength  In  the  soil,  at  Ig 
the  strength  of  the  sand  Is  small  and  large 
strains  can  take  place. 

The  high  g  model  can  correctly  scale  the  pattern 
of  ground  strain  and  the  character  of  a  pore 
pressure  wave,  whereas  the  response  of  a  pile  In  a 
Ig  experiment,  close  In  to  the  charge,  will  be 
determined  by  the  strengt))  of  the  pile  and  the 
proximity  of  the  charge. 


4.0  Conclusions 

4.1  The  time  histories  of  bending  strain  In  model 
tests  have  provided  valuable  data  towards  the 
development  of  analytical  techniques.  In 
particular,  It  Is  clear  that  a  data-base  of 
soil-structure  Interaction  under  blast  loading  can 
be  developed  using  centrifuge  modelling. 

4.2  In  particular,  the  bulld-up  towards  a  peak 
strain  Is  Initiated  as  the  wave  of  pore  pressure 
arrives  at  the  pile,  softening  the  soil  around  tha 
pile.  High  loading  Is  observed  near  tha  pile  top 
with  a  minimum  of  load  around  the  central  section. 
The  peak  strain  Is  reached  as  the  pore  pressure 
wave  decays,  both  radially  and  with  depth.  Larger 
bending  strains  are  then  Invoked  as  the  deeper 
soil  consolidates  more  rapidly  and  provides 
restraint  against  lateral  movement. 
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ELEVATION  THROUGH  TUB.  Y-Z  PLANE 


Table  1 


Quantity  Ratio  of  model  to  prototype 


Length 

Velocity 

Acceleration 

Force 

Stress 

Energy 

Frequency 

Tine 


1/n 

I 

n 

l/n^ 

1 

1/n^ 

n 

1/n  tor  inertial  events 
Ittfi  for  diffusion 


’2- ELEVATION  THROUGH  TUB.  X-Z  PLANE 


2  nun  mbtef  mat 


Table  2 


sand 

sand 

void 

relative 

Exp't 

piles 

mass 

volume 

ratio  density 

V 

kg 

litres 

RSS.iSO 

2 

124. S 

84.6 

0.8 

604 

60 

RSS.131 

6  (line) 

126.5 

86.4 

0.62 

554 

60 

RSS.140 

2 

109.5 

■71.47 

0.74 

754 

60 

RSS.141 

2 

111.0 

71.97 

0.72 

824 

1 

Fig. I  Plan  and  cross- 
sections  through 
the  60g  centrifuge 
model  test  RSS.140 


Saturated 
Sand  bed 
Breeze  block 
Rubber  mat 
Innor  llnar 
Rubbar  mat 
SSOmm  tub 


PLAN  VIEW  OF  MODEL  RSS.140.  LOOKING  ALONG  Z  AXIS 
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